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HAUPTAUFSÄTZE 


Ausbreitung der Kraft in der Halbebene und im Halbraum 
bei anisotropem Material. 
Von K. Wolf in Wien. 


ie Untersuchungen, die in der Erdbaumechanik über die Verteilung der Spannungen und 

Verschiebungen im Erdboden bei der Fundierung von Bauwerken angestellt wurden, 
haben ziemlich beträchtliche Abweichungen gegenüber dem in einem elastischen Körper unter 
gleichen Bedingungen herrschenden Spannungszustand ergeben. Dies rührt wohl in erster 
Linie davon her, daß ein solcher Erdkörper nur sehr angenähert die Bedingungen eines 
elastischen Mediums erfüllt. Zum Teil mag aber dabei auch der Umstand eine Rolle spielen, 
daß der Erdboden sieh nieht isotrop verhält, daß etwa bei einem Boden, wo sandige und 
tonige Schichten aufeinanderfolgen, die Elastizitätsmodule in horizontaler und vertikaler 
Richtung stark voneinander verschieden sind. Daher soll für derartiges anisotropes Material 
im folgenden die Ausbreitung der Kraft, also das zum Boussinesqschen analoge Problem, im 
ebenen und im achsialsymmetrischen Fall dargestellt werden. Das letztere ist allerdings schon 
von Michell') behandelt worden; da aber die Art der Darstellung und die Methoden der 
Lösung etwas andere sind, so dürfte auch die Behandlung dieses Falles nieht überflüssig sein. 


Das ebene Problem: Wir wollen unser Koordinatensystem so legen, daß die positive 
y-Achse vertikal nach abwärts, die positive w-Achse horizontal nach links gerichtet 
(vel. Abb. 1). In der «-Richtung sei der Zugelastizitätsmodul #E,, in der Richtung E,: das 
Hookesche Gesetz können wir dann in der Form ansetzen 


0x Gy Gy 0x 
E, E, 
WO #x, &, und 0%, 6, die übliche Bedeutung haben — es sind die Dehnungs- und die Normal- 
spannungskomponenten —- », und »v, entsprechen den Poissonschen Zahlen für beide 


Riehtungen. Da wir mit einer Vermehrung der Zahl der Elastizitätskonstanten um eine 
einzige, also im ganzen mit drei auskommen wollen, so läge es nahe, v, =», zu setzen, da 
kleine Abweichungen in der Poissonschen Zahl keinen großen Einfluß auf das Spannungs- 
und Verschiebungsbild haben. Das geht aber nicht so ohne weiteres, da immer die Bedingung 
erfüllt sein muß, daß die Spannungskomponenten sich als Ableitungen der Verzerrungs- 
energiefunktion nach den Verzerrungskomponenten ergeben. Die elastischen Konstanten sind 
bekanntlich die Koeffizienten in dieser homogenen, quadratischen Form. Unter Berück- 


1) Proc. Lond. Math. Soc., Bd. 32, 1901, p. 247. 
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siehtigung der Symmetriebedingungen in unserem Ansatz — für den Übergang von +x zu 
x darf sich die Funktion nicht ändern — bekommen wir 


und daher 


Wenn man die Gl. (1) nach den Spannungskomponenten auflöst, so müssen demnach die 
Koeffizienten von &, und #, in den beiden Gleichungen übereinstimmen. Das gibt 
a 

nie E, E, 

wenn man an Stelle von », jetzt v» schreibt. Für den Schubelastizitätsmodul setzen wir 

E,+ E,(1-+2»)’ 
nenten die Gestalt haben 


(i so daß die Beziehungen zwischen Verzerrungs- und Spannungskompo- 


Die kubische Dilatation e hat die Form 
1—rv | | 3) 
Wenn die Airysche Spannungsfunktion F in der bekannten Weise eingeführt wird und 
die Spannungskomponenten daher die Form haben 


e 


so erhalten wir für F aus der Kompatibilitätsbedingung 
die Differentialgleichung 


wo I? F den Laplaceschen Operator bedeutet. In Verbindung mit (2) ergibt sich daraus 
für die Volumsdehnung die Beziehung F’e=0, e ist also auch jetzt eine harmonische Funk- 


. . . . . 
tion. Setzen wir 5 k?, so läßt sieh (4) in der Form schreiben 
0” 
(4a) 
x 
(xy) 


Abb. 1. 
Wir wollen zunächst den Fall betrachten, daß über ein Stück der Halbebene von der 
Breite 2a eine konstante Belastung p, verteilt ist, der übrige Rand sei spannungsfrei (siehe 
Abb. 1). Für isotropes Material hat die entsprechende Spannungsfunktion die Gestalt?) 


F: 


Po 


°) Love: Lehrbuch der Elastizität, $ 151. 


| 
| 
/ | 
| 
/ & 
| 
Y 
4 
Bit 


f.angew. 
und Mech. 


zu 


ich die 


en wir 


kompo- 


(4), 


araus 


"unk- 


(4a). 


der 
iehe 


(5). 


Band 15, Hast ” Wolf, Ausbreitung der Kraft in der Halbebene bei anisotropem Material 351 


Die Bezeichnungen sind dieselben wie in der Abbildung. Für den äolotropen Fall wollen 
wir zunächst die dem Wert r?d allein entsprechende Lösung suchen. 


Setzen wir in (4a) für y den Wert ky,, so erhalten wir 


0” 


da? 


Da V?(r? 9) 40 ist, so werden wir als Lösung von (6) den Wert 


nehmen. Führt man Polarkoordinaten ein 


®F 19F 1®F 


2 W- u - 
or? r roöor 


und setzt als Lösung F=r?© ein, wo © allein Funktion von 9 ist, so erhält man als 
Differentialgleichung für © 


Wendet man die Methode der Variation der Konstanten auf die Lösung der homogenen Gleichung 
an, wo also A und B als Funktionen von 9 angesehen und aus den Gleichungen 
A'’sin29 + b’cos2) = 0 2cos29 A’+2sin29 bD’ | 


bereehnet werden, so erhält man 


dA dB 
A 19 dd > ”sın?d (S) 


und daraus 


Beim Übergang zum isotropen Medium für k? 1 geht der Ausdruck für © in den richtigen 


Wert über, da 


ni 


lim 

und 


ist. Setzt man €, und €, gleich Null und führt wieder kartesische Koordinaten ein, so 
nimmt F die Gestalt an 


Dabei ist mit 


(il) 
und 
0.4 3,97 9A, gr 
oB oB 2 | 


. (6). 
. (8). 
d und 
_ 
| 
_ 
| 
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Subtrahiert man nun entsprechend der Formel (5) von der Lösung mit = -+a, y=0 
als Zentrum die mit dem Zentrum #== a, 9-0 (siehe die Abb. 1), so erhält man als end- 
gültigen Wert 


F aA, +le . (13), 


wo für A,B,, A,B, die Ausdrücke aus den Gl. (11) eingesetzt sind, das eine Mal mit dem 
Weiser I für 9, 9’ und r, r’ und mit = — a statt », das andere Mal mit dem Weiser 2 und 
statt 


| ‚A+k 
Die Spannungskomponenten ergeben sieh daraus, wenn man noch e = Pol r ) setzt, zu 

1 . . . . . . . . 

Po 


Für y 0, alor, = rn, =x+a,r, kr, r,=kr, verschwindet und nimmt 
außerhalb von „= -+ a und außerhalb von .r a den Wert Null an, da dort 9, 9,9%, 
— W', gleich Null beziehungsweise gleich x ist. Auf der Strecke zwischen +a und — a haben 
wir aber o,=p,, da dort 9, 9’, = -+-nz und 9,—=%,--0 ist. Man erkennt also, daß die 
Randbedingung befriedigt ist. Wenn man auch das Eigengewicht des Materiales berück- 
siehtigt, so hat man in bekannter Weise zu o, noch den Summand — y.= hinzuzufügen; y ist 
dabei das spezifische Gewicht des Erdreiches. 

Um nun die Lösung für eine Einzelkraft zu erhalten, hat man nur (13) nach a zu differen- 
tiieren und zur Grenze a gleieh Null überzugehen. Es ist 


2yA,- 2(2-+ a) B, — — a)?’ Sal 
und 


wobei 9, A,=A,=A und B=B,=B gesetzt ist. 
Bezeiechnet man den Betrag der im Ursprung vertikal nach abwärts wirkenden Kraft 


> 


mit 7%, so ist e - zu setzen und man erhält durch Differentiation unter Berücksichtigung 


w 


der Gl. (12) die Spannungskomponenten 


2Pkay 2P kyf 2Pkxy? 


(16). 
Fürk=1,also r’  r, gehen diese Ausdrücke in die bekannten Werte der Boussinesqschen 
Lösung für die Halbebene über. 
Die Werte der Verschiebungen & und », ergeben sich aus den Gl. (2) zu 
2P 1 
.T E, ( E 
TER 


‘2 


3 


Für E,— E, bekommt man in Übereinstimmung mit der Boussinesqschen Lösung 
Für die Oberfläche y„=0 ergibt sich für die Einsenkung », der Wert 


2Pp 


Inr. 


| 
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Der axialsymmetrische Fall. Bei Verwendung von Zylinderkoordinaten 0,9 und z sind 
‚lie Verzerrungskomponenten bei Berücksichtigung der axialen Symmetrie 
u ou dm 


Yz S), 


wobei « die Verschiebung in radialer und w die Verschiebung in axialer Richtung bedeutet. 
Die spezifische Verzerrungsenergiefunktion hat die Form 


‘ 
die Glieder &0 y:o... verschwinden, weil für +z und — z der Wert der Funktion sieh nicht 
ändern darf. Berücksiehtigt man jetzt die Querisotropie in der Weise, daß man ähnlich wie 
im vorher behandelten Fall das Hookesche Gesetz in der Form anschreibt: 


00 OÖ; 00 Öz 00 


E, ist dabei der Elastizitätsmodul in der zur Achse senkrechten Richtung, FE, der in der 
Richtung der Achse, », und », sind die analogen Poissonschen Zahlen. Will man wieder 
die Zahl der elastischen Konstanten nur um eine vermehren, so muß man, ebenso wie früher, 
die Bedingung berücksichtigen, daß sich die Spannungskomponenten als Ableitungen der Ver- 
zerrungsenergiefunktion nach den Verzerrungskomponenten darstellen; man erhält daraus die 
Beziehung v, E,—=r,E,. Da nun bei Erdkörpern die Werte von vr sehr klein sind und ihre 
experimentelle Bestimmung unzuverlässig ist, möge für das Weitere v, und »v, gleich Null 
gesetzt werden, obwohl sich die Rechnung auch allgemein durchführen ließe. Dann bekommen 
wir aus (19) 

E, E, 


E, en, 0, = E, E, To: E, E, ’o: (20). 
Setzt man dies unter Berücksichtigung von (18) in die Gleichgewichtsbedingungen 
Ö 
00 0 T0:) 0,0 0 To; AR 


ein, so erkennt man durch Differentiation der beiden Gleichungen nach o und z und darauf- 
folgende Addition, daß die kubische Dilatation e=- ‚.eW+,, wieder der Gleichung 
000” 02 
genügt. 


Drückt man nun analog dem bekannten Verfahren bei Isotropie’) « und »w= durch eine 
Funktion P in folgender Weise aus: 


so nimmt e die Gestalt an 
| 
o do E,02°)' 
Aus Gl. (22) ergibt sich dann für ® die Differentialgleichung 
ı107/ 
)+ 
oder wenn man für E wieder A? setzt 
10 Ö > 


Bei dem entsprechenden Problem für den isotropen Fall ist P, abgesehen von einer 
multiplikativen Konstanten durch den Ansatz 


P= In(r-+2) 


3) hLove-Timpe: Lehrbuch der Elastizität, 188. 
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1 
gegeben — r ist dabei yo?’+z° — und daher Y’P = En Wenn man in (24) an Stelle von z 


den Wert 2,=kz einführt, so nimmt der Klammerausdruck, geschrieben in o und z,, die 
Form des Laplaceschen Operators an und es ist demnach 


eine Lösung dieser Gleichung, wobei r’=yo?+ k’z’ gesetzt ist. Es ist also nur mehr die 
(41. (25) zu integrieren. In räumlichen Polarkoordinaten hat bei Berücksichtigung der axialen 
Symmetrie (25) die Gestalt 

Or’ ror 00 ed), 
+ bedeutet dabei, wie üblich, den Winkel zwischen r und der positiven z-Achse. Setzt man 
PD gleich r ©, wo © eine Funktion von W allein sein soll, so erhält man 

do l 

„tr eotg d +209= 
du" reotg + cos?) 
Das Integral der homogenen Gleichung ist bekanntlich durch die Kugelfunktion erster Ord- | 
nung gegeben*) und mit der Methode der Variation der Konstanten ergibt sich das allgemeine 
Integral von (27) zu 


os», 1+cos) 


2 
wo 


gesetzt sind, wobei und 9 durch die Beziehung cos 1 eos verknüpft sind. Führt 


‚_yVi—kz 
man für und wieder z ein, gemäß den Formeln cos — ‚so be- 
kommt man 
l r-+2z +kz 
P—= 
2 r —kz 


Um für 2=1 den Wert für den isotropen Körper zu erhalten, hat man noch den Ausdruck 


2 


In(»® der ja ebenfalls eine Lösung von ist, hinzuzufügen und alle 


Yı 4 
Glieder mit #° zu multiplizieren. Dann ergibt sich als endgültiger Wert: 


Aus (23) erhält man mit 0° r? 2°’, wenn die im Zentrum nach abwärts wirkende Kraft 
mit 2’ bezeichnet wird: 


PE, | 2/E, E, | PE, | | | E, 1 | 
PE,]1 
Auf der Oberfläche ıst Die Spannungskomponenten sind: 
PE‚, z L/E, E, | £, 
E, 2 E, | (30 
PE, | E PE, | | | 
Beim Grenzübergang für k==1, also FE, = E,— E bekommt man die bekannten, für 
Isotropie gültigen Werte. 493 
4) Jahnke u. Emde: Funktionentafeln, 2. Aufl., S. 176. 
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Numerische Durchrechnung zu der Abhandlung 
über Grundwasserströmung'). 


(Abhandlungen zur Hydrodynamik IV.) 
Von Georg Hamel und Erich Günther in Berlin. 


er zweite Verfasser hat sich der mühsamen Arbeit unterzogen, zu der Abhandlung «des 

D ersten Verfassers die nötigen numerischen Rechnungen auszuführen. Wir veröffentlichen 
hier von diesen Rechnungen einen Auszug, und zwar geben wir die Methoden und Zahlen- 
werte zur Berechnung des inneren Integrals 

ausführlich, für das Intervall in 11,2, für die anderen Intervalle in VI,1.2.3. 
Dieses Integral wird immer wieder gebraucht werden, wenn jemand weitere Fälle ausrechnen 
will. Dagegen geben wir die Schlußreehnung, das äußere Integral nur für je einen Wert der 
beiden Parameter, stark gekürzt. Wir geben Interessenten natürlich gerne unser volles Zahlen- 
material. Das Ergebnis ist in VII in der Schlußfigur dargestellt. 

Die Gelegenheit soll benutzt werden, hier eine wesentliche Beriehtigung zu der 
Arbeit von 1934 anzubringen. 4. Die Randwertaufgabe, Fall2, Wassergrenze ent- 
hält einen Zeichenfehler. In der letzten Zeile von Seite 135 muß es links statt rechts heißen. 
Infolgedessen muß es in Zeile 4 von Seite 136 heißen: 


., dw 
In den folgenden Zeilen ist er >0 zweimal zu streichen, in Zeile 7 muß grade in un- 


grade verbessert werden und hinter also hinzugefügt werden: da 


In Zeile 10 muß umgekehrt ungrade in grade verbessert werden und hinzugefügt werden: da 


A 


Die unterstrichenen Hauptergebnisse sind richtig, und deshalb ändert sich an den weiteren 
Ergebnissen nichts. Hamel. 


I. Formelzusammenstellung. 
Das Strömungsbild des Grundwassers ist bestimmt durch das Integral 


Dabei gilt 


> 
g" wog die Variable in der Ebene der Hodographenfigur ist (Abb. 1). » be- 
4 deutet eine positive Materialkonstante, 9 die Anziehungsbeschleunigung 
der Erde. 
+4” Durch die Modulfunktion /(g) wird die q-Ebene auf die obere Halb- 
g-Ebene ebene abgebildet, und zwar so, daß 


i(0)=0, il)=1, Al+ie)=o 
wird (Abb. 2). Der 4-Wert von A4’”” wird a, der von E’”’ wird b genannt, 


wobei a>b. Die für einen bestimmten Strömungsfall konstanten a und b 
0 7 
D” pm gm >b" 
Abh. 1. Abb. 


1) In dieser Zeitschrift, Bd. 14, Heft >. 


| | 
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(27). 
| | | 
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wurden in der vorliegenden Arbeit willkürlich angenommen. y ist der Zentriwinkel eines 
Bogens des Halbkreises DC” (Abb. D). y und 4 sind miteinander verbunden durch 
Die Konstante © bestimmt das Größenverhältnis der Strömungsfigur, 
Nach Schwarz, Artikel 42, oder „Abhandlung“, Seite 141, gelten die Entwicklungen : 


in der Nähe von 


in der Nähe von /=1, q=1: 
g=1+ 


| h 
In, +1n16 4) (1 


in der Nähe von vw, gqg=w: 


‚‚„1,121 
II. Die Berechnung des inneren Integrals vi ai für 0.301. | 
1. Die Berechnung von y(M). Aus Formel (1) folgt 
> aretgil q 
Weeen (2) 
| | 3 -zi+1n16 [In + At +. 
Mithin ist 


Nun ist (Schwarz, Artikel 42) der aus der Theorie der elliptischen Funktionen bekannte 
Quotient 


Es wird also 
RK’ 


‘ 


‚k) und sind vollständige elliptische Integrale, die mit oder als 


Argument in der Tafel von Hayashı°), tabuliert sind. 
Bei «dem Erörterungen über das Poissonsche Integral (Abhandlung S. 139) ist f als 
Integrationsvariable in der Ebene gewählt worden. 4 selbst dagegen ist in dem inneren 


Integral \, ; dt ein Parameter. y(f) läßt sieh also punktweise leicht aus der genannten 


Tafel für jedes f berechnen, als Funktion von f dagegen wird sie beim Integrieren recht un- 
handlich. 
Die folgende Tabelle ist mit Hilfe der Tafel von Hayashi, Seite 61, berechnet worden. 


Havashı benutzt andere Bezeichnungen. Unser f entspricht bei ihm dem 4°, Ma K dem 


loxy. Die Beziehung zwischen und unserem 4 ist gegeben durch 4° 

Abb. > veranschaulieht den Verlauf der Funktion y=y(f). Die Kurve ist symmetrisch. 
Ihr Wendepunkt liegt bei #=0,5. Die Symmetrie der Kurve erklärt sich aus den Abbildungs- 
eirenschaften der Modulfunktion. Den Punkten des Kreises im Moduldreieck werden die 


-, KR. Hayashi: Tafeln der Besselschen — Theta Kugel- und andern Funktionen. Verl. Springer 1930, 


| 
| 
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K IN N } 
t K K arctg t M»n K K arctg 
(1). 0.00001 | 7 0.790588 | 4,547 1.354 0.013 091267, 2.268 1.155 
0.000026 — 0,09691 4,321 1.344 0,016 3 — 0,00350 | 2,196 1.143 
0.000083 6 0,27301 4.197 1.337 0,020 3— 0,10129 2,125 1.131 
uneen: 0.000046 0,39795 4.106 1.332 0,025 3 — 0.199350 | 2.053 1.117 
1.000056 0,48 4,035 1.328 | 
0.030 3.027960 1.994 1.106 
000006 6 — 0,57404 3.977 1.324 0,035 3 0,34766 | 1.044 1.096 
(2) 6 — 0.640099 3.028 1,321 0.040 -.0,40677 1.901 1.086 
6 — 069899 3,885 1.319 0,045 3 — 0.45905 1.862 1.078 
0.000096 — 0,67501 3.848 1.316 0,050 3.050595 1.829 1.070 
00010 0,67590 3.814 1.314 
0,06 3 — 058739 1.768 1.056 
(3), 0.0002 0,09695 3,594 1.299 0,07 3 0,6566 1.717 1.043 
1,0003 — 0.275307 3.165 1.289 3 0.716953 1.073 1.032 
0.0004 5.039805 3.9319 1.285 0,09 0,77042 1.734 1.022 
0.0005 0.49496 3.302 1.277 0,10 3081853 1.509 1.012 
(4). 0.0007 0,64113 3,195 1.267 012 3 -0,90249 1.537 0.004 
0.0009 0,75032 3.115 1.260 0.14 3 — 0,97429 1.485 
| 0,83751 3,051 1.254 0.16 2 — 0,03724 1,438 1.06: 
| 0.0013 091010 2.998 1,249 0,18 2 0, 09344 1.397 0,950 
0.0015 5 — 0,97230 2.952 1.244 0,20 2 — 0,14435 1.360 0.057 
0.002 0,09734 2,860 1.234 | 0,83 0,34786 1.211 0.880 
0.003 4 — 0,27365 2.7531 1.219 0,4 2 — 0,50356 1.097 0.832 
4 — 0,39881 2,639 1.209 05 3 — 0,63562 1.000 0.785 
0005104959 2,568 1,199 
2510 1,192 
0.007 4 0.64250 2,461 1.185 
4 .0,70071 2,418 1,179 
0,000 4 — 0,75208 2,380 1.173 
000 0,79806 2,347 1.169 
Kandwerte y zugeordnet. Bei der Abbildung des Dreiecksrandes auf die reelle Achse geht 
innte der Halbkreis €” D” in die Strecke €’"” D"’ der Länge I über (Abb. 4). Auf der Strecke 
bis 1=0,5 fällt y von bis . Nun gehen aber die Punkte des Bereiches II 
durch die Transformation 
aus den Punkten des Bereiches IV hervor. Daraus folgt die Symmetrie der Kurve y y(d): 
"als 
' als Er 
eren | 
nten 
(den. I 170 
dem 
| 
sch. 
die [45413] 0,5 1 p” 


| Abb. 3. Abb. 4. 
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2. Berechnung von \ ‚at für Die punktweise Berechnung von y(f) bietet 
nun zwei Wege für die Bereehnung des Integrals. Einmal kann zwecks graphischer Inte- 

eration der Integrand n ; für ein bestimmtes 4 punktweise berechnet werden. Diese Methode 
wird für und 1=7=% benutzt. 

Es läßt sich aber auch y = y(f) dureh eine bequemer zu handhabende Funktion ersetzen. 
Der komplizierte Aufbau von y = y(f) erfordert zur Approximation ganze rationale Funktionen 
mit vielen Potenzen, wenn die Annäherung in einem hinlänglich großen Intervall gelten soll. 
Innerhalb derselben Fehlergrenzen läßt sich im selben Intervall aber auch die Annäherung 
mit wenigen Sehnen erreichen. Dabei ist als rechnerischer Vorteil bei der späteren Integration 
noch zu beachten, daß außer In f nur f in der ersten Potenz auftritt. Im Intervall A=0 bis 
;—=1 wird die Kurve (Abb. 3) dureh ein Sehnenpolygon approximiert und zwar so, daß die 
Differenz zwischen Kurve und Sehne nie größer als 0,01 wird. Nur im ersten Intervall 
bis #,.0,00001 7, 2,709) wird der Approximationsfehler gleich y, — = 0,45 
gesetzt. Das verwendete Polygon besteht aus 12 Sehnen, die zu den Intervallen mit den 
Endpunkten 0,00001: 0,00005 ;.0,0001 5 0,0005; 0,0015; 0,005: 0,01: 0,025: 0,05; 0,1; 
0,2: 0,5 gehören. 

Die Gleichungen der Sehnen sind 


| 13293.277D + 0,4: 10° 
Il 1323,6421 -t+2,72159  +025-10° 
II 535.4427  +02-10 
Vy 65.082841 02:10 
VII y= 12.60035 1-- 246198 
y= 6,73467 :-+2,40332 
IN y 37542 +02-10°.4 
3 y: +02-10°.1 
15018 -+2,17394  +02-10 
All y= 1.0003 207545 +02-10 


Wegen der Symmetrie von y=y{M) braucht die Approximation zunächst nur für das 
Intervall #0 bis t=0,5 durchgeführt zu werden. Die Symmetrie der Kurve überträgt sich 
auch auf die approximierenden Sehnen, so daß in Anlehnung an (5) gilt: 

Unter Berücksichtigung dieser Beziehung läßt sich die gesamte Anzahl der zu integrierenden 
Gleichungen im Intervall #0 bis #= 1 um die Hälfte vermindern. Die folgende Tabelle (S. 260) 


| 1 
v(t) 
gibt als Ergebnis die Abhängigkeit der beiden Funktionen \; ; dt und e 27, !=*vondAan. 


3. Fehlerbereehnung. Die angegebenen möglichen Fehler enthalten die bei der numerischen 
Reehnung auftretenden Fehler neben den Fehlern, die durch die Sehnenapproximation entstehen. 
Ist der mögliche Fehler nicht besonders erwähnt, so kann die letzte Dezimale höchstens 
um 5» schwanken. 
4. Graphische Integration. Aus später zu erörternden Gründen (S. 259) wird das innere 
Integral für 4=0,5: 0,8; 0,9 auch noch auf graphischem Wege bestimmt. Wir zerlegen nämlich 


(yıd)—y(A) \ (A) 
„dt \ 
Aus der 
Durch 
folgt durch  Aus- 
zählung der Fläche: 5: 
) 
; 
0 
05 1.010 + 0,012 0,000 1,040 + 0,012 
0.8 1.178 0,003 1.758 — 2.936 40,003 
1,363 + 0,001 2,457 - 3,820 +0,001. 


Der angegebene Fehler gestattet in der Zeichnung eine Schwankung der Kurve um '/, mm. 
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III. Berechnung der freien Oberfläche. 
Nach der „Abhandlung“ (S. 144) ist 


7 7 

I 


Tor 
sind) 


Wir lösen die beiden Integrale durch graphische Integration. Um 


1 
b)la— 


zu bilden, wählen wir willkürlich die Konstanten 
a=10, b=5. 
Damit ist im wesentlichen die Gestalt des Dammes festgelegt. In der folgenden Tabelle 


i—1 
S. 260 gibt Spalte (VI) die Werte A-5)(0 NA an. Die Spalte (VII) liefert 


1 
ra 


3 
(A — 5) (10—7)4 


Sie setzt sich zusammen nach der Gleichung 


(VID = (IN: (IV)yıVD. 


In Abb. 5 ist die Abhängigkeit der Funktion e’-cos 5 von y veranschaulicht. 


e’. sin ist durch Spalte VIII bekannt, die sich zusammensetzt nach der Gleichung 


(II) - (V)-y(Vl). 


Abb. 6 zeigt den Integranden e*’- sin ; in Abhängigkeit von y. 


92 
SEN 
8 | 
/ 


AN 
0,7 01 
\ \ 
U 


7 3 Y 0 1 2 
Abh. 5. Abb. 6 


In den beiden Zeichnungen sind die durch graphische Integration gefundenen Werte des 
inneren Integrals für = 0,5; 0,8; 0,9 wegen ihrer Genauigkeit unentbehrlich. 


Die graphische Methode liefert nämlich numerische Werte des Integrals für nicht zu 
singularitätennahe #4 müheloser und genauer als die Berechnung des inneren Integrals mit 
Hilfe der Sehnenapproximation. Für diese A sind aber die Anforderungen an die Genauigkeit 


besonders groß, weil an diesen Stellen für et. cos I und e’sin 2 ein Maximum vorliegt. 


(serade hier verursacht ein geringes Schwanken für \e’ cos, dy und \e’sin{, dy einen 


beträchtlichen Fehler. 


Für singularitätennahe % dagegen verliert die graphische Methode an Einfachheit. Wie 
aus den Bildern 5 und 6 und aus der Tabelle S. 260 hervorgeht, werden die Kurven aber 


| 

_ 


- f.angew. 
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F + Y98z EOECO eel9 000 - SIT 06 
F | 2100| | 22000000°0 SFIEO 0916“ W600 — | 
0 2600 22000000 E296°0 SI 800 — | 200 198 66660 
rIEO SOTO VEHED 0086 ro I 900 IO8°0 660 
2910 1% OSTO 11910°0 LELIO 00 — | 900 cc] 090 
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0,2632 + 0.0 


0.1932 + 0, 


‚000 


0.1795 +0 


004 


0,2866 + 0.00 


0.00647 
0.002098 


0,004 


0.013 


4.059 


6.123 


0.001 


- 1.401 


0.003 
0,001 


— 2,936 


1,268 
3,820 


1.118 


0,80 
0,00 


0.8057 0,5924 


0.8478 0,5303 


> 


06 


— 1.825 0,002 
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sehr wesentlich durch Werte von 4 bestimmt, die äußerst dicht bei 4-1 oder 4=0 liegen. 
Wegen der Steilheit der Kurven kann hier der mögliche Fehler die Integrale nur wenig 


heeinflussen. 


Durch Auszählung der Fläche werden die Integrale \e cos /, dy (Spalte (IX)) und 
\e’sin 5 dy (Spalte (X)) bestimmt. Die den Zahlen angefügte Korrektion läßt die mögliche 
Schwankung der Kurve zu. Die Kurve, die die freie Oberfläche darstellt, bestimmt sieh aus 
diesen beiden Integralen, die für ein bestimmtes y die Koordinaten eines Punktes der freien 
Oberfläche sind. 

Die ausgeführten Rechnungen sind zum großen Teil auch für ein anderes Konstanten- 
paar a und b zu verwenden. Erst Spalte (VD), die wesentlich dureh die Konstanten a und b 
bestimmt wird, und die folgenden zwei Spalten (VII) und (VIII), zu deren Berechnung Spalte (VI) 
gebraucht wird, müssen für jedes Konstantenpaar a und b neu berechnet werden. Anschließend 
erfolgt dann die graphische Integration. 


IV. Das innere Integral für 
Für die beiden angegebenen Intervalle wurde das innere Integral auf graphischem Wege 


bestimmt. Der Integrand n . wurde als Funktion von # für die verschiedenen 4 gezeichnet 


und das Integral durch Auszählung der von der Kurve und den Achsen umgrenzten Fläche 
erhalten. Das Ergebnis der Integration zeigt die folgende Tabelle. Die angegebenen Fehler 
berücksichtigen die mögliche Schwankung der gezeichneten Kurven. 


l 
»(t) 3 v(t 
0 
1.05 — 391 +0,01 + 1.870 + 0,005 6,488 
1.10 — 3,20 0.01 1.529 0,005 4.616 
2.00 -1.04 0,01 0,496 0,005 1.643 
3.00 — 0,620 0.008 0.296 0.004 1.344 
4.00 0.444 0,006 0.212 0,008 
5.00 — 0,343 0,005 0.164 0.0083 1.178 
6.00 0,280 0,005 0,134 0.003 1.143 
7.00 0.239 0.005 0,114 0.003 1.121 
8.00 -0,208 0,005 0,099 0,008 1.104 
12,00 0.136 0,0008 0.065 0,0004 1.067 
15.00 0,108 0,002 0,052 0,001 1.05: 
20,00 0,080 0,0002 0.038 0,0001 1.039 
— 0,01 902 0,06 -4305 0,08 0.0135 
0,10 4,30 0,02 — 2.053 0,01 0,1283 
— 0,20 3.12 0,01 — 1,490 0,005 0,225 
— 0,40 2.09 0,01 0,999 0,005 0,368 
0,50 1,84 0,01 — 0,879 0,005 0.415 
0,70 147 0,01 0,702 0.005 0,495 
1.00 1.14 0,01 0,545 0,005 0,582 
2,00 0,65 0,01 — 0,310 0,005 0,733 
— 3,00 048 0,01 0.229 0,005 0,796 
5.00 0.29 0,01 0,138 0.005 0,870 
— 8,00 0,187 0,005 0,089 0,003 0.914 
- 15.00 0,101 0,002 — 0,048 0,001 0,953 


20,00 0,077 0,002 0,037 0,001 0,964 


| 
| | 
| 
| 
% 
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@—b)(a—/)) 4 
reduziert sich bis auf einige Konstanten der Ausdruck (A) auf Seite 355 
Y(t) 
auf das oben angegebene Integral, denn für große A strebt e 27,14 


Für große |/ 

gegen 1 (Abhandlung 

S. 142) und d(w -iv) wird nach Gl. (B) (S. 255) zu ag. Wegen Gl. (4) aber ist 
| 

E (Abhandlung S. 113). Das Integral läßt sich durch Umformung auf eine gebräuch- 


liehere Form bringen 


aresin | a(dı) 
/ di _ „ı/1-—a do 


VI. Berechnung der Höhe AB des unteren Wasserstandes. 


Ks ist zu berechnen („Abhandlung* S. 143) 


/ 


Das Intervall 10 °.4°. wird zerlegt. Je nach der Größe von 4 werden verschiedene Be- 


reehnungsmethoden herangezogen. 


1. Das Intervall 10.4. 11. Für große |, gilt nach „Abhandlung“ (S. 143) 


1 


di 


Unter Benutzung der Tafel von Hayashi findet man für 71 in der Umgebung von 4 = 10 

den Wert 0,0335. Für = 10 ıst 7+ 0,105. Daraus folgt = R = 0,0002 


Das obige Integral läßt sich also im angegebenen Intervall zerlegen in drei Integrale 


(unter Weglassung der Konstanten): 


11 

10 10 10 

1 


a) e' 


10 


e' wird in zwei Faktoren zerlegt 


1 
| 1 | | 


Der zweite Faktor, noch durch 4 dividiert, erweist sich in der Umgebung von A = 10 nahezu 
als Konstante im Vergleich zum ersten Faktor. Für 18 =4=- 11 kann also 


1 
/ 


dureh 0,0972 + 0,0055 ersetzt werden. 
Der erste Faktor in (6) liefert das Integral 
1 


10 


— 


| 
262 | 
| 
| 
| E 


W. 
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Within ist 


11 
— = 0,5674 (0,0972 + 0,005) — + 
10 
11 1; 
b) 
10 


Wir verfahren hier genau wie beim letzten Abschnitt und setzen 


1 17.2718 4 


Mithin ist 


11 
(0,006 0.001285 + 0.000155. 
11 
e) R\eran. 
1 


Wir setzen 


0 


1 
3 
37,14 1,019 + 0,001 


| 
für und erhalten, da 0,0002 | 


11 
dä 0,867 (1,019 + 0,001) 0,0002 5 — O,OOOLTT 


10 


Die Zusammenfassung der drei Integrale führt zu 


11 


2. Das Intervall 11 4=-20. Die Integration wird für dieses Integral graphisch aus- 
geführt. In der Tafel von Hayashi ist 5 in Abhängigkeit von %*” tabuliert, wobei nach 
Seite 256 


Damit ist die Beziehung zwischen 5 und 4 bekannt und in der folgenden Tabelle angegeben. 


0,1 10 1.5088 14.286 1.7175 
0.05 20 1.8294 0,065 15.384 1.7420 
0,095 10.526 1.6161 0,060 16.667 1.7683 
0,090 11.111 1.6341 0.055 15.015 1.7968 
0.080 12,500 1,6733 0,050 2) 1.8204 
0,075 13.333 1.6947 
(t 


Die folgende Tabelle zeigt die Art der Berechnung des Integranden % 7 


t | 
| 
| 
_ 
l | 
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(4 (A —- 10) 4 | (4-5) A-1W)%A 
| 
11 0,38925 1.075 0,4180 031 
12 1.067 0,2727 T 
13 0,19612 1.062 0,2085 e DET N 
14 0,16060 1.0750 0.1698 
15 0.13662 1.053 0,1439 
0.11920 1.050 0.1251 
17 0.105855 1.047 0,1108 
IS 0.095290 1.044 0.0995 = 
10 0,08671 1.041 0,0903 Abb. 7. 
1.039 0.0827 


In Abb. % ıst der Integrand als Funktion von gezeichnet. Die Auszählung der Fläche ergibt 


wenn eine mögliche Schwankung der Kurve berücksichtigt wird. 
1 
3. Das Intervall 20°. Wir setzen e =14008 für denn 
die Potenz strebt nach I für > vw und hat für 420 den Wert 1,04. Demnach reduziert 
sich \e’dg auf 


i= 


(140,03) 1 4009 


20 


23 
(A—b)(k 


wenn dgq ij (140,03) angenommen wird. Man muß bei dieser Annahme bedenken, daß 


arte 


und daß für 20°-2°- x die Klammer höchstens den Wert 1-+0.03 annimmt. Das Integral 


| 


läßt sich nach einer auf S. 262 angegebenen Umformung berechnen. Es ergibt sich 


\edg (140,08) 14000 
Durch Addition der Teilintegrale (7), (8), (9 folgt 
+ 0,0055 AB. 
W 
VII. Berechnung der Sickerstrecke BC. 
Zu berechnen ist \ e' dp. 


Wir zerlegen das Integral in zwei Teile. 


1. Das Intervall 20°-.7°-0. Wir verfahren hier wieder 
wie beim Intervall I1°-2°-20, Seite 263. In Abb. S ist der 
Integrand 


ar 
—10)4 
als Funktion von 9 gezeichnet. Die Art der Berechnung des Pe 
Integranden ergibt sich aus der folgenden Tabelle. 


1 | 
et | 
1 
# 
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; ) 
0.01 0,426 2.011 0,0135 0,0272 
0.1 0,611 0.214 0,128 0.0547 
0.2 0,703 0.1133 0,225 0,0749 
0.4 0,813 0.0624 0.368 0,0921 
0.5 0,857 0,0519 0,415 0,0045 
0,7 0.920 0.0308 0.496 0.0990 
1.0 1.000 0.0303 0,582 0,1014 
2.0 1.165 0.0176 0,733 
3.0 1,279 0,01282 0,796 0,0002 
— 50 1.420 0,871 0.0779 
1.566 0.004850 0,914 0.0033 
15.0 1.760 0.002138 0,953 0.0440 
20,0 1.45 000140 0.064 0.0361 


Die Beziehung zwischen 4 und 5 in den beiden ersten Spalten liefert die folgende Tabelle: 


k? . k? B 

0.01 0,010 0,426 0,65 1.860 1.151 

0,04 0,0416 0.526 0,70 — 2,335 1.210 

0.039 0.0406 0.524 0,75 3.000 1.279 

0.038 — 0,0395 0,521 0,80 4.000 1.360 

0.1 0.111 0.626 0,832 5.01 1.120 | 
0,831 — 4,98 1.119 

0,15 0.1765 0.684 

0,20 0,2500 0,735 5.669 1.161 

0.25 - 0,3332 0,781 0.90 3,000 1.597 

0.30 0.428 0.825 0.05 19.000 1.821 

0.35 0.538 0,872 

0,40 0,666 0.011 

0.45 - 0,818 0.0955 

0.50 — 1.000 

0.55 1,221 1,045 

0,60 - 1,500 1.097 


Aus Bild 8 folgt 0,109 + 0,002. 


2. Das Intervall © 4.20. Hier gelten dieselben Betrachtungen wie bei V1, 3:5. 264 


1 
e wird jetzt gleich 1 gesetzt. Dann ist 
— 10 — 
[AB: 
- .03 — 0.0445 0.014875 0.0002 BC: Sicherstreche 
(1 +0,03) (1 — 0,04) 145 01484 — CD: freie Oberflache 
B 
Durch Addition der beiden Teilintegrale folgt 
Schicht 
+ 0,W2i- BC. — 
A 
VII. Schlußergebnis. Abb. 9. 


Aus BC, AB, DC ergibt sich Abb. 9. Für einen Spezialfall, bei dem zwei Konstanten 
(@a=10; b-5) willkürlich gewählt wurden, ist nun die wirkliche Gestalt bestimmt worden. Die 


drei Größen BC, AB und DC sind dabei bis auf den Faktor € 3 bekannt. Durch Wahl der 


Konstanten € erhält man also die obenstehend angegebene Abbildung in der gewünschten Größe. 
Die Längen von BC und AB sind auf S. 265 und 264 angegeben. Die freie Oberfläche DE 
bilden bekannte Punkte, deren Koordinaten dureh die Spalten (IX) (Abszissen) und (X) (Ordinaten) 


- 
bis auf den Faktor -, bestimmt sind. Sal 
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On the Transition Curve Considering the Fffect of Variation 
of the Train Speed. 


By Y. Tanaka in Tokio. 

Zusammenfassung. Die vorliegende Arbeit beschäftigt sich mit der Aufgabe, die 
Überzangskurve aus der Geraden in eine Krümmung für ein Eisenbahngleis unter den 
folgenden Voraussetzungen zu bestimmen. Es wird angenommen, daß die Zuggeschwindigkeit 
von einem höheren Werte (in der Geraden) zu einem niedrigeren Werte (in der Krümmung) 
linear mit der Bogenlänge abnimmt, und daß gleichzeitig die Überhöhung der äußeren Schiene 
linear mit der Bogenlänge zunimmt. Dann liefert die Forderung, daß die Gleise keinen seit- 
liehen Schub erfahren sollen, eine Differentialgleichung für die Übergangskurve, welche sich 
integrieren läßt, wenn man annimmt, daß (bei Einführung eines geeigneten Koordinatensystems) 
die Krümmung der Übergangskurve durch den zweiten Differentialquotienten ersetzt werden kann. 

I. Introduetory. The objeet of the present artiele is to submit a mathematical solution 
of a new transition eurve, eonsidering the effeet of the linear variation of the train speed, 

We know indeed that the eubie parabola is the theoretieal transition eurve so far as the 
train speed is uniform and the outer rail has a linear gradient over the same curve, while 
after my opinion, it might not be of no use to investigate the effect of the speed variation 
to the form of the transition curve, and we should like point out that by the present 
eonsideration we can get one of the more general solutions of the transition curve, which 
mieht bring something new to the practical applications. 

Il. Fundamental Equations. In the first place, let us assume that the speed variation of 
an axle with a pair of wheels as shown in Fig. 1. "Then, we may write the speed of the axle 
runing over the eurve as below. 


where Mm 


In this expression, we may note that *a* is positive in case 1, while it is negative in 
ease 2, but it makes no difference in the following mathematical operations. 

In the next, we shall make two usual assumptions in regard to the gradient of the 
outer rail and the amount of the superelevation as shown in Fig. 2 and Fig. 3, whence we 
may write the following equations. 


where  superelevation of the outer rail, length of the transition curve, k, super- 
elevation on the eireular eurve. 


L 
Cırcular curve /ransıflor curve (IrCular curve 


Fig. 1. Pie. 2, Fig. 3. 
And as the superelevation of the outer rail at any point is to be assumed as to make 
the direetion of the resultant force *R* normal to tlıe track surface as shown in Fig. 3, we 
may write without any praetical error as 


| 

where  gauge of track, IN total weight under consideration, eentrifugal force. And 

as the centrifugal force Ü = —, in which m = mass and o -radius of eurvature; we may 
rewrite (9) as G 

g9=gravity acceleration . . 2.2.2.2... 


Therefore, putting (1) into (4), we get 


(", =. Ad 
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I dy _ 
In the present case, as we may generally put 7 = $ from (2) and (5) we shall have 
dır dp? x 
de Gi (v,-ax) d.e? (1—-ax) 
a 
where = | - and 
v, Giv,? 


This is the fundamental equation for the present problem, and the ordinates of the 
new eurve will be given by the solution of this equation. 


III. Solutions of the Fundamental Equation. By integrating the equation (6), we get 


dy Kf \ 
Ä 
in which the integration constant e, will be given as ce, = — I by the eondition, Ze Vater 0. 


And by integrating once more, we have 


5 


where as y- Dat 2=0, we get c,—=2 or 
K 


— (2 ar) —2ax} (S), 


r 4 1 
where we know that K= and a: 
Giv,? 

This is the equation of the new transition eurve and its tangent will be given by 
equation 

In the foregoing solution we have made eonsiderations only for a single axle for the 
sake of simplieity, while some practical cases as for single car operation ete. will be also 
covered by the same assumption. 

In the present, it will not be needless to note that, if we put «= 0 after the expansion 
of the equation (S), we get the following expression of a cubie. parabola, 

y= 4 
And in this eubie parabola, K is a function of v, and not of v,, therefore, if we should like 
caleulate the eurve for v,, we have to rewrite this as 


IV. Numerical Example. To find the numerical ordinates of the new transition eurve 
of equation (8) by the following data: 
?, 40 km/hr, #,=50 km/hr = 50000/3600 m/see, 
m, G=1067 m and k,= mm 0,0HS m. 


In the present case, we have 


and we shall have the following results. 
(m) — (2 Dax y(m) 
0.500027 0.05 0.001 
10 0.1000219 0.10 0.009 
15 0.1500761 0.15 0.052 
20 0.2001550 — 0.20 0.079 
25 0.2503712 — 0.25 0.158 
30 0.3006602 0.30 0.251 
30 0.351075 0.39 0.459 
40 0.4016583 -0.40 0.706. 
The max. ordinate of the new eurve given by tlıe foregoing example are about 79.7"), of 
that given by the usual curve =5;K „a. Ihe tangent angle at the end of the new 


0 
curve for this case is 3° 17’ and that of the usual ceurve is 3° 47’, so that the former ıs about 


6°, of the latter in this case. 33) 


| 
8, 
(8), 
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may 
(4). 
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Die Knickfestigkeit rechteckiger Platten. 


Von ©. H. Faxen in Göteborg (Schweden). 


1. Einleitung. Für die Berechnung der größtmöglichen Belastung, welche rechteckige 
Platten ohne Gefahr für Auskniekung tragen können, gibt es verschiedene neue Methoden, 
die aber alle zu nieht miteinander übereinstimmenden Ergebnissen führen. Da keine der 
numerischen Methoden streng bewiesen wurde, habe ich es unternommen, nachzuprüfen, welche 
Methoden anscheinend riehtig und welche unrichtig sind. Die Ergebnisse sind in Kapitel II 
zusammengefaßt. 

Vom rein mathematischen Standpunkte aus genauere Lösungen des Problems kann man 
allerdings leieht finden. Da sieh diese aber, wie ich anfangs glaubte, für numerische Rech- 
nungen nieht eigenen, so sollen sie diesmal des mangelnden Platzes wegen unerwähnt bleiben. 

Die Schwierigkeit des Problemes hängt davon ab, wie die Ränder der Platte befestigt 
sind. Wenn zwei parallele Ränder gestützt sind, ist die genaue Lösung einfach, sei es, daß 
die übrigen Ränder frei, gestützt, eingeklemmt oder elastisch gelagert sind (Bryan, Reißner, 
Tımoshenko usw.) 

Um die Sachlage in einem Falle zu studieren, wenn die genaue Lösung kompliziert ist, 
wollen wir im folgenden das Beispiel behandeln, wo alle Ränder eingeklemmt sind. 

Für ein solehes Problem wurde meistens das Ritzsche Verfahren angewandt. Dieses 
Verfahren besteht bekanntlich darin, daß man für die Form der beginnenden Ausbeulung der 
Platte einen mathematischen Ausdruck wählt, der die Randbedingungen vollkommen und die 
Ditferentialeleichung der Platte so genau erfüllt, wie dies mit einem leichthantierlichen Aus- 
druck möglieh ist. Man arbeitet dann allerdings nieht mit der Differentialgleichung selbst, 
sondern mit einem entspreehenden Minimumprinzip. Der Erfolg hängt davon ab, ob der er- 
wähnte Ausdruck gut gewählt ist, ob man also betreffs der Form des Ausdruckes eine glück- 
liche Vermutung gemacht hat. 

Einen anderen Weg hat K. Sezawa') betreten, indem er den Ausdruck so wählte, daß 
die Differentialgleichung und einige Randbedingungen genau erfüllt wurden. Der Ausdruck 
enthält einige Konstanten, die so bestimmt werden können, daß die noch übrigen Rand- 
bedingungen in zweckmäßig gewählten Punkten der Ränder genau erfüllt werden. Es zeigte 
sich dann, daß diese Randbedingungen überall ziemlich genau erfüllt wurden, und daß man 
so eine für technische Zwecke recht brauchbare Methode für numerische Rechnungen erhält. 
In Vergleich mit der Ritzschen hat diese Methode den Vorteil, daß man beim Auffinden 
der zweekmäßigsten Ausdrücke nieht auf eine glückliche Vermutung angewiesen ist, sondern 
daß diese sieh bei der Auflösung der Differentialgleichung darbieten. 

(. 1. Tavlor?) ergänzte den Ausdruck von Sezawa zu unendlichen Reihen mit un- 
endlieh vielen Konstanten und zeigte an einem Beispiel durch numerische Rechnungen, wie 
man einen guten und anscheinend riehtigen Wert für die höchstmögliche Belastung erhalten 
kann, welche quadratische Platten ohne Gefahr für Auskniekung tragen können. 

Wenn dieser riehtig ist, sollte also die Taylorsche Methode die Lösung eines anerkannt 
schwierigen mathematischen Problems ermöglichen, eine Lösung, die allerdings nur vom tech- 
nischen Standpunkte befriedigend ist. Zur Nachprüfung seiner Ergebnisse soll seine Methode 
im folgenden abgeändert und auf verschiedene Probleme angewandt werden. Es wird sich 
zeieen, daß die auf verschiedenen Wegen erhaltenen numerischen Ergebnisse miteinander 
übereinstimmen. Der Verfasser beabsichtigt. die Untersuchungen fortzusetzen. 

Die Rechenmethoden werden in Kapitel I und II angegeben. Der größte Teil der nume- 
rischen Ergebnisse ist in Kapitel II zusammengestellt. 

Durch die Unterstützung des CO halmerschen Forschungsfonds ist die vorliegende Unter- 
suchung ermöglicht. 


Kapitel I. 


$ 1. Darstellung der Funktion »" durch trigonometrische Reihen. | bedeutet im folgenden 
den Differentialoperator 
0x 


1), ZAMM, S. 227 bis 220, 19932. Vgl. auch Nadais Dissertation, die Verf. leider nicht gesehen hat. Ferner 
Inglis: Inst. Naval Architeets, 1925, wo weitere Literaturangaben zu finden sind. 


*), ZAMM, S. 147 bis 152, 1955. 
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ie Differentialgleichung der Überschrift lautet also in ausführlicherer Schreibweise 


; ist eine Konstante. 

Jede Lösung dieser Differentialgleichung läßt sich in eine Summe von zwei Funktionen 


aufspalten, so daß und » Lösungen der Differentialgleichungen Iu 0, Ir sind, 
und zwar sieht man leicht, daß 


+72 


ist. 

Damit die Formeln nicht unnötig lang ausfallen, wird vorausgesetzt, daß ww (we, ) eine 
verade Funktion sowohl in bezug auf = wie in bezug auf y sei. Für jeden Wert von = und y 
oilt also 


Ähnliche Formeln gelten dann auch für und (w, 

Da u die Differentialgleichung Ju 0 erfüllt, d. h. eine harmonische Funktion ist, liegt 
der Gedanke nahe, diese Funktion als eine trigonometrische Reihe zu schreiben. 

Wir machen den folgenden Ansatz, durch den jede im Bereich 


reguläre Lösung der Gleichung Au =-0 dargestellt werden kann: 
\ "a ,„eoshn.r cosn y—+ y 


n=] 


(>). 
Ebenso gilt: Jede im Bereiche (4) reguläre Lösung der Differential- 


sleichung Ar +/’v=0 kann so 


A,cosi# B,c0siy 


\ 42) cosn y-+B,cosn.wcosh (yyn?-- (6) 


seschrieben werden. 

Da die Funktion »w= nach (2) eine Summe einer Funktion « und einer Funktion » mit 
der eben angegebenen Eigenschaft ist, haben wir die allgemeine Form der Lösung der ın der 
Kapitelübersehrift stehenden Differentialgleiehung (1) gefunden. 

Das Ergebnis ist eine alternative Form zu der von Taylor gegebenen Lösung derselben 


Differentialgleichung. 

$ 2. Das Eigenwertproblem der Differentialgleichung (1). Ks wird eine im Inneren des 
(rebietes (4) reguläre Lösung 9° (x, y) der Differentialgleichung (1) gesucht. 8 Gr, y) soll ferner 
die Randbedingungen 


w(tr,y)=d, 
8) 


erfüllen. 

Dieses Problem besitzt nur für einzelne Werte von /, welche man die Eigenwerte nennt, 
eine oder mehrere Lösungen, die Eigenfunktionen. Wir wollen diese zu finden versuchen. 

Nach der Gl. (2) kann man ww in eine Summe «+» aufspalten, wobei « und » die in 
dem vorangehenden Paragraphen angegebenen Eigenschaften haben. Wie die Funktionen « 
und »= sich am Rande verhalten, wissen wir aber nieht. Es bleibt also eine durch die folgen- 
den Rechnungen zu beweisende Annahme, daß « und » von der Form (5) und (6) sind, und 
zwar wird dieser Beweis diesmal nicht vollständig durchgeführt werden. 


| 
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Die Koeffizienten der Reihen (5) und (6) und die Konstante 4 sollen so bestimmt werden, 
daß die Randbedingungen (7), (5) und (9) erfüllt sind. Die Randbedingungen (8) und (9) er- 


fi rdern, daß Ir um A, sin / B, sin in, A, 


ist. Die Randbedingung (7) gibt zwei Gleichungen, welche zur Bestimmung der Koeffizienten 
dienen. Um möglichst kurze Formeln zu erhalten, beschränken wir uns ferner darauf, nur 
solehe Funktionen « und » zu suchen, welche unverändert bleiben, wenn x und y vertauscht 
werden, d.h. wir suchen nur solche Lösungen, in denen u (@,y) =u(y,#), v =ev(y,x) ist. 
Dann ist 


coshnxcosny-+ cosnzcoshny 


(9), 
nsınh n 
cosh(e yn? -4)cosn y+ cosn.wcosh (yyn? 
== sinh (a — 4°) 


Setzt man diese Reihen in die Randbedingung (7) ein, so erhält man zwei Gleichungen, welche 
mit der folgenden gleichbedeutend sind 


\ ( coshn.x + cosn.r cosh nz 
n sinhn 


n—1 
( eosh Yn? — 7°) + cos n cosh (a Yn? — 7°) 

Diese Gleichung muß für jeden Wert von x erfüllt sein. Jedes Glied kann nieht gleich 
Null sein, denn dann wäre =0, Da die Funktion y) nieht überall gleich 
Null sein soll, muß also 4 von Null verschieden sein. Man soll also diejenigen Werte von 4 
aufsuchen, für welche die Gl. (10) erfüllt sein kann, ohne daß die Zahlen €, a,, «a,, 
Ay, ... alle gleich Null sind. Diese Werte von / sind die gesuchten Eigenwerte des am An- 
fang dieses Paragraphen gestellten Problems. 

Um die Eigenwerte numerisch zu bestimmen, kann man drei naheliegende Wege ein- 
schlagen. 

Der eine ist eine Verallgemeinerung von Sezawas Methode. Man berücksichtigt nur 
eine endliche Anzahl Glieder der Reihe in (10) und bestimmt / so, daß die Gl. (10) für 
möglichst viele zweckmäßig gewählte Werte von = genau erfüllt wird. 

Die zweite, der Taylorschen nachgebildete, Methode ist, die rechte Seite der Gl. (10) 
mit cos m. zu multiplizieren und von — 7 zu + zu integrieren. Indem man m der Reihe 
nach gleich O0, 1. 2,3, ... setzt, erhält man eine unendliche Anzahl von Gleichungen für die 


Bestimmung von a,, d,, ... nämlich: 
(11) 
i 
worin (- yn + 
mn (mn?) (m?+ — 4°) (m (13), 
m 2 m? (2 m?’ — 4°) 


ist. Gelingt es, aus den Gl. (12) die Koeffizienten a,, d,, @,, ... zu bestimmen, so erhält 
man (C aus der Gl. (11). Der nächste Paragraph enthält unsere Methode, eine für praktische 
Zwecke hinreichend genaue Lösung der Gl. (12) zu finden. 


$ 3. Die numerische Berechnung der Eigenwerte. Um die Eigenwerte, d. h. diejenigen 
Werte von 4. für welehe die Gl. (12) nieht verschwindende Lösungen «a,, a,, A,, ... haben, 
zu suchen, wird man zunächst diejenigen 4#Werte aufsuchen, für welche das entsprechende 
endliche Gleichungssystem 
N 
n—| 
nicht verschwindende Lösungen besitzt. Falls nun diese 4-Werte Grenzwerte zustreben, 
wenn N ins unendliche wächst, wird man erwarten, daß diese Grenzwerte die Eigenwerte sind. 
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Aus numerischen Rechnungen erhielt G. I. Tavlor*) den Eindruck, daß es sich wirk- 
lich so verhielte. Er hat die Reehnungen anders gemacht und hat anstatt (12) ein anderes 
Gleichungssystem erhalten. Er suchte den kleinsten 4-Wert auf, für welchen das entsprechende 
System von N Gleichungen mit N unbekannten Lösungen hat. Er fand, daß sich bei wachsen- 
dem N ) rasch einen Grenzwert anzunähern schien. 

Unsere Rechnungen sind, wie gesagt, anders als bei Taylor geführt. Wir können 
laher sein Ergebnis prüfen. Zu dem Zwecke untersuchen wir erstens, ob der kleinste 4-Wert 
auch hier bei wachsendem N einem Grenzwert zustrebt. Wenn es sich so verhält, wollen 
wir zweitens prüfen, ob dieser Grenzwert mit Taylors Ergebnis übereinstimmt. 

Ingenieur A. Anderfeldt hat numerische Rechnungen ausgeführt, deren Ergebnisse in 
der folgenden Tafel zusammengestellt sind: 


2? 1000 D, 1000 D, 10000 D, 10000 D, 100.000 D, 
1.3 167.28 | — 28,00 973 14 

1.323 1.857 | — 2,678 | 15 — 
1.325 15.588 | — 0,609 | — 8,447 0,5128 - 0,3462 
1328 3697 2.61 81 25 7.61 
1,3228 | 13255 1,3259 1,3252 1.3252 


Die Bedeutung von Dy ist die Determinante der Gl. (15), d. i. 
DIN 
&N 
Dx= 


CN1, CNN 

Der kleinste Wert von 4° wird ausgewählt, für welchen diese Determinante gleich Null wird, 
denn für diesen 4°-Wert besitzen die Gl. (15) die gewünschte Lösung. Aus jeder Vertikalreihe 
der Tafel kann man durch graphische Interpolation den Wert von 7° für N=1,2,3,4 und 5 
berechnen. Dieser ist in der letzten Horizontalreihe der Tafel angegeben. Man wird zunächst 
geneigt sein zu erwarten, daß eine nachfolgende Vertikalreihe immer genauere Werte als eine 
vorangehende ergeben soll. Der Verfasser fürchtet aber, daß die vierte und fünfte Vertikal- 
reihe nieht bessere Werte als die dritte ergeben haben, weil nur fünfstellige Tafeln der 
der trigonometrischen und hyperbolischen Funktionen benutzt wurden. Die unvermeidlichen 
Fehler in der letzten Stelle haben einen vielleicht größeren Einfluß in den viel längeren 
Berechnungen der Determinanten D, und D,. Der Verfasser möchte den Wert 47° 1,326 als 
Eirgebnis der Berechnungen betrachten. 

Das Ergebnis stimmt befriedigend mit Taylors 442°=-5,30. Verschiedene Methoden 
führen also zu übereinstimmenden Ergebnissen. Es gibt keinen Grund, die Richtigkeit dieser 
Werte zu bezweifeln. | 

Sezawa rechnete nur in der ersten Annäherung und fand den Wert 4#°=5,61. 

Für eine weitere numerische Prüfung wurde 4? auch mittels einer vierten Methode aus: 
gerechnet. In Vergleich mit der voranstehenden schließt sich diese Methode den Taylor- 
schen Entwieklungen getreuer an. Sie weicht jedoch darin ab, daß die trigonometrischen 
Reihen durchgehend nach ungeraden Vielfachen von « und 9 fortschreiten, währenddem 
Taylor in seinen Gl. (7) und (8) zu geraden Vielfachen übergeht. Unsere Methode wird im 
nächsten Kapitel für ein anderes Problem näher entwickelt werden (Kap. II, $ 2). 

Wenn nur die erste Gleichung und die erste Unbekannte berücksichtigt wurde, fand ich 
für 2° den Wert 1,29. 

Kapitel Il. 
Die Differentialgleichung 4.11. +4 


$ 1. Die Randbedingungen. Die Differentialgleichung soll für einen rechteckigen Bereich 
mit dem Seitenverhältnis k gelöst werden. Die Größe des Rechteckes kann man beliebig 
wählen, weil die Lösungen des zu behandelnden Problemes für ähnliche Rechtecke sich nur 
dureh einen Proportionalitätsfaktor unterscheiden. Um möglichst kurze Ausdrücke zu erhalten, 
wählen wir den Bereich 


_ 


3), A.a.0. 


) 
| 
2 
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Die Lösung mw der Differentialgleichung soll an den Rändern des Bereiches folgende Bedingungen 
erfüllen: 


w (a, für + , 9 beliebig und für y= + 5, beliebig. . . (17), 
al 


$2. Die in bezug auf . geraden Lösungen. Unabhängig von den Werten der Koeffizienten 
4, und db, erfüllt jedes Glied der folgenden Reihen die Differentialgleiehung und die Rand- 
bedinzuneen (17) 


cosh p, .r eosh p, .r - feoshq,y  coshq,+ 
m \ \+ \ h,eos Kır (20). 


_ 


In diesen Reihen nehmen m und n» nur die Werte der ungeraden ganzen Zahlen an. Die 
Bezeichnungen haben die Bedeutung 


Die Randbedingungen (15) ergeben die Gleichung 


Y 
ak 
\ 
. na fcoshy,y  cosh 
2 
Die Randbedingungen (19) ergeben 
. ma fecosh p, cosh p, 
0 
cosh ecosh 
v 
4,7 
\ h,eos (RK th | 


Die Gl. (23) und (24) dienen zur Bestimmung der noch unbekannten Koeffizienten 
BB rer Die Gleichungen müssen für jeden Wert von = und y im Bereiche 
(16) erfüllt sein. Sie entspreehen der Gl. (10) des ersten Kapitels. Man wird erwarten, 
daß (25) und (24) auf dieselben Weisen gelöst werden können. Diesmal haben wir uns aber 
auf die Berechnung numerischer Werte beschränkt und haben nur die folgenden Methoden 
näher verfolet. 

Für einen beliebigen Wert von 4 gibt es keine andere Lösung der Gl. (23) und (24) als 
a, Der kleinste Wert von wird gesucht, für den diese Koeffizienten 
nicht alle gleich Null sein müssen, damit die erwähnten Gleichungen erfüllt werden. 

Betrachtet man immer noch in den Gl. (23) und (24) a,„ und 5b, als die Unbekannten, so 
sind diese mit veränderlichen, von = und y abhängigen, Zahlen multipliziert. Um konstante 
Vorzahlen zu erhalten, kann man die Gl. (23) mit cos m,y und die Gl. (24) mit cos K, x 
multiplizieren und dann jede Gleichung von 4° 0 bzw. &=-0 bis /2 integrieren. Dabei ist 
N, n,/k. m, und», sind ungerade ganze positive Zahlen. Läßt man nachher die Beizahl 1 
aus und führt die Bezeichnungen’) 


. m? — 2° 


7 — 4? (25) 
Sm? cos (ki) + ecosh m’ 2°) 
10, 
) teh 4,teh 26) 


') Griechisches @ und lateinisches a haben verschiedene Bedeutung! 
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m-+n 

(27) 

A m k? m ] m’ B, h,, . . . . . . . . . . . (25) 

ein. so erhält man nach elementaren Rechnungen das Gleichungssystem 

dm +2 Cm 7 (m +), Pn +, ( mn A (n 3, .. .) (29). 
a=1,3,5,.:° 


Die Aufgabe ist nun, diejenigen Werte von 4 zu bereehnen, für welehe dieses System 
von Null verschiedene Lösungen A,.4A;...., B,PB;,... besitzt. Um numerische Werte zu er- 
halten, haben wir als eine erste Annäherung die kleinsten Werte von 4 ausgerechnet, 
für welche das abgekürzte Gleichungssystem 

a, A, +c,B,=0, ER. 
eine Lösung A, B, besitzt. Als zweite Annäherung wurden diejenigen kleinsten Werte von 4 
bezeichnet, für welche das Gleiehungssystem 

aA, +c,B tes B+%,A, +0, B=0 


eine Lösung besitzt. Dieses Gleichungssystem erhält man aus (29), wenn man alle Zahlen 
mit höheren Beizahlen als 3 gleich Null setzt oder diese Zahlen nieht berücksichtigt. 
Bezeichnet man die Determinante des Systemes (30) mit D, und des Systemes (31) mit 
D,. werden also diejenigen Werte von 4 gesucht, für welche folgende Ausdrücke gleich Null 
werden: 
in erster Annäherung: D, 
in zweiter Annäherung: 
D, (a, Pi (a; Ps es) (a, Ps Cz + Az 6, Cz 2 Cı3 Ca 
D, und D, wurden für verschiedene Werte von 4 berechnet, und dann wurde der kleinste 
Wert von 4 ausprobiert, für den D, und D, gleich Null werden. Die numerischen Ergebnisse 
sind in Kapitel III zusammengestellt. 
$ 3. Die in bezug auf x ungeraden Lösungen. ‚Jedes Glied der Reihen 


sinh p, sinh 
m=]1,3,5, sınh sinh | 
\ b„sın Kır 


erfüllt die Differentialgleichung und die Randbedingungen (17). Die Bezeichnungen sind die- 
selben wie in $ 2. Ein Unterschied liegt darin, daß » die Reihe der geraden statt der 
ungeraden ganzen Zahlen durchläuft. 

Die Randbedingungen (15) und (19) ergeben in ähnlicher Weise wie im vorigen Para- 
graphen ein Gleichungssystem von der Form (29) mit dem eben erwähnten Unterschiede, daß 
n eine gerade ganze Zahl ist. Die Bezeichnungen unterscheiden sieh nur dureh einige Vor- 
zeichen von denen in $ 2: 


(rk ym?’ 2) sin(ak%) 


a — / 92 a 
"Sm? cosh (ak 4°) —cos (ak 
m+-n—] 


(34) 
P„ behält die Bedeutung (26) unverändert bei. 

Gesucht wird der kleinste Wert von 4, für welchen das in Rhede stehende Gleiehungs- 
system der Form (29) von Null verschiedene Lösungen in den Unbekannten a ,d,,..., ba, bis... 
besitzt. Die numerischen Rechnungen wurden in derselben Weise wie im vorigen Paragraphen 
ausgeführt. Die Ergebnisse sind in Kapitel III zusammengestellt. 

Die allgemeine Lösung der Differentialgleiechung mit Berücksichtigung der Rand- 
bedingungen (17), (15) und (19) erhält man durch eine einfache Verallgemeinerung der Aus- 
(drücke (20) und (32). Diese zerfällt aber im Laufe der Rechnungen und gibt uns daher nichts 
uenes, das für unsere Zwecke zu berücksichtigen wäre, 


| 

| 

| 
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$ 4. Eine andere Art Reihenentwicklungen für die Lösung der Differentialgleichung. Um 
die Riehtigkeit der numerischen Ergebnisse zu prüfen, habe ich für eine quadratische Platte 
auch die folgende Form der Lösung  (, y) aufgestellt und durchgerechnet. Die Absicht war 
auch nachzuprüfen, welche Form den numerischen Rechnungen am besten angepaßt ist. 


cosh p, cosh p, 
+, b, eos | cosh _ cosh 
\y, 9, sinh 


Die Bezeichnungen sind dieselben wie früher. p,, pP, haben also die in (22) angegebene 
Bedeutung. 4% ist gleich eins, weil das Gebiet rechteckig ist. Die Gl. (21) nehmen daher 
die Form 

1=Vn’—2in, ,=yVn’+2in 


an. m und » laufen die Reihe aller ganzen positiven Zahlen durch. 
Jedes Glied erfüllt, wie ersichtlich, die Differentialgleichung. Jedes Glied der Reihen 
erfüllt die Randbedingungen 
für z=r wadfür ze—n, =0 für y und für y=-—n. 


Diese Form der Lösung ist also dem Bereiche 


nSeSn, "nsysa 
angepaßt. Dieses ist ein Quadrat mit der doppelten Seitenlänge in Vergleich mit dem Bereiche 
der vorigen Paragraphen, das ja für k=1 auch gleichseitig ist. 
ls ist keine Schwierigkeit, die Lösung auch für ein Rechteck mit ungleichen Seiten durch- 
zuführen. Die numerischen Rechnungen sind aber nur für ein gleichseitiges Bereich gemacht. 
Die noch übrigen Randbedingungen (x, + und w(+7,9y)=V ergeben nach ähn- 
liehen Rechnungen wie die in den früheren Paragraphen vorkommenden das Gleichungssystem 


An = Con = 0, 
m==] n=] 
Y 
C+ Im An +2 Cmn u, (nm 1,2, .) 
Y 


Die Bedeutung der Bezeichnungen ist 


113 2n? 

D,: 2 a4 cot(274))+ 15 je 

sin (274) +4 sinh (27 — 7°) 

Im / 2 ») / 2 12' 2 

ym’ cosh (27 ym? — 4°) cos (27%) 
[eothag, 
(_ 
— 44’ n? 


e, 4, und B,„ sind die zu bestimmenden Koeffizienten. 

Da der Bereich, wie erwähnt. hier die doppelte Seitenlänge hat wie in den vorigen 
Paragraphen, wird erwartungsgemäß der Wert von 4 halb so groß wie früher ausfallen. Die 
in Kapitel III angegebenen numerischen Ergebnisse bestätigen diese Erwartung. 


Als eine erste Annäherung können wir denjenigen Wert = 1,36 betrachten, für 
welchen der Ausdruck D, gleich Null wird. Die zweite Annäherung ist dann der Wert 
1,5817 , für den die Determinante D, des Gleichungssystems 
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sleich Null wird. Für die dritte Annäherung 4, 1,5875 verschwindet die Determinante D, 


des Systemes von fünf Gleichungen 
D,C —.c,, 4, — 6, Bı — As — B, 


+c0,4,+Pß, B.=D, 


Kapitel II. 
Die Knickfestigkeit rechteckiger Platten. 


$ 1. Bezeichnungen. Die rechteckige Platte soll durch die Geraden y=+b, = +kb 
begrenzt sein. Auf den Rändern y=+b wirkt je Längeneinheit ein Druck P,, welcher der 
y-Achse parallel ist. Der Druck P’, auf den Rändern «= + kb ist der #-Achse parallel. Be- 
zeichnet w die Ausweichung senkrecht zur Plattenebene, so hat die Differentialgleiehung des 
Gleichgewichts im kritischen Zustande der Stabilität die Form 


N ist der Plattenmodul 21) E ist der Elastizitätsmodul, o die Poissonsche 


Zahl, h die Dicke der Platte. 

Was die Ränder der Platte betrifft, können verschiedene Voraussetzungen gemacht 
werden. Wegen der in der Einleitung angeführten Gründen wurde das Beispiel durehgerechnet, 
daß alle vier Ränder eingeklemmt sind. In Formeln ausgedrückt: 


für =+kb muß w=ß0, =0, 
für y=t+b muß 
Neue Veränderliche werden entweder durch die Gleichungen 
€ n=7 ’ =, (37) 
b b n? N 1? N | 
oder durch die Gleichungen 


eingeführt. In den folgenden Paragraphen werden spezielle Voraussetzungen über die Kräfte 
P, und P, gemacht werden. Die obenstehende Differentialgleichung und die Randbedingungen 
reduzieren sich dann auf die schon untersuchten. 


$S 2. Der Druck ist auf alle vier Ränder gleichmäßig verteilt. Für eine quadratische 
Platte mit gleichmäßiger Belastung auf alle vier Ränder ist P, = P, und k=1. Bei der Be- 
rechnung der Kniekfestigkeit wird man auf das Eigenwertproblem des ersten Kapitels geführt. 
Als kleinster Eigenwert ergab sieh dort A = y 1,326 in naher Übereinstimmung mıt K. Sezawas 
und besonders G. I. Taylors Werten. Aus (37) findet man dann den Druck 


1809 N 


b? 


(39). 


Gefahr für Kniekung tritt für eine quadratische Platte ein, wenn sie 
dureh Kräfte zusammengedrückt wird, die aufalle vier Ränder gleich ver- 
teilt sind und den Wert P= 13,08 N/b? je Längeneinheit erreichen. 


S$S 3. Der Druck wirkt nur auf zwei Rändern einer rechteckigen Platte. Der Druck 7,, 
der in $ 1 erklärt wurde, wird gleich Null angenommen. Die Platte hat die Länge 2kb in 
(ler Richtung der x-Achse und die Breite 2b. Der Druck P,, den wir kürzer P benennen, 
ist der x-Achse parallel. 

Die Berechnung der Knickfestigkeit führt auf die Eigenwertprobleme des zweiten 
Kapitels. Nach (35) ist der Druck, bei dem Gefahr für Kniekung eintritt, gleich P- N, 


| 
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Um Annäherungswerte von 4 zu erhalten, haben wir diejenigen Werte berechnet, für 
welche gewisse s-reihige Determinanten, die in Kapitel II$ 2,3 und 4 erklärt wurden, gleich 
Null werden. 

Die zeitraubenden numerischen Rechnungen hat der Verfasser begonnen. Seine Frau, 
Elsa Faxen, und Architekt L. M. Giertz haben diese vollendet, wobei der letztere 
methodische Verbesserungen einführte. Aus den Gleichungen der Paragraphen 2 und 3 des 
Kapitels II wurden die in der folgenden Tafel angegebenen Annäherungswerte für — 
gefunden. 

In Kapitel II$ 2 wurde die in bezug auf = gerade Lösung des Problems gesucht. Das 
entsprieht bei Kniekung einer ungeraden Anzahl Ausbeulungen, oder wenn man Sezawas 
Ausdrucksweise den Vorzug gibt, einer ungeraden Anzahl halber Wellenlängen. Die Ergebnisse 
des Kapitels II$ 5 werden demgemäß in der Tafel durch die Worte „gerade Anzahl Aus- 
beulungen bei Kniekung”“ angedeutet. 

Ph: 


Annäherungswerte der Zahl = 
w Anzahl der Ausbeulungen bei Knickung 
gerade 
| | | | DB =6 
| 0.5 44.55 45.20 
0,75 28.1 
1.0 23.99 24.86 28,16 28,62 
| 7/6 23,71 
25 
| 1,25 22,77 24.55 
4/3 21.18 
20\ 1,5 19.82 22,74 19,79 20,49 
| | | 13 13.16 
\ 1.75 19,94 19.965 18,52 20,18 
| 1.8 14,40 
2,0 16,37 | 18,755 1943 | 1746 | 1998 


m 125 15 2m 23 


Die letzte Ziffer in jeder Zahl ist unsicher. Einen Stützpunkt für die Beurteilung der 
Genauiekeit dieser Zahlen kann man in der folgenden Reihe erhalten. Wendet man auch 
die Methoden des Paragraphen 4 Kapitel Il an, findet man für P6°/N die Zahlenreihe: aus D,: 18,3, 
aus D,: 23,99, aus D,: 24,69, aus D,: 24,56, aus D,: 24,57. Die obenstehende Tafel enthält 
die zweite und die vierte dieser Zahlen. 

Je größer die Anzahl s der Reihen der Determinante D, ist, je größer sind, wie man 
sieht, die Annäherungswerte. Man hat jedoch aus diesen Zahlen den Eindruck, daß die Reihe 
der Annäherungswerte rasch konvergiert. 

Für technische Zwecke dürfte in dieser Reihe der aus D, = 0 berechnete Annäherungs- 
wert 23,09 schon hinreichend genau sein. Insbesondere sollten die Werte genügen, die aus 
D, 0 berechnet werden. Leider reicht die Genauigkeit der Gleichung D,=0 aus nur für 
kleinere Werte des Verhältnisses % der Länge der Platte in der »-Riehtung zu der Länge in 
der y-Riehtung. Aus der obenstehenden Tafel sieht man nämlich, daß D, und D, sich mehr 
unterscheiden für erößere Werte von k als für # gleich Eins. Für ?=1,75 und k=2 kann 
man sorzar bezweifeln, ob die Gleichung D, — 0 hinreichend genaue Werte liefert. Daher wurde 
die Gleichung D, 0 für diese A-Werte gelöst. Das Ergebnis wurde in die Abbildung als 
stricehpunktierte Kurve eingezeichnet. 

Daß die immer besseren Annäherungswerte eine monoton wachsende Folge bilden, dürfte 
darin eine physikalische Erklärung finden, daß die steigende Annäherung einer immer festeren 
Bindung der Ränder entspricht, so daß die Steifheit der Platte dadureh erhöht wird. 

Gefahr für Kniekung tritt für eine rechteckige Platte ein, wenn sie 
dureh Kräfte zusammengedrückt wird, die auf zwei entgegengesetzte 
Ränder wirken und die den Wert P erreiehen. Für P erhält man einen An- 
näherungswert aus obenstehender Tafel, welehe die Annäherungswerte 
des Bruches PW/N enthält. Die Werte aus der Gleichung D,= 0 sind genauer als die 
aus D, 


_ lech. 
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Mehr anschaulich als durch die Tafel wird die Kniekgefahr durch die Kurven der vor- 
stehenden Figur dargestellt. Die gestrichelten Kurven bedeuten Kniekung durch eine ungerade 
Anzahl Ausbeulungen. Die ausgezogenen Kurven bedeuten Knickgefahr durch eine gerade 
Anzahl Ausbeulungen. 

Dünne Kurven bedeuten das Ergebnis der Gleichung D,=0. Die starken Kurven be- 
deuten das Ergebnis der Gleichung D, =0, sind also genauer. Der strichpunktierte Kurven- 
zug wurde aus der Gleichung D,=0 ausgerechnet. Dieser ist also noch genauer. In der 
Abbildung wird aber diese Verbesserung unmerklich, wenn k <1,75 ist. 

Die dünne gestrichelte Kurve geht von links bis k=y 5. Sie bedeutet eine Ausbeulung 
und endet dort. Der nae :hfolgende gende entspricht Gefahr für Kniekung 
dureh drei Ausbeulungen. 

Die starke gestrichelte Kurve fällt von 70,50 an anfangs stark, dann sehr langsam. 
Das Ende bei 1,25 rührt daher, daß die Kurve nicht weiterverfolgt wurde. Dieser Kurvenzug 
bedeutet Kniekung durch nur eine Ausbeulung. 

Die wirkliche Kniekgefahr erhält man in unserer besten Annäherung, 
wenn man die Ordinate Pb?N für diejenigen der starken Kurven abliest, 
welehe für den betreffenden Wert von %k den kleinsten Wert gibt. Wenn 
» 1,75 ist, gibt jedoch die stricehpunktierte Kurve die beste Annäherung. 

Bei einem kleinen Wert von k findet man auf diese Weise, daß eine einzige Ausbeulung 
in bezug auf Kniekung am gefährlichsten ist (starke gestriche ie Kurve). Wird % größer als 
1.15, werden zwei Ausbeulungen gefährlicher, denn die ausgezogene starke Kurve schneidet 
die starke gestrichelte und geht dann unterhalb dieser. Die Mittellinie der Platte hat die 
Form eines länglichen s. 

Bei noch größeren Werten von k, d. h. noch mehr länglichen Platten,. schneidet die 
starke gestrichelte Kurve nochmals und fällt unterhalb der starken ausgezogenen. So ver- 
läuft auch die Verbesserung, die strichpunktierte Kurve. Das bedeutet Knickgefahr bei drei 
Ausbeulungen. 

Unsere starken Kurven stimmen mit Sezawas Werten recht gut, wie man aus der 
folgenden auf unsere Bezeichnungen umgerechnete Tafel ersieht: 


Zahl 

Wellen- 1.0 1.5 

längen 
50,7 26,3 32.9 
2 29.1 20.2 17.1 
225 14.9 
4 22.5 


Diese Tafel stellt nur eine Art erster Annäherung dar. Darin dürfte die Erklärung der 
weniger guten Übereinstimmung bei k=2 liegen. 

H. Grauers°’) behauptet, die genaue Lösung der Schwingungs- und Knie kungsproble me 
der rechteckigen Platten gefunden zu haben. Seine Darstellung leidet aber an inneren Wider- 
sprüchen, auf die wir hier aus Platzmangel nieht eingehen können. Die numerischen Ergeb- 
nisse stimmen, soweit geprüft wurde, auch nicht. 

Wir möchten endlich bemerken, daß die dureh (25), (26), (33) erklärten «a„ und pP, 
eine selbständige Bedeutung haben. «a, 0 bestimmt die Knickbelastung, wenn nur die mit 
der „Achse parallelen Ränder der Platte eingeklemmt sind und die übrigen zwei gestützt. 
Aus 5,0 kann man dagegen die Kniekbelastung berechnen für eine Platte, deren mit der 
x-Achse parallelen Ränder gestützt und die übrigen eingeklemmt sind. Die Kraft ist der 
Achse parallel und wirkt auf die mit der y-Achse parallelen Ränder. 

Der Verfasser hat die Hoffnung, u. a. die Rechnungen auch für den Fall durehzuführen, in 
dem die Befestigung der Platte keinerlei Ausdehnung nach den nicht belasteten Rändern hin 
erlaubt. 

Die Verwendbarkeit der obenstehenden Methoden auf manche schwierige Randwert- 
probleme bei den Differentialgleichungen Ju —+...=0 oder I + 0) möge am Ende er- 
wähnt werden. 324 


(söteborg, im Mai 1934. 


5), Ingenjörsvetenskapsakademiens handlingar N:r 98, Stockholm 1929. 
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Das Ausknicken eines Plattenstreifens unter Schub- und 
Druckkräften. 


Von €. Schmieden in Rostock. 


Einleitung. Die vorliegende Arbeit hat das Ziel, die Knicklasten eines unendlich-langen 
Plattenstreifens bei gleichzeitiger Beanspruchung durch Schub- und Druckkräfte in seiner 
Ebene zu bereehnen. Die zur Lösung führende Methode wurde erstmalig von Southwell 
und Skan') angewandt und seitdem in einer ganzen Reihe von Arbeiten verschiedener 
Autoren auf ähnlich gelagerte Probleme übertragen. Methodisch bringt die vorliegende Arbeit 
daher nichts wesentlich Neues; dagegen dürften ihre Resultate einiges praktische Interesse 
bieten, da sich ganz außerordentlich einfache Schlußformeln ergeben, die mit einem Minimum 
von Arbeitsaufwand die Knicklast für jede beliebige Beanspruchung zu berechnen gestatten. 


1. Differentialgleichung und Lösungsansatz. Der Plattenstreifen erstrecke sich in „-Richtung 
von — so bis + und besitze die Breite 2b, so daß also „== + b die Gleichung seiner Ränder 
ist. Der Streifen werde in seiner Ebene am Rande durch einen konstanten Schub S/em und 
einen konstanten Druck P,/em in y-Riehtung, ferner durch einen konstanten Druck P,/em 
in „-Riehtung beansprucht. Unter diesen Voraussetzungen lautet die Differentialgleichung für 
die Durechbiegung w 


\ 


N ist dabei die Plattensteifigkeit, 2h die Plattendieke und » die Poissonsche Konstante. 


Durch Einführung von 


wird (1) dimensionslos und geht über in 


Randbedingungen sind 


für &: »w muß in & eine rein-periodische Funktion sein; 


für „7: a) eingespannter Rand: w = Mr 
b) frei-aufliegender Rand: =0 für n=tl. 


Das Wesen der Methode von Southwell und Skan besteht nun darin, daß nicht w 
direkt, sondern eine komplexe Funktion w’ bestimmt wird, die der Differentialgleichung 
genügt (#9’ kann natürlich nieht analytisch sein). Da (la) eine lineare Differentialgleichung 
ist, muß dann sowohl der Realteil von »»’ wie auch der Imaginärteil für sich dieser Differential- 
eleichung genügen. Als Durchbiegung w kann man dann entweder den Realteil oder den 
Imaginärteil von »’ benutzen, da beide bis auf eine Parallelverschiebung des Koordinaten- 
anfangs in Richtung dasselbe geben. 


Für ein Partikularintegral »»’ machen wir jetzt mit festem, aber noch freiem z den Ansatz 

der beim Einsetzen in (la) für A die Gleichung vierten Grades mit reellen Koeffizienten liefert: 
"+22? . . 


/u einem x gehören vier Wurzeln 4/;. so daß man vier unabhängige Partikularintegrale und 
damit, bei gegebenem x, das vollständige Integral erhält, das die Randbedingung für & erfüllt. 


!) Proc. Royal Soe. Ser. A. Bd. 105 (1924), S. 582. 
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Da der Koeffizient von 4° in (3) verschwindet, ist die Summe der vier Wurzeln Null, 
so daß man ansetzen kann 


Hier ist @ immer reell, $ und y sind, unabhängig voneinander, entweder reell oder rein 
imaginär, je nachdem die Wurzeln von (3) alle reell, oder zwei reell und zwei konjugiert 
komplex, oder je zwei konjugiert komplex ausfallen. Auf Grund von (4) kann man durch 
Koeffizientenvergleich (3) ersetzen durch 


2. Diskussion der Knickgleichung. Setzen wir die Randbedingungen für mw’ an, so er- 
halten wir in beiden Fällen vier homogene lineare Gleichungen für die vier Integrations- 
konstanten, die nur dann eine von Null verschiedene Lösung besitzen, wenn die Determinante 
der Koeffizienten verschwindet. Diese Determinante, gleich Null gesetzt, ist die Kniekgleichung. 


Sind alle Wurzeln 4; verschieden voneinander, so erhält man 


a) eingespannt: 2 — eos ta] = [4a® — — y*|sin 25 sin 2y, 


| (m). 
b) frei-aufliegend: Sa? [eos 2p cos 2y ta] = [H a? y?) y’P]sin2/p sin 2y 
Sind zwei Wurzeln gleich (y==0), so erhält man auf dieselbe Weise 
a) eingespannt:  P(cos2P — costa) = P?)sin2P, 
dm). 
b) frei-aufliegend: (cos 2P - costa) = — p?)sin2p oder P=V 
Für zwei Doppelwurzeln (die notwendig reell sein müssen, s. u.) (= y=0) 
a) eingespannt: sın’Z2a, 
| u”). 
b) frei-aufliegend: da’(l costa)=V 
Und schließlich für eine dreifache und eine einfache Wurzel (y 2a) 
a) eingespannt: 4dacosta=sinta, 
(11°). 


b) frei-aufliegend: at acos4a+>3sin da) 
Eine vierfache Wurzel tritt nur für «0 auf, also in dem trivialen Fall des Nicht- 
auskniekens. 


Bei der Diskussion der Gl. (II) können wir uns auf den Bereich «>0 und p,y positiv 
reell oder positiv imaginär beschränken, da sich bei Vorzeichenumkehr einer dieser Größen 
(die Gleiehungen nicht ändern. Dann besitzen die Gl. (II) zunächst die trivialen Lösungen 


2. 3. P+y=2a, 4. 


—9da. 


Alle diese Lösungen sind unbrauchbar, da jedesmal zwei der 4; gleich werden, also Doppel- 
wurzeln auftreten, für die (Ib) nieht gilt. Ebenso ist es bei den Gl. (11); hier sind =0 
und 5==2a triviale Lösungen, denen der Fall zweier Doppelwurzeln bzw. der einer dreifachen 
Wurzel entsprieht. (11°) besitzt im Falle a) nur die Lösung @«==0, im Falle b) außerdenı 


NT 


die Wurzeln «a - . Für einen frer-auflieeenden Streifen können also mörlicherweise zwei 


Doppelwurzeln auftreten. Schließlich kann eine dreifache Wurzel für beide Arten der Auf- 
lagerung auftreten. 

Wir beweisen jetzt zunächst, daß die Kniekgleichung im Bereich kleiner «a, p, y keine 
Lösung besitzt. Entwickeln wir (II) in der Umgebung des Nullpunktes, so erhalten wir in 
beiden Fällen bis auf Terme achter Ordnung 


—- ka — 
so daß für kleine a, ß, y nur die auszuschließenden Lösungen auftreten 


5 
| 
— ; | 
| 
| 
. (8). 
»rfüllt. 
| B=0; y=0; 2a=+t(ß+y); 2a=+(ß—y). 
4 | | | 
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Weiter beweisen wir, daß nicht alle 4, komplex sein können, oder mit anderen Worten, 
daß $ und y nicht gleichzeitig imaginär sein können. Setzen wir =ip" und y=iy’ und 
lassen die Akzente gleich wieder fort, so geht (lla) über in 


2 By[&of 2P 2y —cos4a]=[4a?+ P? y?] Sin Ein 2y, a,ß,y>VO 
oder 
2 2% 2y - 1+2sin2a]l= [4a + 
Nun eilt zunächst 
t Sin2PSin2y, 


da für 


sin’ 2a Zin2Pp Sin2y 


ls bleibt zu zeigen, daß auch 


2 By[Eoj Eof2y—1]< (A? Sin2P Sin2y. 


‘ 
Zum Beweise dieser Ungleichung leiten wir aus ihr die folgenden Ungleichungen ab: 


2 By [Eof 2 y) + —y) -2]< + [Eof + y) — — y)] 


(B+ < Sin’ 7) 


Sin?’(d—y) _ Sin’(P-+y) 
Diese letzte Ungleichung ist aber sicher erfüllt, da beide >O und mit wachsendem 
monoton zunimmt. 
Ebenso läßt sich zeigen, daß beim Auftreten einer Doppelwurzel die beiden andern 
Wurzeln notwendig reell sein müssen. 
Im Falle (IIb) schließlich sieht man unmittelbar, daß P und y nicht beide imaginär 
sein können, da hier beide Seiten der Gleichung immer entgegengesetzte Vorzeichen bekommen. 
Im folgenden setzen wir demnach immer 5 als reell voraus, während y sowohl reell 
wie imaginär werden kann. 


3. Reiner Druck in Längs- und Querrichtung (7 0). Ist der Schub nicht vorhanden, so 
folgt aus (die andere Möglichkeit 5 y scheidet aus, wie sich im folgenden ergibt). 
Damit vereinfachen sieh (ID) und (ID) erheblieh zu 


2 . 2.2.2 2... Alla), 


Aus (1) folgt, daß > im ganzen Bereich imaginär sein muß: für o,=0 ist das klar. 
Würde für ein endliehes o, y reell, so müßte es ein o, geben, für das y=t ist. Für y=0 
eilt aber (11’), das jetzt lautet: 


uy eehört jedenfalls z=0, also o, yeider ällen die 
/u gehört jedenfall al bzw l. h. ın beiden Fällen d 
Kulerlast für reinen Druck in Querrichtung. Die andere Möglichkeit ==? liefert für 
0,9 +22°> 9°, scheidet also aus. Somit ist y im ganzen zulässigen Bereich rein imaginär. 
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Für den frei-aufliegenden Streifen kann man jetzt die Lösung sofort hinschreiben: 


‚nn 7 
aus (I1b) folgt also = und hieraus 
n° 4 
also x N 4 


/u jedem 
eehören zwei Werte von z, wobei für z>0 und z>x beide Male 0,>x geht. Unter allen 
o, gibt es ein kleinstes, das sich leicht bestimmen läßt. Man erhält 


erhält man also ein kontinuierliches Eigenwertspektrum: zu jedem ©, 


Interessant ist, daß man für ein o,, das beliebig wenig kleiner ist als —, die Eulerlast, zum 


Ausknieken ein > braucht. 


Nicht ganz so einfach erledigt sieh der Fall eingespannter Ränder. Die Gl. (Ila) schreibt 


sich wegen des imaginären y (y=iy’ und der Akzent wieder fortgelassen) 


2 ßBy(1— cos2P Euvj2y) = (P? — y?)sin2P&in2y. 
Wie man leicht nachrechnet, läßt sich dieser Ausdruck zerlegen in 
cos Ziny — ysinp &ofy) (P + yeosp Ziny) 0. 


Zusammengehörige Werte von p,y folgen also aus den beiden Gleichungen 


‘ 


ptgPp: oder petgeß=yltgy. 


Bei festem y liegen die kleinsten Werte von 5 bei der ersten Gleichung in dem Bereich 


> 
„ <P<a, bei der zweiten in dem Bereich a <P<5y7, also ist, um die kleinsten Knick- 


) 
lasten zu erhalten, die erste Gleichung zu benutzen. Sind 5 und y bestimmt, so erhält man 
bei gezebenem o, und aus (]), 


Die numerischen Ergebnisse sind in Tabelle I zusammengestellt. Die unterstrichenen 
Werte gehören zu den numerisch sehr genau bestimmten Minima (für 0,95 und 0, 9,8578 


nicht berechnet). Für o,=7° erhält man als Minimalwert und 0. Das ist 
+) 
2 
wieder in dem Sinne zu verstehen, daß für beliebig wenig kleiner als 0x > 
muß, um ein Ausknieken der Platte zu erreichen. 
Tabelle]. 
% x % x Y x % oO x x x 
2,616 38.36 3.894 22,53 2.0695 15,12 1.805 10.51 1.666 163 1,519 8,78 
3.696 21,84 3.009 17.50 2,278 13.79 1.50 10.05 1.159 914 1.315 8,23 
2,749 17.67 2.400 15.66 1.994 13,32 1.571 3.85 1.243 S,80 1.088 7,13 
2.319 17.20 2,174 15.19 1.902 19.51 1.514 1.003 s.71 0.825 1,27 
2.064 17.65 1.921 17.78 1.627 13.09 1.135 9.95 1,706 0,437 6,96 
1.480) 22,97 1,497 18,2% 1.330 15.15 0.869 11.01 


| | 
| 
| 
| 
| IM] 
| | 
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4. Druck in Querrichtung und Schub (os, 0). Es ist klar, daß man sich auf Werte von 
2 


a . 
(eingespannt) und 


(frei-aufliegend) beschränken kann. Die Gl. (I) lauten jetzt 


2a — 2x? — 0,; (a? — P?) (a — y?)=x°; a(® —y’)=xt. 


Für reinen Schub (o,° 0) haben Southwell und Skan die Lösung angegeben, y wird rein 
imaginär. Für verschwindenden Schub wird aber y reell, so daß bei kombinierter Beanspruchung 
für gewisse o, und 7 y verschwinden, also Doppelwurzeln auftreten müssen. Für diese Doppel- 
wurzeln muß außerdem a’ > und o,< 05 sein. Die Auflösung von (11’) ergibt, daß bei ein- 

gespanntem Rand Doppelwurzeln für #,> 9,15 auftreten, und bei frei-aufliegendem Rand für 
6,» 2,145. In beiden Fällen ist der Bereich, in dem Doppelwurzeln auftreten können, sehr 
klein, da die zugehörigen Werte von o, sich nur wenig von den größtmöglichen o, unterscheiden. 


Das numerische Material wurde so gewonnen, daß o, ein fester Wert <oy, erteilt wurde 
und die zugehörigen Werte von z und r aus (Il) und (Ila) bzw. (IIb) numerisch berechnet 
wurden. Die Ergebnisse sind in der Tabelle II zusammengestellt, wobei das Zahlenmaterial 
für o, VO der Arbeit von Southwell-Skan entnommen wurde. 


Tabelle II. 


a) Eingespannte Ränder: 0x, =. 


oy—b 
a % T 3 % T a y % T 
2.079 1.787 4,70: 23,07] 2.008 1,802 4,0395. 2 19.63 | 2,036 1.800.3,6787 2,0. 17,07 
1.977 1.804 4,3340 2218| 1,916 1.816 3.5463 19.36 | 1,920. 1.815 3,0285 15 12.9 
1.003 1.817 3.934 2275| 1.850 | 1,829 3.0315 1 23.18| 1.850 1,822 24885 1.0 17,54 
1.858 | 1,826 3,6087) 25,78 1.823 1.821 2,1135 0,5 | 28,00 
oy=95 0, 9.8578 
1.938 | 1.803 125025 1,5 | 12.29] 1.940 | 1.798 | 2,4005 1.5 | 11,66| 1,830 | 1,778 1,415{ | 1,0 2,45 
1.840 | 1.797 1.7325 1,0: 11.46| 1,835. 1.787 1558: 1.0 | 10,32| 1.748 1.735 8,64 
1.775 | 1.771 05 1445| 1,752 1.746 0,428i 05 11,32] 1.699 | 1,695 0.944 | 0,4 8,42 
1.734 1.731 0.4105 04 12,26] 1,646 | 1.645 1285 | 02 8,69 


b) Frei-aufliegende Ränder: 


Oy 
3 % a y % T 
2205 1.479 15.6%6i| 11531 | 1395 3214i| 15 292,58 
1.S50 1.452 4.479 17.42 | 1.440 | 1.367 12.846i 12 | 11.96 
1.556 1.410 3.4300 13,63 | 1.388 | 1,347 2,6217 1.0 | 12,06 
1445 1.383 | 2,9770 13.21 | 1.2096 1.292 2175i | 05 1659 
1.366 1.370 2.6091 1.97 
1.6 2,4 
1.529 1.381 3.0705 1.5 11.61 11.526 1.366 2.9365 | 1,5 | 10,67 
1.430 1.349 2,6825 1.2 10,74 | 1.415 1.323 | 2497i| 12 9.412 
1.370 1.322 24295 1.0 1048 11.345 | 1.288 2.2095 | 1.0 8.79 
1.248 1243 1.8877. 05 1274 [1.159 1.151 1453i| 05 7.97 


1.100.1,099 1.1865 0,3 


| 
| 
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Aus der Tabelle sieht man, daß bei eingespannten Rändern a, ß, y desto mehr gegen 


schen, je mehr sich o, a? nähert und je kleiner z ist. An Stelle dieses Wertetripels kann 
man auch das Tripel a,y klein, 5 annähernd = einführen, was für die 4, dieselben Werte 
ergibt. Für die folgende Rechnung ist dies letztere Tripel angenehmer zu behandeln. 


Wir verschaffen uns jetzt die Nachbarlösungen zu indem 
wir setzen 0, =? —0',, B=a+p, a=a, y=y', Alle gestrichenen Größen sollen 
dabei als klein betrachtet werden. 


Für den Fall eingespannter Ränder, den wir allein ausführlich durchreehnen wollen, 
hekommen wir beim Einsetzen in (lla), wenn wir die Akzente gleich wieder fortlassen , 


2 (a+Bß)yl[eos2 —cos4a]= [4 a’ — — — y’]sin2 2y 
2°? 
(a n?—2nß - P)(a — a ("+22 — 
Eintwiekeln wir die erste Gleichung, so bekommen wir bis auf Terme fünfter Ordnung nach 
Kürzen durch 4y 


In erster Näherung folgt pP = y?- 4a?, p ist also klein zweiter Ordnung. Setzen wir dem- 
gemäß P 0°, so wird bis auf Terme sechster Ordnung 


a (y? —4 0?) 16a — 30) +5 —4da)= 


| 
In den Termen vierter Ordnune können wir 0° -—- (v? - 4a?) setzen und erhalten somit ın 


derselben Annäherung 


| 3 


während die beiden anderen Gleichungen bei demselben Grad der Näherung liefern 


2 +) A 

6a? —3y? —(y? — 4a?) - 2x’-+0,, a’ — y? 


Die Elimination von a,y aus diesen Gleichungen ergibt zwischen 7 und x die folgende Be- 
2 
17 
ziehung | "er 0,002) 


Diese Näherungsformel gilt nach ihrer Ableitung für kleine o, und z. Sie gibt für 


/ 


und 2>0, T>x, aber so, daß also klein erster Ordnung ist, wenn ©, 


als klein zweiter Ordnung vorausgesetzt wird, wie es aus den obigen Gleichungen folgt. Der 
Koeffizient von 7? ist im ganzen zulässigen Bereich so nahe Eins, daß man ihn genau genug 
(durch Eins ersetzen kann, also 


2 
=— +72°x%° 
3 3% 
und für 
= 


| 
15.95 
| 
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Auch in diesem Falle erhalten wir also wieder das Resultat, daß zu einem Querdruck, 


2 
der beliebig wenig kleiner als 7° ist, ein endlicher Schub r > V „ nötig ist, um den Streifen 


zum Auskniecken zu bringen. Die Erklärung für dieses zunächst überraschende Verhalten der 
Platte ist in beiden Fällen dieselbe: Wirkt ein Querdruck allein, so kniekt die Platte als 
(Ganzes ohne Wellenbildung aus, bei kombinierter Beanspruchung tritt dagegen beim Ausknicken 
immer eine Wellenbildung ein, und die potentiellen Energien in beiden Zuständen unterscheiden 
sich um einen endlichen Betrag. 


5. Druck in Längs- und Querrichtung mit zusätzlichem kleinen Schub. In derselben Weise 
wie am Ende des vorigen Absehnitts verschaffen wir uns jetzt die Nachbarlösungen zu dem 
Falle reiner Druckbeanspruchung, setzen also r als klein gegen o, und o, voraus. Weiter 
wollen wir o, fest vorgeben und nach dem Zusammenhang zwischen 7 und o, fragen. Aus- 
führlich behandeln wir nur den Fall frei-aufliegender Ränder. Die Minimalwerte von 0,.,%,P, y 
als Funktion von o, bezeichnen wir mit Oxg Zu Po, y. und entnehmen dem Abschn. 2 


Po 5 ‚vo 4 | 4 Au, .) | 4 Gy; = ») +2 


Aus der Gl. (IIb) des Abschn. 2 folgt zunächst wegen 9, = wenn wir wieder 


und y=y,+y’ setzen (a klein erster Ordnung = «'’) 
(y? + PB) Sin2y, 


oder wegen P = 0x0 0y 


Die Gl. (I) liefern, wenn wir 0x, 0x: 2° 2, + setzen und gleich die endlichen Terme 
fortlassen 


3a’+2yy, -P 


Die Elimination von a’, pP’, y’ aus diesen Gleichungen gibt eine Beziehung zwischen 7 
und o’,., da sieh die mit #° behafteten Glieder dabei herauskürzen, wie es sein muß, da wir 
an der Stelle des Minimums von o, als Funktion von x entwickelt haben, z sich also in erster 
Näherung nieht ändern darf. Diese Gleichung zwischen 7 und 0’, lautet 


M 

y 
( Oxo 


Speziell folgt für 0 und damit 


Ta 


ay» 
v3” 2 
In dieser ersten Annäherung bleibt z unverändert. Der hinzutretende Schub bewirkt 
eben in erster Annäherung nur eine Verzerrung des Knickfeldes, eine Änderung der Wellen- 
länge dagegen erst in zweiter Ordnung. 


“ur Io OÖ 5; Y uni . 
xo ‚vo ‚iv x 4 (32 3 7”) 


Diese spezielle Lösung gilt streng (vel. das Ende von Absehn. 3). 


2a?’ y,n*Cotgy, 
. 
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Eine entsprechende Näherungslösung läßt sich auch bei eingespannten Rändern gewinnen. 
Ihre Ableitung ist allerdings wesentlich komplizierter, da man die Minimalwerte von o, nicht 
in geschlossener Form angeben kann, sondern nur ihre Zahlenwerte. Prinzipiell verläuft die 
Rechnung aber genau wie bei frei-aufliegenden Rändern, so daß wir auf eine Wiedergabe hier 
verziehten können. In beiden Fällen erhält man eine Beziehung von der Form 


"= — 0); 


webei 5, jetzt den wirklichen Druck in Längsrichtung bedeutet. Der Faktor A hängt nur 
von o, ab. In der Tabelle III ist er für eine Reihe von Werten von o, eingetragen. Mit 


Tabelle 


a) eingespannte Ränder b) frei-aufliegende Ränder 
0 128,802 17,200 | 2,379 | 22,26 | 2,00 | 22,0 0 4,2127 9,8696 | 1,5708 13.23 | 1,42 | 13,00 


3024251 15486 2,174 1938 1.74 19.2 | 0,8 4.0047 8,9916. 1,4242 123,01 1.25 | 11,95 
619468 13,313 1.902 16.10 154 159 | 1.6 3.7502 7,8606 1.2095 10,51, 1.02 | 10,48 
9% 13,491 9813| 1,314 11,50 1.14 11.3 | 24 |3,3294 5,7504 0,6386 | 7,98 | 0,50 | 7.97 


a? 9,870 6,580 0 0 8.061 — 3.1912 4,0318 0 7,09 


7" und z#* sind die Werte von r und 2 bezeichnet, die sieh auf Grund der Näherungsformel 
für 6,=0, also für verschwindenden Längsdruck ergeben. In den Nachbarkolonnen sind die 
Werte von r und x eingetragen, die sich auf Grund einer Interpolation aus den im Abschn. 3 
berechneten Werten als kleinster Schub mit zugehörigem z ergeben. Diese Werte sind 
natürlich auch noch mit den unvermeidlichen Interpolationsfehlern behaftet, die besonders bei 
den Werten von z ziemlich erheblich ausfallen dürften. Die Übereinstimmung der beiden 
Werte des kritischen Schubes ist so gut, wie man es nur wünschen kann, die Übereinstimmung 
der #-Werte bei Berücksichtigung der Interpolationsverfälschung, die bei den #-Werten wohl 
um eine Größenordnung höher ist als bei den r-Werten (das Minimum verläuft außerordentlich 
ach) ist auch noch erträglich. 


Aus diesem Ergebnis können wir schließen, daß unsere Näherungsformel im ganzen 
Gebiet der kombinierten Druck- und Schub-Beanspruchung die kritischen Lasten mit praktisch 
mehr wie ausreichender Genauigkeit zu berechnen gestattet, da sie selbst in dem äußersten 
Grenzfall reiner Schubbeanspruchung noch die kritische Last innerhalb einer Fehlergrenze 
von 2°/o liefert. Um ganz sicher zu gehen, wurde übrigens auch noch der Fall Längsdruck 
und Schub, also o,=0, streng durchgerechnet und die Resultate mit den Werten auf Grund 
der Näherungsformel verglichen. Das Resultat entsprach den Erwartungen: die Abweichungen 
nahmen mit wachsendem 7 allmählich zu, blieben aber immer unter der obigen Fehlergrenze 
von 2 °Jo. 


Die Knicklasten des Plattenstreifens bilden eine untere Grenze für Platten endlicher 
Länge, deren Längsseiten ebenso gelagert sind wie der Streifen. Die wirklichen Knicklasten 
solcher Platten müssen natürlich größer sein, immerhin dürften aber die Abweichungen beı 
nicht zu ungünstigem Seitenverhältnis nicht allzu groß sein, so daß unsere Knickformel auch 
hier brauchbare Resultate geben dürfte. 536 


\ 
| 
| 
y 2 
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Zur Berechnung des Reibungsschwingungsdämpfers. 
Von Fritz Söchting in Wien. 


ur Berechnung des Reibungssch wingungsdämpfers liegt das Verfahren von J. Ormondroyd') 
/ bzw. (1. Jendrassık?) vor. Bei diesem Verfahren wird das Reibungsglied konstant und 
die Schwingungsbewegung sinusförmig angenommen. Sie untersuchen nur die Dampfung auf 
der Kraftseite und erhalten um 40% höhere Werte für das Dämpfungsmoment, als theoretisch 
erforderlich ist. 

Das folgende Verfahren ermittelt für ein Zweimassensystem ohne diese Vereinfachungen 
bei gegebener Massenverteilung das Reibungsmoment für den Kleinstwert der Ausschläge und 
damit der Beanspruchungen. 

Gegeben sei ein Zweimassensvstem mit den Trägheitsmomenten I,, I,, der Feder- 
konstante e und dem auf 7, wirkenden erregenden Moment Psino#. Wenn der Dämpfer | 
mit dem Trägheitsmoment /, durch das Reibungsmoment M mit dem Zweimassensystem ge- | 
koppelt ist, so sei » die Eigenfrequenz der neuen Massenanordnung. | 

l. Bestimmung des Dämpfungsmomentes. 

ls sind zwei Fälle möglich, je nachdem ob der Dämpfer auf der Kraftseite (der Stelle I, 
vel. Abb. la) oder auf der Lastseite (Stelle Il, vel. Abb. 1b) wirkt. 

a) Der Dämpfer wirkt auf der Kraftseite. 


Fall a Fall b 


Für diesen Fall lauten die Differentialgleichungen der Bewegung 


FM=V0 | | 


Aus diesem Gleiehungssystem (1) erhält man durch Differentiieren 

=0 | 


Wie man aus Gl. (2) ersieht, fällt die kritische Drehzahl der gedämpften Schwingung 
mit der der ungedämpften zusammen °). 


4 
| 
Aus Gl. (2) ist 9, eliminierbar, und unter Berücksichtigang von Gl. (3) erhält man 
Zu 
und nach zweimaligen Differentiieren 
( (Vi m? == cos en . 
12 


Aus diesen beiden Gleiehungen kann man das Kosinusglied eliminieren 
Die eharakteristische Gleichung lautet dazu 


') J). Ormondreoyd: Frietion Dampers and their Applieation to Engines. Verhandlungen des 3. Internationalen 
Kongresses für Technische Mechanik in Stoekholm, III. Teil, S. 221, Stoekholm 193», 

?\,G. Jendrassik: Theorie des Reibungsschwingungsdämpfers, VDT 1953, S. 1009, 

°»), MH. Frahm: Neue Untersuehungen über die dynamischen Vorgänge in den Wellenleitungen von Sehiffs- 
maschinen mit besonderer Berücksichtigung der Resonanzschwingungen, VDI 1902, S. 797 und S. 880. 
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Daher lautet die Lösung für 9, 
Die Konstanten €,, F, und H, bestimmen sieh aus der Differentialgleichung (4) zu 


Es ist also 
Mt: Pt 
2(1,+1 
Ist 9, bekannt, so bestimmt sich 9, aus einer der Gl. (1), und y, wird durch Integration 
aus der ersten Gl. (1) gefunden. Von Interesse ist aber nur die Differenz £ der beiden Aus- 
schläge 9, und 9,, von der die Beanspruchung der Welle abhängt. Es ist nach Gl. (1) 


9,=4A,+B,t 3 ‚+ coswt Ba). 
2 


Pont 
 \ 2(1,- 


(6). 


Die Werte der Konstanten E, und @, bestimmen sich aus den vorgegebenen Anfangs- 
bedingungen und gelten nur solange, als das Reibungsmoment sein Vorzeichen nieht wechselt. 
Der Vorzeichenwechsel des Reibungsmomentes fällt mit dem der Relativgeschwindigkeit { 
der Dämpfermasse /, gegenüber der Masse I, 

vo 3 
zusammen. Es sei nun angenommen, daß der Vorzeichenwechsel von Z immer nach wf= a 


stattfindet), daher gelten die Werte der Konstanten in (6) nur während einer Zeitdauer 
7 


i=— und sind dann jedesmal von neuem zu bestimmen. Es ist bei f, = 0 
eh 


und bei of, 7 
I, 


2 I, 2 


Ira). 


- 


Sollen nun die Beanspruchuneen und damit die Relativverdrehungen nieht zu groß 
g 
werden, dann darf für den Grenzfall, bei welchem die Ausschläge gerade noch endlich bleiben, 


E, &, sein und aus den Gl. (7) und (7a) erhält man für das Reibungsinoment den Wert 
JT 
M: 


(der ein Zugroßwerden der Beanspruchung verhindert, die Ausschläge bleiben endlich, trotz- 
(dem die Frequenz der erregenden Kraft gleich der der Eigenfrequenz ıst. 


4) Diese Annahme ist nieht nötig. Tatsächlich ist nach der ersten Gl. (1) 


412 
und nach Gl. (1) und Gl. (a) 
Mt: „ Plz lila.” Pt 
damit wird 
Yı Bs e(h +19)? | coswt A 3ollı +1) (ba) 


Die Konstanten der Gl. (6a) bestimmen sieh aus den vorgegebenen Anfangsbedingungen. Es sei bei fo 0, So. ®, 
dann lautet Gl. (sa) 


1 


Der nächste relative Stillstand, bis zu welehem die Werte der Gl. (6), Gl. (6b) gelten, ist aus Gl. (6b) zu er 
mitteln. Für das nächste Intervall sind nun zufolge des Zeichenweehsels von M die Werte der Konstanten neuerlich 
zu bestimmen. 

Will man untersuchen, wie groß Es, Ga und die Zeitdauer von einem relativen Stillstand bis zum nächsten, bei 
gegebenen Anfangsverhältnissen nach sehr langer Zeit werden, dann müßte man den Grenzwert bestimmen, der sich 
nur sehr umständlieh ermitteln läßt. Diese Grenzwertbestimmung kann umgangen werden, wenn die Anfangswerte 
so gewählt werden, daß der Vorzeichenweehsel immer nach einer Zeitdauer or=7 stattfindet. Bei der besonderen An- 


nahme bei or— 7 wird 


os 
(1—eosot) 4 I, 30 1») 


Es (1— + \osinot (6b). 


n 
e (I, + Ia)? 


| 

| 
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Der Wert für M nach Gl. (5) bedingt die kleinste Rtelativverdrehung; denn für jeden 
Wert, der kleiner als 7 P ist, wachsen die Größtausschläge nach Gl. (7) und (7a) immer 


mehr und mehr an. Aber auch für jeden größeren Wert von M sind die Relativverdrehungen 
erößer; denn ist die Frequenz k=++@ und ist £{=4, 5,0, das heißt, die Bewegung von I, 
ist die von /,, so wird bei Vergrößerung des Koppelmomentes M auch der Ausschlag 9, ver- 
erößert, bei welchem der Dämpfer zu wirken beginnt, weil bei der ungedämpften Schwingungs- 
bewegung 5, proportional dem Ausschlag ist. 

b) Der Dämpfer wirkt an der Lastseite. 

Für diesen Fall gelten die Differentialgzleiehungen 
I, 9, - 9)=FPsint 


F M=0 


Differentiiert man Gl. (9) nach der Zeit, so erhält man Gl. (2) und für die kritische 
Drelzahl wieder Gl. (3), das ist die der ungedämpften Schwingung. Als Lösung für y, erhält 
man ebenfalls Gl. (5). Die Werte der Konstanten €,, F,, H, bestimmen sich aber jetzt mittels 
der Gl. (9) zu 

M 
:= 
Dies sind aber die gleichen Werte wie bei der Dämpfung auf der Kraftseite, so daß man als 
Lösung des Gleichungssystems (9) wieder Gl. (a) bekommt. Für die Relativverdrehung £ 
erhält man 
P 


MI 
Pt | 


Die Konstanten der Gl. (10) gelten nur von einem Zeichenwechsel der Relativgeschwin- 
ddiekeit Ö bis zum nächsten und bestimmen sich aus den vorgegebenen Anfangsbedingungen. 
Es sei nun wieder wie früher die Annahme getroffen, daß Z bei ®t-0,=,... zu Null wird. 
Dann bestimmen sich die Relativverdrehungen bei f,—=0 

MT, ? 
und beı ot, : 
MI I, Pn 


Sollen die Ausschläge endlich bleiben, dann darf &, höchstens gleich — &, werden. Dies ist 
dann der Fall, wenn 


M ist. 
I, 4 
2. Berechnung des Dämpferträgheitsmomentes T,. 
a) Der Dämpfer wirkt auf der Kraftseite. 
Um /, zu bestimmen, sei angenommen, daß Z 0 sei und daß die Frequenz k der er- 
regenden Kraft nieht mit © zusammenfällt. Dann gelten 
(I, -@,) Psın kt, 
Daraus bestimmt sich’) bei Vernachlässigung der Eigenschwingung 
A,- Bt-+e.. 
Dabei ist die Pulsation e, 
Pic— 


Dies gilt aber wegen der Reibungskopplung nur so lange, als 


sinkt. 


ist. Nun ist 9, daher gelten auch die Ungleichungen 
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Nimmt man nun an, daß das Gleichheitszeichen gerade "noch bei At=-, gilt und die 


Pulsation 2,=a vorgeschrieben ist, bei welcher der Dämpfer zur Wirkung kommen soll, so 
bestimmt sich nach Gl. (12) das /, aus der Gleichung 


MI, 


A 


Ist /, gegenüber /, und /, vernachlässigbar, so wird 


, LM 


al, +Il,ce (13). 
b) Der Dämpfer wirkt auf der Lastseite. 
Dann gelten unter den obigen Annahmen die Differentialgleichungen 


Y,) 
Dann ist wieder unter Vernachlässigung der Eigenschwingung 
VE: B,t-+e, 
Pesinkt 
un«d 
Ks gilt wieder Gl. (11), nur ist jetzt 


Soll nun die Relatıvverschiebung von /, gegenüber /, bei kt, erfolgen und sei dann die 


Pulsation &, 9, so bestimmt sich das zugehörige /, aus der folgenden Gleichung 


M MI (I,-+ I,) 
I, +1,+1)+Pel,' 
Ist nın /, wieder gegenüber /, und /, vernachlässigbar, so erhält man die einfache Beziehung 


3. Vergleich der Ergebnisse des obigen Verfahrens mit dem von J. Ormondroyd. 


Ormondroyd untersucht nur den Fall der Dämpfung auf der Kraftseite und erhält 
für das Dämpfungsmoment 


ny2 


M= 


also einen um ea. 40°, größeren Wert als nach Gl. (S), und für das Trägheitsmoment des 


Dämpfers 
M [n\? 


ebenfalls einen größeren Wert als nach Gl. (13). 


Wird der Dämpfer nach den Angaben von J. Ormondroyd entworfen, so Ist die Ge- 
währ gegeben, daß die Ausschläge endlich bleiben. 


Zusammenfassung. 


ls wird für ein Zweimassensystem, auf das ein Reibungsschwingungsdämpfer wirkt, 
jenes Dämpfungsmoment berechnet, bei dem die Relativverdrehungen und damit die Bean- 
spruchungen endlich bleiben, wenn eine erregende Kraft Psin ot mit der Frequenz © gleich 
der der Eigenschwingung auf das System wirkt und keine Flüssigkeitsdämpfung vorhanden 
ist. Außerdem wird der Zusammenhang zwischen dem Dämpferträgheitsmoment und dem 
Winkelausschlag angegeben. Untersucht wurden die beiden möglichen Fälle (der Dämpfer 
wirkt auf der Kraftseite oder auf der Lastseite. S21 


| 
| 
| | 
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Der Einfluß von Schwingungen des Staurohrs oder von 
Geschwindigkeitsschwankungen auf die Staudruckanzeige. 


Von Georg Kiel. 


Bericht der Deutschen Versuchsanstalt für Luftfahrt. E. V.. Berlin-Adlershof. 


I. Einleitung. 


s treten zuweilen Versuchsbedingungen auf, bei denen die Staudruckanzeige entweder durch 
F Geschwindigrkeitsschwankungen oder — was auf das gleiche herauskommt durch 
Schwingungen des Staugeräts gefälscht werden kann. Bei genauen Staudruckmessungen 
außerhalb der Spannweite eines Flugzeugs, die der Verfasser im Fluge durehführte, traten 
unter ungünstigen Versuchsbedingungen am äußersten Ende der Meßeinrichtung (in einem 
Abstand von 29 m vom Flügelende) Schwingungen von erheblicher Amplitude und Frequenz 
auf. Der hierdureh entstehende Fehler in der Staudruckanzeige erhielt hierbei besondere 
Bedeutung. Da nun die Frage nach dem Einfluß von Schwingungen des Staurohrs oder von 
Geschwindigkeitsschwankungen auf die Staudruekanzeige öfter auftaucht und ganz allgemein 
interessiert, soll sie im folgenden theoretisch und versuchsmäßig näher untersucht werden 


II. Formelzeichen. 


Geschwindigkeit, 
v, Im/s] zeitlich unveränderliche Geschwindigkeit, bzw. zeitlich unveränderlicher 


Mittelwert der Geschwindigkeit, 
m/s] Geschwindigkeitsschwankung in Strömungsriehtung von 
Im/s]| Geschwindigkeitsschwankung senkrecht zur Strömungsriehtung von ®,, 
|kg/m’| Staudruck, 
Ike/m?] der aus der zeitlich unveränderlichen Geschwindigkeit », resultierende 


Staudruck, 

Ya Ike/m?] angezeigter Staudruck-Mittelwert, 

Ike/ın?| bei Schwingungen senkrecht zur Strömungsriehtung resultierender 
Staudruck, wenn der Einfluß der dabei entstehenden Schräganströmung 
vernachlässigt wird, 

Pu Ike/m’] Gesamtdruck, 

p statischer Druck, 

a Im] Amplitude der Schwingung, 

[Hz] Kreisfrequenz, 

f [Hz] Frequenz der Schwingung, 

Is] Zeit, 

Is] Schwingungsdauer, 

[Grad]  Schräganströmwinkel. 


III. Theoretische Betrachtung über die Staudruckanzeige bei Geschwindigkeits- 
schwankungen oder Schwingungen des Staurohrs. 


1. Staudruckanzeige bei Schwingungen in Strömungsrichtung. Es sei zunächst angenommen, 
daß ein Staurohr, das sieh in einer Strömung von unveränderlicher Geschwindigkeit befindet, 
Schwingungen in Strömungsriehtung ausführt. 
Relativ zum Staugerät schwankt dann die 
(Geschwindigkeit » um einen zeitlich unver- 
änderlichen Mittelwert ®,. Abb. I gibt die 
Geschwindigkeitsverhältnisse wieder. Die Ge- 
schwindigkeitsschwankung v, sei als Funktion 
der Zeit gegeben, also »,(f). Es ist dann die 
Geschwindigkeit 


AD 
A 5357 
Abb. 1. Zeitlicher Verlauf der Geschwindigkeit bei Trägheitsfreie Geschwindigkeitsmeßgeräte 
in zeigen den wahren zeitlichen Geschwindigkeits- 


verlauf an. Das gleiche gilt bei trägheitsfreien Staudruckmeßgeräten beim Staudruck. Sind 
nun aber die Trägheitskräfte nieht vernachlässigbar klein, so ist die Anzeige bei Geschwindig- 
keitsschwankungen dureh die Abhängigkeit des Meßgerätes von der Geschwindigkeit bedingt. 
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Die Anzeige von Staudruckmeßgeräten nach dem Differenzdruckverfahren (Staurohre) 
ist infolge der Beziehung 


proportional dem Geschwindigkeitsquadrat. Trägheitsbehaftete Staugeräte können also im 
eünstigsten Fall nur den zeitlichen Mittelwert des Geschwindigkeitsquadrates anzeigen. Unter 
/ugrundelegung von Gl. (1) ist 


| 
Für den vom Meßgerät angezeigten zeitlichen Mittelwert des Staudrucks ergibt sich nun: 
0 


Das Integral 7 RO dt, der zeitliche Mittelwert der Gesechwindigkeitsschwankung v,, 


über eine genügend lange Zeit T genommen, hat den Wert Null, wie auch aus Abb. I er- 
siehtlieh ist. Es ist also der angezeigte Staudruckmittelwert 


r 
Der Wert des Integrals 7 | Iv,(M)] dt stellt den Fehler der Staudruckanzeige dar. Dieser 


0 
Fehler ist stets positiv, weil der Integrand nur quadratisch vorkommt. Dies gilt jedoch nur 
unter der einschränkenden Voraussetzung, daß das Druckanzeigegerät die Mittelwertbildung 
für das Geschwindigkeitsquadrat auch richtig durehführt. Durch Leitungskapazität oder 
Drosseln kann diese Mittelwertbildung gefälscht werden. Derartige Einflüsse sollen aber bei 
der vorliegenden Untersuchung vernachlässigt werden. Es ist also unter dieser Voraus- 
setzung stets 


da > 4% bzw. > 
d. h., der bei Geschwindigkeitsschwankungen vom Staugerät angezeigte Staudruck 4. ist 
immer größer als der aus dem zeitlich unveränderlichen Geschwindigkeitsmittelwert re- 
2 


0 


sultierende Staudruck %,=", . Der hier behandelte Fall tritt auf, wenn ein Staugerät um 


seine Mittellage Schwingungen in Strömungsriehtung ausführt. 


Zu 


Die Aufgabe besteht nun darin, den Wert des Integrals n | lv,(DP dt zu ermitteln. 
0 
Den wirklichen Verhältnissen wird man am nächsten kommen, wenn für die Geschwindigkeits- 
schwankung ein sinusförmiger Verlauf angenommen wird. Unter der Annahme einer Sinus- 
schwingung läßt sich die Funktion »,(f) in der Form schreiben: 
Hierbei ist 
a Im] die Amplitude der Schwingung, 


[Hz] die Kreisfrequenz, (f[Hz] Frequenz der Schwingung), 


die Zeit. 
Unter Benutzung des Ausdrucks (4) ist dann das Fehlerglied 
T T i 
0 


| 

| 
| 
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\ 
Nimmt man an, daß bei der Staudruckanzeige der Mittelwert einer Schwingung la. h, F= 


oder ein Vielfaches davon angezeigt wird, so ergibt sich 


| 
T | ao. 
Wird nun an Stelle der Kreisfrequenz @ die Frequenz der Schwingung f= ,_ eingeführt, so 


nimmt das Fehlerglied die Form an: 


| 


Somit ist die Staudruckfehlanzeige 


| 
Unter der Voraussetzung der riehtigen Mittelwertbildung für das Geschwindigkeitsquadrat ist 
der absolute Betrag der Staudruckfehlanzeige demnach unabhängig von der Größe der je- 
weiligen mittleren Geschwindigkeit #, bzw. des daraus resultierenden mittleren Staudrucks 
| 
fo 50’, , er Ist vielmehr nur abhängig von der Amplitude und Frequenz der Schwingungen. 


Abb. 2 gibt die Staudruckfehlanzeige q, 4, bei verschiedenen Schwingungsamplituden a in 
Abhängigkeit von der Frequenz f wieder. Hierbei ist die Luftdichte o- r Ikgs’/m’'| ange- 


nommen. Um ein invariantes Schaubild zu erhalten, ist in Abb. 3 94 — 1, als Funktion von 
af dargestellt. Bezieht man den Fehler in der Staudruckanzeige auf den aus der mittleren 
Geschwindigkeit resultierenden Staudruck 4, und setzt 


“do 0 ‘fo 


kg/m‘ 


76 


Sraudruckfehlanzeige 
S 


| 
| 


— 
4 
60 60 WO 120 MOHz 0 7 2,3 ymss 
schwingungsfrequenz EEE) 
Abb. ? Staudruckfehlanzeige 44 - in Ab- Abb. 3. Staudruckfehlanzeige 
hängigkeit von der Schwingungsfrequenz / bei ga — 40 in Abhängigkeit von a-f. 


verschiedenen Amplituden «a. 


so ist klar ersichtlich, daß der Eintluß der Staudruckfehlanzeige bei gleichbleibender Amplitude 
und Frequenz um so geringfügiger wird, je größer die mittlere Geschwindigkeit v, ist. 


| 
22 221 | 
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20 20 | | 
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2. Staudruckanzeige bei Schwingungen senkrecht zur Strömungsrichtung. Es sei nun der 
Fall betrachtet, daß ein Staurohr senkrecht zur Strömungsrichtung, bzw. — was dasselbe ist — 
senkrecht zu seiner Symmetrieachse schwingt (vgl. Abb. 4). Zur wirklichen, zeitlich un- 
veränderlichen Strömungsgeschwindigkeit ®, kommt die zu dieser senkrecht gerichtete Ge- 
schwindigkeit ®, hinzu. ”, sei wiederum als Funktion der Zeit gegeben [”,,(f)]|. Abb. 4 zeigt 
schematisch die Geschwindigkeitsverhältnisse. Es ist die am Staugerät vorhandene Ge- 
schwindigkeit 


Ein trägheitsfreies und vom Anströmwinkel unabhängiges Geschwindigkeitsmeßgerät würde 
den resultierenden Gescehwindigkeitsverlauf von » nach Gl. (9) anzeigen (vgl. gestrichelte 
Kurve in Abb. 4). Das Mittel aus diesem Geschwindigkeitsverlauf entspricht dem zeitlichen 
Mittelwert der Geschwindigkeit relativ zum Staugerät (d. h. dem Mittelwert der resultierenden 
Geschwindigkeit aus », und ®,) und ist infolge der zusätzlichen Geschwindigkeit »,, stets 
orößer als die tatsächlich zu messende Geschwindigkeit ®,. Ähnliches gilt für den Staudruck- 
verlauf, der sich unter Benutzung von Gl. (9) ergibt zu: 


Ein trägheitsfreies und von der Anströmriehtung unempfindliches Staugerät würde den durch 
Gl. (10) dargestellten Staudruckverlauf anzeigen. Ist das Staugerät aber trägheitsbehaftet, 
so ist bei Vernachlässigung des Einflusses der Schräganströmung der vom Meßgerät günstigsten- 
falls angezeigte zeitliche Mittelwert des Staudrucks 


Die nun erhaltene Gl. (11) ist identisch mit der für Schwingungen in Strömungsrichtung 
erhaltenen Gl. (3). Sieht man von der Empfindlichkeit des Staugeräts von der Schräg- 
anströmung ab, so wird man, unter der Voraussetzung der gleichen Schwingungsform, dieselbe 
Fehlanzeige erhalten, gleichgültig, ob die Schwingung in oder senkrecht zur Strömungsriehtung 


1 
erfolgt. Da or, der zu messende Staudruck ist, stellt: der Wert des Integrals 


at 


(die Staudruckfehlanzeige dar. Dieser Fehler ist unter Voraussetzung der früher getroffenen 
Annahme stets positiv. 


70 
Ta | 
7 
| 
08} | 
| 
| 
0,6 
| | | 
I 
| | 
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Messung von Kumbruch 


| 
| 
G2\ nach Formel 126 berechn. 
x | 
| 
| : | 
| | 
| 
0 20° 40° 60° 
25354] Schröganstromung 0° 
Abb. 4. Die Geschwindigkeitsverhältnisse beim senk- Abb.5. Staudruekanzeige des Prandt| 


recht zur Strömungsriehtung schwingenden Staurohr. rohrs bei Schräganströmung. 


- 
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Die Schwingung des Staurohrs senkrecht zur Strömungsrichtung bringt nun aber gleich- 
zeitig auch eine sich periodisch ändernde Schräganströmung mit sich. Staurohre und Stau- 
düsen zeigen bei größerer Schräganströmung einen zu kleinen Wert an. Soll die Fehlanzeige 
bei Schwingungen des Staurohrs allgemeingültig erfaßt werden, so ist Gl. (11) noch in der 
Weise zu erweitern, daß auch der Einfluß der Schräganströmung des Staugeräts miterfaßt 
wird. Im folgenden soll der Einfluß der Schräganströmung für das normale Prandtlrohr 
berücksichtigt werden. Abb. 5 zeigt die Abhängigkeit der Staudruckanzeige des normalen 
Prandtlrohrs von der Schräganströmung nach Messungen von Kumbruch'). 


Um den Einfluß der Schräganströmung formelmäßig in Gl. (11) einführen zu können, 
soll die von Kumbruch gemessene Kurve durch eine Fourier-Reihe angenähert werden. 
Die abzuleitende Formel braucht die Verhältnisse nur über einen verhältnismäßig kleinen 
Schräganströmbereieh zu erfassen. Bezeichnet g° den Staudruckwert, der unter dem Einfluß 
der Schwingungen, aber bei Vernachlässigung des Einflusses der Schräganströmung erhalten 
wird, 4. den bei Berücksichtigung des Schräganströmeinflusses angezeigten Staudruckmittelwert 
und o den Anströmwinkel, so kann man setzen: 


Ya 


r r ( ® r Ya 
Werden der gemessenen Kurve =f(0) je 3 einander zugehörige Werte für = und o ent- 


nommen, so ergeben sich nach Auflösung der 3 Gleichungen für die einzelnen Konstanten 
folgende Werte: 

Durch Einsetzen dieser Werte in Gl. (12a) und Berücksichtigung der trigonometrischen 
Beziehungen 


eosBo coso(deos’o 3) und cos (16 cos! — 20 cos’ 6 +9) 


ergibt sich 
cos o cos’ a 4cos’ o 
In Abb. 5 ist die nach dieser Formel berechnete Kurve gestrichelt eingezeichnet. Die 
Übereinstimmung der gerechneten mit den von Kumbruch gemessenen Werten ist bis zu 
30° Schräganströmung recht vollkommen. Größere Schräganströmungen interessieren praktisch 
nicht mehr. 


= s. Unter der Voraussetzung, daß in jedem Augenblick, 
+ 

Ja 


wo der Anströmwinkel o einen bestimmten Betrag o, erreicht, das Verhältnis a denselben 


Nun ist cos © 


Wert hat, den es bei über längere Zeit konstant gehaltenem o, besitzt, wird hiermit unter 
Berücksichtigung von (10) der angezeigte Staudruckmittelwert 


und 

2,41 +7,44(°°) 1601 


Hierbei ıst 


Wird für »,(M) wieder, wie weiter oben, eine Sinusschwingung angenommen und in Analogie 


zu Gl. (6) 


eesetzt, so wird 


Forschungsarbeiten des VDI, Heft 240, 


MH. Kumbruch: „Messung strömender Luft mittels Stangeräten.“ 
19?1 (Abb. 18). 
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Für den Fall, daß z =0,2?), bzw. — 0,045 ist, verschwindet, wie aus Abb. 5 zu ersehen 


", 
ist, der Einfluß der Schräganströmung. Es gilt dann in gleicher Weise wie bei der Schwingung 
in Strömungsrichtung die einfache Formel 


fa 


f? 


Nach den bekannten Regeln der analytischen Geometrie läßt sich auf Grund von Gl. (14) 


leicht der erreiehbare Höchstwert von Ja Yerechnen. Durch Differentiation von Gl. (14) nach 


. und Nullsetzen des für Tat) gefundenen Ausdrucks ergibt sich: 


Die Auflösung ergibt die zu (2a) gehörigen Werte: 


do 
0,8572 bzw. + 


Durch Einsetzen dieser Werte in Gl. (14) ergibt sich 


1.08. 


Der bei Schwingungen des Prandtlrohrs senkrecht zur Strömungsrichtung erreichbare 
Fehler kann also bei Zugrundelegung von Gl. (14) bei allen Amplituden und Frequenzen im 
ungünstigsten Fall höchstens S"/o q, betragen. Gleichzeitig läßt, sich die Bedingung angeben, 


die beim Erreichen des Wertes (°) erfüllt sein muß. Nach dem Vorauszehenden ist der 


Io IN 


zu | ie) gehörige Wert von 


‘fo IMAX 


0.8572 


af 
Durch Auflösen dieses Ausdrucks nach / ergibt sich 


0,02, 


af 


Man erreieht also den Größtfehler 


af 
‘Yo 

Ist. 


Zur Veranschaulichung der Größe der zu erwartenden Staudruckfehlanzeige sind in 
Abb. 6 die bei verschiedenen Geschwindigkeiten (», 20; 30; 40: [m/s]) für eine Amplitude 


von «003m nach Gl. (14) berechneten Fehler 100 am { in Abhängigkeit von der Frequenz 
0 

der Schwingung aufgetragen. Zum Vergleich sind gleichzeitig auch die einer Schwingung in 

Strömungsrichtung entsprechenden Kurven nach Gl. (S) gestrichelt eingezeichnet. Die Fehler- 

kurven weisen, was bereits nach Formel 15 vorausgesagt wurde, bei allen Gesehwindigkeiten 

", als größte Staudruckfehlanzeige q, auf. 


11,59, 


3») Die übrigen Wurzelwerte haben keinen Sinn, weil sie außerhalb des Geltungsbereiehes der Gl. (1% liegen. 


(12a). 
| 
— 
| 
| 
| 
} 
| 
| 
” 


296 Kiel,.Der Einfluß von Schwingungen des Staurohrs auf die Staudruckanzeige Math. und Mech. 


Kin allgemein gültiges’Schaubild der Staudruckfehlanzeige, das gegen eine Anderung 
der Einzelgrößen invariant ist, wird erhalten, wenn die Staudruckfehlanzeige in Abhängigkeit 


von dargestellt wird. Abb. 7 zeigt dieses Schaubild, in dem der Verlauf der Staudruck- 


20 717717 
/ 
I / 
sol #0, %% | | 
c 74 + / / 4 S / 
> / / / S 1 emessene Werte: 4 
/ / / I / Schwingung senkrecht 
S / / I H zur Strömungsrichtung 
> NW 
Schwingung 
R S | / in Stromungsrichtung 
S 
S S6 
4 
2 
2 60 68 120 760 
Frequenz f 002 00% 0066 008 00 0% 
20 A567] af 
Abb. 6.  Standruekfehlanzeige in Ab- Abb. 7. Standruckfehlanzeige 
/o 1) 


hängigkeit von der Freqnenzf bei einer Amplitude 

von m und verschiedenen Geschwindig 
keiten vn. 

Schwingung senkrecht zur Strömungsriehtung. 

Schwingung in Strömungsriehtung. 


als Funktion von 


fehlanzeige [100 als Funktion von sowohl für die Schwingung senkrecht zur 
0 0 

Strömungsriehtung als auch für die in Strömungsriehtung dargestellt ist. Sind die jederzeit 

leicht meßbaren Größen Amplitude «af[m]|, Frequenz f|Hz] und Geschwindigkeit [m/s] 


a 
bekannt, so läßt sich ohne langwierige Rechnung dureh Bildung des Ausdrucks ne sofort die 


0 
zu erwartende Staudruckfehlanzeige aus Abb. 7 ablesen. In Wirklichkeit werden wohl selten 
reine Schwingungen in Strömungsrichtung oder senkrecht dazu auftreten. Die beiden Kurven 
der Abb. 7 würden demnach Grenzkurven darstellen, zwischen denen die praktisch auf- 
tretenden Fehler zumeist liegen dürften. 


IV. Ergebnis der Staudruckmessung mit einem senkrecht zur Strömungsrichtung 
schwingenden Prandtlrohr und Vergleich der gemessenen Staudruckfehlanzeige mit 
den berechneten Werten. 


Die vorangegangenen theoretischen Betrachtungen setzten Sinusschwingungen des Stau- 
rohrs voraus und gingen weiterhin von der Annahme aus, daß bei der Staudruckanzeige in- 
folge der 'Trägheit und Dämpfung der Mittelwert der Staudruckschwingung, also ein dem 
Mittelwert des Geschwindigkeitsquadrates proportionaler Betrag angezeigt wird. Um nun ein 
Bild über die Größe der tatsächlich praktisch auftretenden Fehlanzeige des Staudrucks zu 
bekommen, wurde diese für den Fall der Schwingung senkrecht zur Strömungsrichtung ver- 
suchsmäßig im Instrumenten-Windkanal der Deutschen Versuchsanstalt für Luftfahrt (DVL) 
bestimmt. 


1. Versuchsanordnung und Versuchsdurchführung. In Abb. S ıst die Versuchsanordnung 
schematisch dargestellt. Abb. 9 zeigt die in den Windkanal eingebaute Versuchseinriehtung 
im Lichtbild. Die zur Strömungsriehtung senkrechten Schwingungen des Prandtlrohrs wurden 
dureh einen von einem Gleichstrom-Nebenschlußmotor angetriebenen Kurbelbetrieb erzeugt. 
Durch Einschalten entsprechender Widerstände in dem Nebenschluß konnte die Drehzahl und 
hiermit die Frequenz der Schwingungen weitgehend verändert werden. Als schwingendes 
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Staurohr wurde ein normales Prandtlrohr von 6 mm Durchmesser verwendet. Um einer- 
seits ein Abbrechen des Prandtlrohrs infolge des bei den sehnellen Schwingungen auf- 
tretenden Einspannmomentes zu verhüten und andererseits zusätzliche Eigenschwingungen des 
Prandtlrohrs so gut wie möglich zu vermeiden, wurde das Prandtlrohr nach der Ein- 
spannstelle zu konisch verstärkt. 


feststehendes 
Prandtirohr 
schwingendes 
Prandtirohr 
> 
’_a-30mm 


Abb. 9. Die in den Windkanal eingebaute 
Versuchseinriehtung. 


Abb. 8. Schema der Versuchsanordnung. 


Das Meßverfahren bestand darin, daß bei verschiedenen Drehzahlen des Motors, also 
verschiedenen Schwingungsfrequenzen des Prandtlrohrs, die Staudruckanzeige des schwingenden 
Prandtlrohrs q,„ und gleichzeitig die eines feststehenden q, bestimmt wurde. Der Staudruck 
des Windkanals q, wurde während jeder Meßreihe, so gut es möglich war, konstant gehalten. 
Zur Staudruckmessung wurden Alkohol-Mikromanometer benutzt. Vor Beginn der Versuche 
wurde die Staudruckverteilung in dem Schwingungsbereich des schwingenden Prandtlrohrs 
bestimmt. Sie erwies sich als konstant. 


2. Ergebnis der Messungen und Vergleich mit der Rechnung. Das Ergebnis der Messungen 
der Staudruckanzeige des schwingenden Prandtlrohrs ist in Abb. 10 wiedergegeben. Die 
Staudruckfehlanzeige ist in Prozenten des 2 


wirklichen Staudrucks als Funktion der Motor- 90 = hgym 


8 + 278 m/S 
drehzahl bzw. der Frequenz der Schwingung | De ı 
aufgetragen je — — f(n) bzw. — | 


Es wurden zwei Versuche durchgeführt, der DT | 
eine bei einem mittleren Staudruck von 
455 [kg/m?] und 27,8 [m/s], der an- 
dere bei 9,83 [kg/m?] und 37,6 [m/s]. 
Die Amplitude der Schwingung betrug 30 mm. 
Neben den Meßergebnissen sind gleichzeitig 
in Abb. 10 noch die nach den weiter oben ab- 
geleiteten Formeln berechneten Fehlerkurven 
eingezeichnet. Mit zunehmender Schwingungs- 
frequenz werden die gemessenen Staudruck- 
anzeigefehler größer als die berechneten. 0 
Es ist zu bemerken, daß die Meßwerte L | | | | 
wegen der Kleinheit der auftretenden Stau- f 
druckdifferenzen etwas unsicher sind. Die 


So 


| 9583 kg/m? wi 

40.376 

| =30mm \ 
4 


500 7000 7500 2000 T/mın 
Drehzahl 


Sraudruckfehlanzeige 700 


berechn. \ 


| 
35 Hz 


Jeweils zugehörige gleichzeitige Ablesung der 


Abb. 10. Gemessene Standruekfehlanzeige in Abhängigkeit 
vonder Motordrehzahl» bzw. der Frequenz der Schwingung, 
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beiden Staudrücke gestaltete sich etwas schwierig, weil die Geschwindigkeit des Instrumenten- 
Windkanals zeitweise ein wenig schwankte, was gerade bei den vorliegenden Messungen als 
störend empfunden wurde. In einem Windkanal mit besserer Geschwindigkeitsregelung ließen 
sich wahrscheinlich genauere Messungen erzielen. Bei der vorliegenden Aufgabe handelt es 
sich jedoch nur darum, festzustellen, ob bei derartigen Schwingungen überhaupt meßbare 
Fehler auftreten. Für die größenordnungsmäßige Bestimmung der Fehler dürfte die Meß- 
genauigkeit ausreichen. 

Der Grund dafür, daß die berechneten Fehlerwerte kleiner sind als die gemessenen, 
kann einmal in den der Rechnung zugrunde liegenden verhältnismäßig einfachen Annahmen 
liegen. Der Rechnung wurde die Annahme zugrunde gelegt, daß ein dem Geschwindigkeits- 
uadrat proportionaler Mittelwert angezeigt wird. Dieser Fall ist der günstigste, und die 
unter dieser Voraussetzung berechneten Fehler sind die Mindestfehler, die unbedingt auftreten 
müssen. Daß größere Anzeigefehler auftreten, als die Rechnung erwarten läßt, kann auch in 
den Eigenschaften der Versuchseinrichtung, der Anzeige-Instrumente und der Druckleitungen 
seinen Grund haben; denn wie bereits weiter oben erwähnt, wurden bei den theoretischen 
Betrachtungen diejenigen Anzeigefehler nieht berücksichtigt, die durch Leitungskapazität oder 
Drosseln verursacht werden können, Trotz der gegenüber den schwingenden Massen ver: 
hältnismäßig großen Masse der feststehenden Teile der Versuchseinriehtung und des Motors 
und trotz teilweisen Massenausgleichs der hin- und herschwingenden Teile ließen sich Eigen- 
sechwingungen der gesamten Versuchseinriehtung nicht ganz vermeiden. Es muß weiterhin 
dlamit gerechnet werden, daß auch der Kopf des schwingenden Prandtlrohrs selbst zusätz- 
liche Schwingungen ausgeführt hat, so daß die wirkliche Amplitude wahrscheinlich etwas 
größer war als der bei der Rechnung eingesetzte Kurbelradius. Die Abweichung der beim 
Versuch gewählten Schwingungsform von der rechnerisch zugrunde gelegten reinen Sinus- 
schwingung ist praktisch belanglos?). Wenn auch die gemessenen Fehler offenbar größer sind 
als die berechneten, so können die abgeleiteten Formeln trotzdem zur Abschätzung des 
mindestens auftretenden kleinsten Fehlers dienen, so daß in allen Fällen, in denen sich schon 
rechnerisch unzulässig hohe Fehler ergeben, die wirklichen bestimmt nicht mehr vernach- 
lässigt werden dürfen. 


V. Zusammenfassung und Schlußfolgerungen aus dem Ergebnis der Untersuchung. 


Bei Schwingungen des Staurohrs bzw. bei Geschwindigkeitsschwankungen kann die 
Staudruckanzeige gefälscht werden. Die Fehlanzeige ist in der Abhängigkeit des Staugerätes 
von der Geschwindigkeit begründet. Vernachlässigt man die Anzeigefehler, die durch Leitungs- 
kapazität oder Drosseln verursacht werden können, so ist die Anzeige von Staudruckmeß- 
geräten nach dem Differenzdruckverfahren dem Geschwindigkeitsquadrat proportional. Träg- 
heitsbehaftete Staugeräte können deshalb im günstigsten Fall nur einen dem zeitlichen 
Mittelwert des Geschwindigkeitsquadrates proportionalen Wert auzeigen. 


Es wurden zwei Fälle untersucht, und zwar: 
1. Schwingungen des Staurohrs bzw. Geschwindigkeitsschwankungen in Strömungsrichtung, 
2. Schwingungen des Staurohrs senkrecht zur Strömungsrichtung. 


4 Bei den abgeleiteten Formeln wurde für ey=aosinot angenommen. Die genaue Schwingungsform des 
Prandtlrohrs bei den Versuchen würde ergeben: 


Hierbei ist 2 ji — 5 0,231 das Schubstangenverhältnis. Für das Fehlerglied ergibt sieh dann in analoger Weise wie 
weiter oben 
4 / 1 


- 


5) 
Wird nun über eine Sehwingungsdauer T = n integriert, so erhält man: 


T 
| 4] 


Die Abweichung ist also vernachlässigbar klein. 
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Unter der Voraussetzung, daß das Staurohr den Mittelwert einer oder mehrerer Schwingungen 
liefert, erhält man bei der Annahme einer Sinusschwingung von der Frequenz f und mit der 
Amplitude a bei Schwingung in Strömungsrichtung für die Staudruckfehlanzeige 


Der absolute Betrag der Staudruckfehlanzeige ist unabhängig von der Größe der jeweiligen 
zeitlich unveränderlichen mittleren Geschwindigkeit ®, bzw. des daraus resultierenden Stau- 

0 

drucks q, = 9 


Bei Schwingungen senkrecht zur Strömungsriehtung tritt gleichzeitig eine sich periodisch 
ändernde Schräganströmung auf. Die durch die Schwingung hervorgerufene stets positive 
Staudruckfehlanzeige erfährt dureh die Schräganströmung gleichzeitig wieder eine Verminderung. 
Es wurde nun angenommen, daß beim Schwingen des Prandtlrohrs in jedem Augenblick, 
wo der Anströmwinkel einen bestimmten Wert erreicht, der Einfluß der Schräganströmung 
derselbe ist wie bei gleicher stationärer Schräganströmung. Unter den gleichen Annahmen 
wie bei der Schwingung in Strömungsriehtung (Sinus-Schwingung und Integration über eine 
Schwingungsdauer) und Berücksichtigung des Einflusses der Schräganströmung ergibt sich 
dann beim senkrecht zur Strömungsriehtung schwingenden Pr andtlrohr für das Verhältnis 


or 
von angezeigtem Staudruck q. zum tatsächlich zu messenden 4, vo 
2 
za = — | | 
‘fo 
Hierbei ist 


Der bei Schwingungen des Prandtlrohrs senkrecht zur Strömungsriehtung nach Gl. (14) 
erreichbare Größtfehler kann bei allen Amplituden und Frequenzen im ungünstigsten Fall 
höchstens S°/, q, betragen. Man erreicht den Größtfehler 


(de) — 1,08, wenn 
‘fo 0 

z dargestellt, so ergibt sich ein Schau- 

bild, das gegen eine Änderung der Einzelgrößen invariant ist (vgl. Abb. 7). Sind Amplitude a [m], 
Frequenz f[Hz] und Geschwindigkeit "„ [m/s] bekannt, so läßt sich ohne langwierige Reehnung 


Wird die Staudruckfehlanzeige in Abhängigkeit von 


durch Bildung des Ausdrucks sofort die zu erwartende Staudruckfehlanzeige ablesen. Die 


0 
praktisch auftretenden Fälle dürften in der Regel zwischen den beiden Grenzfällen, Schwingung 
in Strömungsrichtung und Schwingung senkrecht zur Strömungsrichtung, liegen. 


Für den Fall der Schwingung senkrecht zur Strömungsrichtung wurde die Staudruck- 
fehlanzeige durch Windkanalversuche für eine Amplitude von a = 30 mm bestimmt (vgl. Abb. 10). 
Die gemessene Staudruckfehlanzeige stimmte nur bei kleinen Frequenzen mit den berechneten 
Werten überein. Mit zunehmender Frequenz wurden die gemessenen Werte größer als die 
berechneten. Die Abweichung der Meßergebnisse von den berechneten Anzeigefehlern dürfte 
teilweise in den verhältnismäßig einfachen Rechnungsannahmen ihren Grund haben, teils 
auch in der ungenauen Erfassung der Schwingungsamplitude und in der Vernachlässigung der 
Anzeigefehler, die durch Leitungskapazität und Drosselung entstehen können. Da das Ziel 
der Aufgabe nur die Ermittlung der Größenordnung der auftretenden Staudruckanzeigefehler 
war, wurde die Ursache der Abweichung der Meßwerte von den berechneten Werten nicht 
weiter verfolgt. 

Zur Abschätzung der Größenordnung der auftretenden Staudruckanzeigefehler dürfte die 
Genauigkeit der abgeleiteten Formeln ausreichen. Im allgemeinen sind die praktisch auf- 
tre tenden Schwingungen der Staugeräte für die Staudruckanzeige belanglos, weil einerseits 
die hierbei möglichen Fehler zumeist bedeutungslos sind und andererseits unliebsame größere 
Schwingungen leicht vermieden werden können. In Sonderfällen, z. B. bei genauen Staudruck- 
messungen am fliegenden Flugzeug und auch unter besonderen Versuchsbedingungen im 
Laboratorium, sind Schwingungen der Staurohre jedoch zuweilen zu beachten, wenn ihre 
Amplitude und Frequenz größere Ausmaße annehmen. Das gleiche gilt natürlich auch von 
Geschwindigkeitsschwankungen, sofern der Staudruck der zeitlich unveränderlichen mittleren 
Geschwindigkeit v, gemessen werden soll. 535 
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Über asymptotische Reihenentwicklungen in der 
Operatorenrechnung. 


Von Ake Pleijel in Stockholm. 


Zuerst will ich in Kürze einige Formeln aus der Heavisideschen Operatorenreehnung 
angeben. 
In der Operatorengleichung 


ist das Zeichen für y(p) ist die bekannte Operatorenfunktion und wird gesucht. 


Die Gl. (1) wird hier am besten als äquivalent mit der komplexen Integraldarstellung von 


hit) definiert: 
Il (ent 


R—n-- . 


wo der Integräationsweg die rechte Hälfte eines Kreises um den Ursprung mit dem ins Un- 
endliche wachsenden Radius A ist. Unter den unten angegebenen Voraussetzungen, denen 
y(p) in physikalischen Problemen gewöhnlich unterworfen ist, wird für positive f-Werte das 
entsprechende Integral über dem linken Halbkreis null, wenn R ins Unendliche übergeht. 
Da wir in dem folgenden hit) für #>0 zu betrachten haben, können wir die Integration über 
den ganzen Kreis erstrecken. 

Über Operatorengleichungen der Form == fi) braucht hier nur gesagt zu werden, 
daß die in Frage kommende Operation unter ziemlich allgemeinen Voraussetzungen definiert 
werden kann, und daß diese Definition so wird, daß die Gleichung 

gültig ist. Weiterhin gilt die sog. Verschiebungsregel 
fiti: eat.g (p+a)-fih 
und für die formelle Behandlung von (1) ist die Operatorenbeziehung 


von Wichtigkeit. Für gılt 


1 Cet 
\ 


(3a), 


wo am ganzen Kreis entlang integriert wird'). 
Im folgenden soll g(p) so beschaffen sein, daß der Grenzwert (2) existiert. Außerdem 
soll gelten, daß 


für 


Endlich wird angenommen, daß g(p) eine, in der ganzen Ebene mit Ausnahme einer end- 
lichen Anzahl Punkte, reguläre analytische Funktion ist. Wenn a eine, im Endlichen be- 
legene, singuläre Stelle ist, soll eine für genügend kleine Werte von |p - a, konvergente 
Iintwieklung 


gelten, wo % rational >VO und a willkürlich ist. 

Die Absicht der vorliegenden Arbeit ist, mit Hilfe von (2) und unter obenstehenden 
Voraussetzungen betreffend g(p) eine formelle Methode zu studieren, die die Lösung von (1) 
in Form einer asymptotischen Reihe ergibt. In den meisten Problemen der angewandten 
Operatorenreehnung, wo es sich um asymptotische Entwicklung der Lösung handelt, ist g den 
obenstehenden Voraussetzungen unterworfen. 


!) Liieratur: H. Jeffreys, Operational Methods in Mathematical Physies. Cambridge 1931. 
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Die Methode, die auf Heaviside zurückgeht, geht aus folgendem hervor: Wegen der 
Verschiebungsregel gilt (a ist eine singuläre Stelle von ' = 


1 eat eat - eat (pra)e al 
Außerdem ist 
p 
at — 


Also ist die Operatorengleichung (1) äquivalent mit 


e 


hit=et!. 


p 
pra 
Nach den Annahmen, die wir zugrunde gelegt haben, gilt für genügend kleine |p| eine 
konvergente Entwicklung 


Wird diese Reihe nun in die Operatorengleichung (4) eingesetzt, und läßt man die Reihe 
Glied für Glied operieren, so erhält man wegen (3) eine Reihenentwicklung MH, nach ab- 
nehmenden Potenzen von #, mit e"* multipliziert. Die formelle, im allgemeinen divergente 
Reihe H„ wird im folgenden die zu a hörende Heavisidesche Entwicklung genannt. Es 
ist leicht zu verstehen, daß weder das gliedweise Operieren noch die Anwendung von () 
zulässig ist. Trotzdem zeigt sich manchmal, daß die erhaltene Reihe e*!.H, eine verwend- 
bare, asymptotische Entwicklung von h(f) für £>0 ergibt. Man erhält jedoch eine solehe 
Heavisidesche Reihe für jede Singularıtät, und man kann sich fragen, ob man mit Hilfe 
dieser Reihen eine asymptotische Entwicklung für A(t) herleiten kann, die für # reell, > 
gültig ist. 

Eine verwandte Frage ist früher von J. R. Carson?) und N. Wiener’) behandelt 
worden. Diese Autoren gehen von der Annahme aus, daß der Ursprung schon anfangs eine 
sineuläre Stelle ist und fragen dann, unter welchen Umständen die daraus sich ergebende 
asymptotische Reihe A(f) darstellt. Die untenstehende Methode ergibt, unter den gemachten 
Voraussetzungen, immer das asymptotische Verhalten von hf) für Freell>0. Man erhält 
also auch den stationären Endzustand. Unsere Resultate werden mit Hilfe einer gewöhn- 
lichen Abschätzung der Integrale des Typs (2) gewonnen, und die Anwendung dieser Ab- 
schätzung ist wegen der nahen Verwandtschaft der Operatorenreehnung mit der Theorie der 
gewöhnlichen Differentialgleichungen natürlich. 

Wir gehen aus von 


RR 


Bevor man integriert, muß man jedoch p(p) dureh Sehnitte in der p-Ebene eindeutig machen 
und einen Funktionszweig auswählen. Das letztere ist in diesem Zusammenhang von 
geringerer Bedeutung. Die Schnitte werden auf folgende Weise gelegt: Der Punkt a soll ein 
Verzweigungspunkt der Funktion g(p) sein. Er wird von einem Kreis e umgeben 


p-al=r, 
in dem r so klein genommen ist, daß die übrigen Singularitäten völlig ausgeschlossen werden 
(so daß sie noch nieht einmal auf der Peripherie liegen). Der Schnitt soll nun in e von a 
zu (a r) einem Radius entlang gehen. Außerhalb von e soll der Schnitt, ohne andere 


Singularitäten von zu treffen, (a — r) mit reellen — verbinden. Dieser letzte Teil, S, 


des Schnittweges wird als eine kontinuierliche Kurve mit kontinuierlich variierender Tangente 
angenommen. Die Richtungskoeffizienten der Tangente sollen gleichmäßig beschränkt sein, 
so daß, wenn p=.x#.-+iy, S entlang eine Abschätzung mit konstantem k, 


dp 


gültig ist. 
Für den Teil S der Schnittkurve gilt der Hilfssatz: Wenn die Funktion Y(p) entlang 5 
kontinuierlich ist und einer Ungleichung 


2), J. R. Carson: Eleetrie Cirenit Theory And The Operational Caleulus, 1926. 
3) N. Wiener: The Operational Caleulus, Math. Ann. 95, 1926. 
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genügt, wo K und x reelle, positive Konstanten sind und wie oben p=x=-+iy, dann ist 


ta 
\ 


unabhängig von gq. 


Die Gültigkeit dieses Satzes geht daraus hervor, daß für t>x 


a—r 
S — 
Das Integral (6) kann nach Cauchys Integralsatz in eine Summe von Integralen auf- 


veteilt werden, die entlang gewissen, die Singularitäten von umschließenden, aber sonst 


willkürlichen Kurven genommen werden. Diese Kurven werden auf folgende Art und Weise 
gelegt: Für einen Pol: ein Kreis um denselben, klein genug, um die übrigen Singularitäten 


von vollständig auszuschließen. Für eine Verzweigungsstelle: auf der einen Seite $, 


von dem dazugehörenden äußeren Schnitt von — x zu (@a—r). Dann um e und wieder zurück 
zu — oo an der anderen Seite des Schnittes, S,, entlang. Alles in positiver Richtung. 


Wir studieren den Beitrag I,„, den man von der Verzweigungsstelle a erhält. Es gilt, daß 


l 
Sı 


Setzt man hier das erste Integral auf der rechten Seite gleich A und setzt man, ent- 
sprechend (5), in das andere 
(p+ a) 4 
pra p 
nachdem man erst durch eine Variabelntransformation p dureh (p-+a) ersetzt hat (wobei 


ec, S,,S, in e’,S’,,S’, verwandelt werden), so erhält man 


"opti 


p 2 


( 


Der Interrationsweg im mittleren Integral kann nach Caucehys Satz durch die trans- 
formierten Wege 5’,, S’, (in entgegengesetzter Riehtung) und einen unendlich großen Kreis 
um den Ursprung ersetzt werden. Der Kreis wird positiv durchlaufen. Das Integral um 
diesen unendlich großen Kreis ergibt, wie man leicht nach (3a) ersieht, gerade den Anfang 
der Heavisideschen Entwicklung MH, um a. Die Integrale längs S, und S, werden zu- 
sammen mit B bezeichnet. Es ergibt sich, daß 


eat ent 


(4-+ Be@at)—+ (Anfang der Heavisideschen Entw.) 


c’ 


In dieser Gleichung kann der Hilfsatz auf A und B angewandt werden, so daß unab- 
hängig von 


lim (A+Be—0, 


Um zu zeigen, daß der Anfang der Heavisideschen Entwicklung eine semikonvergente 
Reihe für, ergibt, braucht man, da die niedrigste Potenz der Entwicklung t-W""+«) ist, 
e 


nur zu zeigen, daß 


. l dar > 
Um Pure, 
2ri)p 
c’ 
Um dies zu beweisen, nehmen wir an, daß wir n so groß gewählt haben, daß 
Das Integral 
I (et 
ar > 
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n ist das längs des in der untenstehenden Abbildung bezeichneten Weges genommen ist, wird nach 
| Cauehy null. Wegen (7) nähert sich das Integral um den inneren Kreis dem Grenzwert 
null, wenn der Radius des Kreises gegen null strebt. Es ergibt sich 


1 
ar p 
ec’ 
1 ep! 
t—x' -r 
len auf- 1 Cent 
2ri)p 
ber sonst 
d Weise (das erste Integral den oberen Weg, das zweite den unteren entlang). 
laritäten Wenn der Maximumwert von 
- zurück auf dem Kreis |Z|=o>r* (o< der Konvergenzradius von ($8)) mit M bezeichnet wird, so gilt 
nach Cauchys Koeffizientenschätzung 
M 
daß 
und daraus erhält man für |p|=r die Schätzung des Restglieds 
n-+ | 
ın, ent- n | + l + n + 1 + —+ pt rl | 
Für das Restgliedintegral ergibt sich wegen der obengeschriebenen Darstellung dieses 
Integrals und wegen der Restgliedschätzung: 
r 
(wobei ent 
» 
(kin «) 
< konst | e ch in 4 de == konst 
trans- 
ie Damit ist bewiesen, daß die Heavisidesche Entwicklung für den Punkt a eine semi- 
um 
Anfang konvergente Reihe für = ergibt. Dasselbe gilt für die Pole, wobei man aber nur endlich 
en zu: viele Glieder #* erhält, wo s ganz und 0 ist. Jede Heaviside- Entwicklung ist asymptotisch 
in dem Sinne, daß 
— Teils — h 


wenn f=-"die höchste, vernachlässigte Potenz ist. 
unab- Wir wollen nun die Summe der Beiträge /,„ für große, reelle und positive #-Werte 
| studieren. Nach Cauchys Satz gilt 


| 
a a 
rgente | wo I. eine Funktion ist, die die semikonvergente Entwicklung H„ besitzt. Wir betrachten 
ist, die Singularitäten @ von die den größten Realteil haben. Die Beiträge von diesen 


werden, wegen ihrer Exponentialfaktoren, für große t-Werte ausschlaggebend. Es ist leicht 
zu zeigen, daß der Fehler, den man macht, wenn man 

Set. 

d a a 

schreibt, wo H’a Teilsummen der entsprechenden, semikonvergenten Reihen sind, von größerer 
Ordnung ist als die Summe der übrigen Beiträge, unabhängig von der Wahl der Teilsummen. 
Dies gilt, wenn unter den Reihen wirkliche Reihen mit unendlich vielen Gliedern sind. Das 
Mitnehmen der übrigen Beiträge ist in diesem Falle ganz ohne Nutzen. Man erhält die 
folgende Regel: 


| 


Ztsehr. f. angew. 
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Die Heavisideschen Entwicklungen werden für die Singularitäten @ von r gebildet, 
) 
die den größten Realteil « haben. Die Summe 
a 


ergibt eine asyımptotische Entwicklung der Funktion A (f) für reelle, positive #-Werte. 
Die Art dieser Entwieklung ist die folgende: Man hat 


hit) 
Hate), 

a 


wo wie oben 41’5 Teilsummen sind. Die höchsten vernachlässigten Potenzen sollen t=-#ı, 
Ih»... sein. Dann gilt 
e()=O(t-R), 
Wo 
Dies ergibt sich ohne Schwierigkeit, wenn man beachtet, daß die Faktoren e“—-«t be- 
schränkt sind. 

In &inigen Fällen ist es indessen richtig und vorteilhaft, auch andere Entwicklungen 
als die von @ herrührenden mitzunehmen. Dies gilt, wenn die Entwicklungen der Beiträge 
von «a nur endlieh viele Glieder enthalten (wie für Pole). Denn dann ist der oben besprochene 
Fehler null. Wie die Regel in diesem Falle zu ändern ist, ist leicht einzusehen. 


Als ein Beispiel möchte ich folgendes vorbringen: Nach J. R. Carson kann man die 

„Heavisidesche Regel” auf der Operatorengleichung 

nieht anwenden. Es handelt sieh bei ihm nur um die Entwicklung um den Ursprung. Nach 
der oben angeführten Regel erhält man eine asymptotische Entwicklung für /(t), wenn man 
die von +7 herrührenden Beiträge zu der von J. R. Carson angegebenen Reihe hinzufügt. 
Diese Beiträge ergeben zusammen: 


In=| 


R 


Die von J. R. Carson angegebene Entwicklung (um den Ursprung) ist 


. 
| 


( (1 1:3-5 
| aRt Pot 
| 1.3.5 \ 
akt Bot 


Man hat also 


was auch den stationären Endzustand ergibt. 
Uber das obenstehende Beispiel und über die Bedeutung unserer Methode weise ich auf 
das zitierte Buch von J. R. Carson, Seite 67, hin. 


ee . 


/um Sehluß möchte ich eine Sache von mathematischem Interesse anführen. Die 
Voraussetzung 


lım 0 


p 


die für = Zare p galt, soll nun für alle Werte des Arguments von p gelten. Dann kann 
man die obenstehende Methode auch für <O und sogar für komplexe f anwenden. Wenn 
wir eine asymptotische Reihe, gültig in Richtung vom Ursprung nach dem Punkt (cos y 
sin suchen wollen, werden die Singularitäten p = --iy von Bedeutung, die den größten 
Wert von 
(„er cos + SIn 


ergeben. Man kann leieht mit Hilfe der konvexen Hülle der Singularitäten von ep Rich- 


tungen angeben, in denen verschiedene asymptotische Entwicklungen gültig werden. Dies 
tritt z.B. bei der Besselschen Funktion ‚/,(f) ein, die durch die Operatorengleichung 


rereben Ist. 233 


f. angew. 


i 15, Heft 5 
ınd Mech. Band 15, Heft 5 


Oktober 1935 Buchbesprechungen 305 


BUCHBESPRECHUNGEN 


ebildet, 
H. SEIFERT u. W. THRELFALL, Lehrbuch Komplexes und der Theorie der Überlagerung ge- 
ter Topologie. \ll u. 353 8. mit 132 Fig. widmet. Im neunten Kapitel werden die dreidimen 
Leipzig u. Berlin 1934, B. G. Teubner. Geb, 20 M. sionalen Räume, im zehnten die mehrdimensionalen 
Räume behandelt. Hier sieh an auf Sehritt uni 
Die um die Mitte des vorigen Jahrhunderts vom wie viele der 
uortinger Mediziner Listin S begründete und kurz Fragen noch immer unbeantwortet sind. Dann folet 
darauf von B. Riemann mit größtem Erfolg auf 5 
5 die Theorie des Abbildungsgrads. die man Brou- 
je Funktionentheorie angewandte Topologie ist 
lanere Zeit als ei sehr spezieller Zweie ler Geo wer verdankt, sowie einILe Fixpunktsätze. Das 
Mathematiker von Interesse zu sein schien, 
‚re Dalze ‚weilel Se St ? 
der Mechanik von grundlegender Bedeutung sind. 
hat sich das Anwendungsgebiet dieser Disziplin stän- die’ Werken bone: 
dig erweitert und ist allmählich in fast alle Teile 
der reinen und der angewandten Mathematik ein- Probleme, die aus Raumer parnis weggwela (sen wer 
be- zedrung den mußten und enthalten viele AÄnrerungeen für 
In den letzten Jahrzehnten ist die Topologie. 
fußend auf den Arbeiten von H. Poineare und Das enthält ein sorrlältizes Literaturvei 
lungen zeiehnis und wird dureh viele Firuren, (die die An 
iträge auszedehnten Disziplin geworden, und es hat sich schauen 
ein riesiges Material von Ergebnissen angesammelt. 
‚chene das in den verschiedensten Zeitschriften zerstreut München, CCaratheodory. 524 


an die 


Nach auftauchen, sind beträchtlich: man muß die ab- 
ı man straktesten Schlußketten (ohne welche nur ein sofort Preis weh 
ufügt zerrinnendes Luftschloß entstehen würde) mit dem 
h lebendigsten AÄnschauungsvermögen verbinden, wel- Mit dem vorliegenden zweiten Bande liegt (die von 
ches allein imstande ist, die gewonnenen Resultate Frank und v. Mises herausgegebene Neubearbeitung 
zu deuten und den Anwendungen zuzuführen. Diese «des Riemann-Weber vollständig in zweiter Auflage 
komplizierte Aufgabe haben die beiden Autoren, die vor. Dabei sind die Änderungen und Erweiterungen. 
übrigens dureh ihre eirenen Untersuchungen im die der zweite Band erfahren hat. wesentlich um 
Gebiete der Topologie einen glänzenden Namen er-  fangreicher als die des ersten. Das ist dadureh be 
worben haben, restlos welöst. Auch die schwieri-  dingt. daß es nieht nur galt, durch die Änderungen 
zeren Teile des Buches können glatt gelesen werden, eine größere Einheitliehkeit der Darstellung anzu 
und fortwährend wird, an sorgfältig gewählten, stets streben. — auch hier ist es, wie oft in Büchern, die 
interessanten Beispielen, die Bedeutung der vor- von mehreren Verfassern geschrieben sind, dab «die 
getragenen Sätze reich illustriert. Ich bin der Über- zweite Auflage ihrem Ideal in diesem Punkte 
zeuzung, daß auch der Anfänger bei zeschiekter wesentlich näher kommt als (die erste, somdern 
Benutzung des Buches, d. h. wenn er (lie axioma- vor allem den Stoff der weiteren Entwieklung der 
tischen Schlußfolgen (die in Wirkliehkeit das Rück- — Physik anzupassen: denn der Frank-Mises ist kein 
grat der Topologie bilden) zunächst flüchtiger liest. Buch. das einen abzeschlossenen Stoff behandelt 
h auf sehr schnell in die Lage kommen wird. bei einem sondern ein Buch, das mit der Wissenschaft lebt. 
zweiten Studium des Buches, ein volles Verständnis Die wiehtieste Tatsache, der aus dem letzten 
auch für die komplizierteren Fragen zu erwerben. (runde Rechnung getragen werden mubte, ist die 
Die Das Buch beginnt mit einem einleitenden Kapitel.  Entwieklung der Wellenmeehanik seit «dem Kr 
in dem eine sehr eeschieckte Auswahl von Proble- scheinen der ersten Auflage. Nicht nur ist (des 
men, die in der Topologie auftreten. rein anschaulich wegen ein von Beck (Kansas) zeschriebener Ah 
behandelt wird. Im zweiten Kapitel wird der Begriff schnitt über die mathematischen Probleme der 
ler axiomatisch definiert, und die ein- Wellenmeehanik. der die Schrödingersche Gleichung. 
fachsten Bausteine der topologischen Räume, die ihre Eizenwerte, ihre wiehtigesten Lösungstvpen, (lie 
kann Simplexe sowie auch ihre Zusammensetzung zu  Störungstheorie und das Mehrkörperproblem sowie 
Venn Komplexen erklärt. In dritten Kapitel werden die (die Diraesche relativistische Wellengleichung hehan 
ersten Invarianten. d. h. die Homologiegruppen, die delt, angefügt worden. sondern es sind auch in 
0S y j Bettischen Zahlen. die Torsionskoeffizienten einge- anderen Abschnitten Umarbeitungen und Zusätze ge 
‚Bten führt. Daß dies Invarianten sind. d. h. daß sie von macht worden, die auf die Wellenmeehanik Rück 
ter speziellen Darstellune eines Komplexes unab- sicht nehmen. Um nur das Wiehtigste zu erwähnen. 
hängie sind. wird im folgenden Kapitel mit Hilfe ist in den Abschnitt über Mechanik zur Vorbereitun« 
ler Methode der simplizialen Approximation,. die ein Kapitel über die Strahlenoptik eingefügt und 
tich- \ Brouwer ersonnen und später verschiedene Auto- die Hamiltonsche mathematische Analogie zwischen 
ren vervollkommnet haben. bewiesen. Im fünften Lichtstrahl und Bahnkurve schärfer herauszearbeitet. 
Dies | Inapitel werden weitere Invarianten, wie z. B. Ferner ist in den Abschnitt über Wärmeleitung uni 
i ie Dimensionen eines Gebildes, untersucht. Hier- Diffusion ein Paragraph über «ie Probleme der 
auf werden (die zweidimensionalen Flächen be-  Wellenmechanik aufgenommen. die auf Differential 
trachtet, die einzigen Gebilde, die man bei dem  gleichungen vom Typus der Diffusionsgleiehung 
heutigen Stande der Topologie erschöpfend be- führen. An kleineren andere Gebiete berücksieh- 
Men, ; handeln kann. Die folgenden Kapitel (7. und 8.) tigenden Zusätzen seien erwähnt: ein kurzer Ab 
333 | sind dem Studium der Fundamentalgruppe eines schnitt über das Elektronenmikroskop und «ie ein- 


und infolgedessen schwer zugänglich ist. 

Es ist daher warm zu begrrüßen, daß Seifert 
und Threlfall ein Buch geschrieben haben. das 
bei nieht zu großem Umfang eine ausgezeichnete 
Übersicht über «die Methoden und über die Haupt- 
resultate dieses außerordentlich anziehenden Teils 
der Mathematik liefert. Die Schwierigkeiten, die bei 
ler Redaktion eines derartigen Buches immerfort 


Dr. PHILIPP FRANK, Prof. a. d. Deutschen Univ, 
in Prar. und Dr, RICHARD v.MISES, Prof. a. Univ. 
Istanbul, Die Differential- und Integral- 
eleichungenderMechanik und Physik. 
Bd. II Physikalischer Teil. 2. verm. Aufl. zuel. 
Ss. Aufl. von Riemann-Webers Partiellen Differential 
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eehendere Behandlung einiger Sätze der statistischen 
Mechanik. Darüber hinaus findet man überall Zu- 
sätze und Erweiterungen und wenn es nur Hinweise 
auf neuere Arbeiten sind, die die hier gegebene Dar- 
stellunz ergänzen. 

Dienen all diese Erweiterungen zum zuten Teil 
auch der gegenseitigen Anpassung der einzelnen Ab- 
schnitte aneinander, so wird die Einheitliehkeit der 
Darstellung vor allem dadurch wesentlich gefördert. 
daß die Bearbeitung der Hydromechanik der idealen 
und der zähen Flüssiekeiten einheitlich von v. Mises 
und Schulz durehzeführt ist. 

Natürlich ist dureh die erwähnten Erweiterungen 
der Umfang und damit leider auch der Preis des 
Bandes beträchtlich gewachsen. Das ist aber das 
einzige, was daran bedauerlich ist, denn alle Än- 
derungen sind wesentliche Verbesserungen, für die 
man den Herauszebern und ihren Mitarbeitern nur 
dankbar sein kann. Trotz aller Erweiterungen hat 
der Frank-Mises den Lehrbucheharakter gewahrt. 
und zwar den eines ausgezeichneten Lehrbuches, das 
aus der Literatur der anzewandten Mathematik und 
der theoretischen Physik nicht mehr fortzudenken 
ist. Auch «dem Leser dieser Zeitschrift. d. h. also 
dem wissenschaftlich interessierten wird 
das Buch in vielen Fällen ein guter Berater sein. 
Möge es auch in diesen Kreisen immer weitere Ver- 
breitune finden, 

Freiberg (Sa.) 


Fr. A. Willers. 537 


R. BECKER, o. Professor a. d. Techn. Hochseh. 
Berlin. Theorie der Elektrizität. (Neu: 
hbearb. des Werkes von M. Abraham.) Band I. Ein- 
führung in die Maxwellsche Theorie der Elektri- 
zität. mit einem einleitenden Abschnitte über das 
Rechnen mit Vektorgrößen in der Physik, 10. Aufl. 
VI #265 S. m. 59 Fig. Leipzig und Berlin 1933. 
Verlag B. 6. Teubner. Preis geb. 14.50 M. 


Es sprieht für die Brauchbarkeit der Neubearbei- 
tung des Abrahamschen Buches durch Becker, dab 
schon nach so kurzer Zeit eine neue Auflage nötig 
geworden ist. Gegenüber der neunten im 13. Jahr- 
eane 10933 dieser Zeitschrift von F. Noether be- 
sprochenen Auflage sind wesentliche Änderungen 
nieht vorgenommen worden. Im einzelnen sind je- 
loch vielfache Erzänzungen so ist z. B. ein Pa- 
ragraph über den Entmagnetisierungsfaktor einge- 


fügt aber auch Vereinfachungen und Kürzungen 
— z. B. in dem Paragraphen über Drahtwellen bei 
endlichem Widerstand auf Leitungen — erfolgt. 


Eine nochmalige eingehende Empfehlung des aus- 
eezeichneten Lehrbuches erübrigt sieh. 


Fr. A. Willers. 504 


Freiberg (Sa.). 


X. HAAS, Prof. für Physik an der Universität 
Wien. Kleiner Grundriß der theore- 
tischen Physik. Kleine, besonders bearbeitete 
Ausgabe der Einführung in die theoretische Physik. 
VII + 183 8. mit 22 Fire. Berlin und Leipzig 1934. 
Verla Walter de Gruyter u. Co. Preis geb. 5.30 M. 


Der kleine Grundriß ist eine Parallel-Ausgabe zu 
der bekannten, viel gebrauchten „Einführung in die 
theoretische Physik“ des gleichen Verfassers, hat 
aber nur den sechsten Teil des Umfanges. Trotz- 
dem gibt er einen Überblick über das gesamte Ge- 
biet der theoretischen Physik unter besonderer Be- 
rücksiehtieunge der modernen Probleme. Der ge- 
ringe Umfang bedingt natürlich, daß die Ableitungen 
und Anwendungen meistens fehlen, so daß in man- 
ehen Gebieten nur das Gerüst gereben wird. Um 
hier wirklich zu einem Verständnis zu kommen, 
wird daher die Heranziehung ausführlicherer Dar- 
stellungen erforderlich sein. Aber infolge der ge- 
schiekten Auswahl des Stoffes und bei der außer- 
ordentlichen Kunst des Verfassers, schwierige 


Dinge so darzustellen, daß sie den Anschein der 


Einfachheit gewinnen, dürfte das vorliegende Buch 
für den ganz zeeigmet sein, der einen ersten Über- 
blick gewinnen oder der vorhandene Kenntnisse 
auffrischen will. 

Freiberg (Sa.). 


Fr. A. Willers. 504 


KARL MARBE, Grundfragen der Wahr- 
scheinliehkeitsreehnung und theore- 
tischen Statistik. V + 177 München 1934, 
Verlaxe C, H. Beck. Preis 8 RM. 

Dort, wo die praktische oder grundsätzliche Unzu- 


länglichkeit unserer Beobachtungen der kausalen Er- 
eründung des Naturzeschehens unüberwindliche 


Schranken setzt, hat sich die statistische Unter- 
suchungsweise eingebürgert und in weitem Um- 


fanze auch bewährt. Das mathematische Rüst- 
zeug (dieser Methode bildet die  Wahrschein- 
liehkeitsreehnung: ihr Hauptbegriff, die Wahr- 
scheinliehkeit, wird durch eine ideale  Ope- 
ration erklärt, nämlich durch die Erhebung 
einer unendlichen Folge gleichartiger Einzelbeob- 
achtungen. Da eine statistische Erhebung nur end- 
lieh viele Beobachtungen erfassen kann, ist man 
eezwungen, bei den Anwendungen die idealen 
Wahrscheinlichkeiten dureh empirische Häufigkeiten 
zu ersetzen: die hieraus zu erwartenden Ab- 
weichungen lassen sich aber nach den Lehren der 
Wahrscheinliehkeitsreehnung abschätzen. Dadurch 
öffnet sich ein Wer zur Prüfung der Frage, 
ob die Beobachtungen mit den Folgerungen der 
Wahrscheinlichkeitsreehnung hinreichend im Ein- 
klange sind. 

Der Verfasser bestreitet das für eine Reihe von 
Erscheinungen an der Hand statistischer Stoff- 
sammlungen, die mit Sorgfalt zusammengetragen 
sind. In der Tat zeiren sich hier mancherlei auf- 
fallende Abweichungen: sie erklären sich aber zu 
einem Teile aus der Beerenztheit des Stoffes, zu 
einem anderen daraus, daß die Erhebungsart die 
Forderungen der Wahrscheinliehkeitsreehnung nur 
unvollkommen erfüllt. Keinesfalls bilden sie eine 
tragfähige Unterlage für Marbes „Lehre vom stati- 
stiichen Ausgleiche“, wonach die empirischen 
Häufickeiten schneller gegen die ideale Wahr- 
scheinliehkeit konvergieren sollen, als es nach der 
Wahrscheinliehkeitsreehnunge zu erwarten ist. — 
Das Buch liefert dureh Stoffsammlung und mathe- 


matische Analyse einen wertvollen Beitrag zum 
Problem der Iterationen. 
Dresden. P. E. Böhmer. 497 


Dr.-Ing. E. h. WILHELM NUSSELT VDI, o. 
Prof. a. d. Techn. Hochsch. München, Tech- 
nische Thermodynamik. (Sammlg. Göschen 
3d. 1084.) 144 S. m. 65 Abb. Berlin und Leipzig 
1934. Walter de Gruyter & Co. Preis geb. 1,62 M. 


Wenn heute eine Technische Thermodynamik im 
Umfange eines Göschenbändchens erscheint. so wird 
man nieht fehlzehen in der Annahme, daß darin der 
bekannte Weg eingeschlagen wird. der dureh die 
Namen Carnot und Clausius gekennzeichnet 
ist, und daß sein Ziel der II. Wärmehauptsatz und 
dessen technische Anwendungen sind. Um so mehr 
muß man es beerüßen, wenn ein führender For- 
scher wie Prof. W,Nusselt, München, sich bereit 
findet, den Band zu verfassen. Denn die richtige 
(sewiechtsverteilung über alle Einzelheiten beim 
Aufbau und die treffendsten Beispiele für jede 
Stelle stehen nur dem zu Gebote, der als Forscher 
immer wieder kritisch an neue Aufgaben herangehen 
muß, wie sie sich bei der Fortentwicklung ergeben. 
Daher wird nieht nur derjenige Nutzen aus dem 
Buche ziehen, der es zu seiner eigenen Belehrung 
dureharbeitet, sondern auch jeder Lehrer der Tech- 
nischen Thermodynamik wird sich die eine oder 
andere Wendung merken. (Nur beim ersten Bei- 
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spiel für nichtumkehrbare Vorgänge auf Seite 27/28 
hätte man dem anspruchsvollen Leser deutlicher 
saren können, daß nicht die Beschleunigung, son- 
dern die Zähigkeit des Gases einen nichtumkehr- 
baren Druckanteil hervorruft.) In der Auswahl des 
Stoffes mußte natürlich dem beschränkten Umfang 
Rechnung getragen werden. So widmet sich der 
jand zu nahezu gleichen Dritteln den Vorbereitun- 
ven (1. Hauptsatz, technische Arbeit, Energiestrom 
und besondere Zustandsänderungen), dem II. Haupt- 
satz und dem technisch so außerordentlich wichti- 
ven Dampfzustand (nasser und überhitzter Dampf). 
Gerade in dem letzten Abschnitt sind die älteren 
und die neueren Zustandsgleiehungen für Wasser- 
dampf bemerkenswert ausführlich behandelt. Bei 
einer Technischen Thermodynamik braucht kaum 
erwähnt zu werden, daß den Arbeits-, Wärme- und 
i-s-Diagrammen geebührender Raum gegeben ist. Die 
zahlreichen Abbildungen und der bekannte auch in 
den Formeln saubere Druck machen den Band an- 
venehm lesbar. 

Dresden. A. Busemann. 498 


Dr.-Ing. HUGO RICHTER VDI, Priv.-Dozt. a. dl. 
Bergakademie Freiberg, Rohrhydraulik. All- 
vemeine Grundlagen, Forschung, Praktische Be- 
rechnung und Ausführung von Rohrleitungen. IX 
+ 256 5. m. 192 Textabb. u. 44 Zahlentaf. Berlin 
1934. Verlag von Julius Springer, Preis geb. 
22,50 M. 


Es war ein recht dankenswertes Unternehmen, 
aus der umfangreichen neuen Literatur der Hydrau- 
lik dasjenige zusammenzustellen, was sieh auf die 
stationäre Strömung in Rohrleitungen bezieht. Der 
Verfasser hat ohne viel eigenes hinzuzutun die Er- 
gebnisse gesammelt, bis einschließlich zu denen der 
neuesten Untersuchungen vonPrandtl.Kärmän 
und Nikuradse. Das Buch ist gegliedert in einen 
theoretischen Teil, der die alleemeineren Sätze und 
Versuchsresultate enthält, und einen praktischen, 
der die Rechnungsvorschriften für die praktische 
jenutzung, auch unter Verwendung von Nomo- 
srammen gibt. Behandelt wird hauptsächlich das 
gerade Rohr von Kreisquerschnitt, doch sind auch 
Andeutungen für den Fall anderer Querschnitte. 
sowie für Rohrkrümmungen zereben. Das reich- 
haltige Zahlenmaterial scheint einwandfrei zu sein. 
Das Buch wird sich gewiß in den Kreisen der In- 
genieure, die mit der Herstellung von Leitungen zu 
tun haben, bewähren. Mises. 503 


Dr.-Inz.O.K.FRÖHLICH, Druckverteilung 
im Baugrunde mit besonderer Berücksich- 
tigung der plastischen Erscheinungen. 183 S. m. 
68 Textabb. Wien 1934, Verlag Julius Springer. 
Preis kart. 15 M. 


Der Verfasser stellt sich in dem vorliegenden 
suche die Aufgabe, dem ausführenden Ingenieur 
jene Unterlagen zu vermitteln, welche zur Be- 
stimmung der Spannungsverteilung und der Form- 
änderungen des Baugrundes notwendig sind. Es ist 
dem Verfasser mit Erfolg darum zu tun gewesen, ein 
möglichst abgeschlossenes Bild des heutigen Standes 
der bisher auf diesem Gebiete errungenen Erkennt- 
nisse zu geben, welches durch eigene Unter- 
suchungen und Gedanken erweitert, dem Leser die 
physikalisch-mechanische Seite des Problems näher- 
bringen soll. Inhalt und Anordnung des Stoffes. 
nämlich zuerst die Spannungsverteilungen im ela- 
stischen und plastischen Bereich, sodann die Form- 
änderungen zu behandeln, weiterhin die vorhandenen 
Versuche kritisch zu besprechen und schließlich aus- 
zeführte Frndierungen als Beispiele für Unter- 
suchungz und Nachreehnung heranzuziehen, lassen 
dieses Buch, das von einem Ingenieur für Ingenieure 
zeschrieben wurde, als einen sehr brauchbaren Leit- 


faden für den Praktiker erscheinen. Aber auch der 
Theoretiker wird aus demselben eine Reihe von An- 
regungen schöpfen und auf die Bearbeitung von 
Problemen geführt werden, deren Lösung für den 
Grundbau von Wichtigkeit wäre. Denn die Theorie 
hat bislanz eine Reihe von Fragen, die von Bedeu- 
tung für die Praxis sind, unbeantwortet gelassen, 
und wenn deshalb der Verfasser in solehen Fällen 
zu Näherungslösungeen gegriffen hat, die einer 
strengen Begründung nieht standhalten, so darf ihm 
hieraus zumindest in jenen Fällen kein Vorwurf ge- 
macht werden, wo er den Leser auf diesen Umstand 
aufmerksam macht. Andere Fälle, wie z. B. die Ein- 
führung eines verfüzbaren Parameters um die 
Radialspannungen im unendlichen Halbraum bei Be- 
anspruchung durch eine Einzellast in der Form 0, = 
z cos dr r ? darstellen zu können, sowie die Fol- 
eerungen, die hieraus zum Teil gezogen werden, be- 
dürfen wohl noeh einer Klärung dureh die strengere 
Theorie. 

Wien. Ernst Melan. 508 


Dr.-Ing. WILHELM FLÜGGE, Priv.-Doz. a. d. Un. 
Göttingen, Statik und Dynamik der Scha- 
len. VIL + 240 S, m. 95 Textabb. Berlin 1934, Ver- 
lag von Julius Springer. Preis geh. 21 M. 


Das Schrifttum über die Theorie der Schalen mit 
und ohne Biegungssteifigkeit ist in den letzten zwei 
Jahrzehnten so angewachsen, daß es kaum mehr zu 
übersehen ist, auch sind viele wiehtige Aufsätze 
schwer zueänglich. Bei den auf den verschiedenen 
Gebieten der Technik immer mehr zunehmenden 
praktischen Anwendungen der Schalentheorie ist 
weiter die Vielgestaltiekeit der einschlägigen Ar- 
beiten, die in ihren Voraussetzungen auch sehr ver- 
schieden sind, ungemein störend. Hier will das vor- 
liegende Buch dureh eine übersichtliche und zusam- 
menfassende Bearbeitung auf einheitlicher Grund- 
lage Abhilfe schaffen. 

Es sei vorwereenommen, daß dies dem Verfasser 
voll zelungen ist. Das Buch ist eine Bereicherung 
unserer Fachliteratur, von der jeder auf dem Ge- 
biete der Schalentheorie arbeitende Ingenieur oder 
Forscher Nutzen haben wird. 

Das Buch gibt nach einer knappen Darstellung der 
alleemeinen Grundlagen der Theorie der dünnen 
Schalen die vollständige Theorie der Membranen. 
Es werden behandelt: die besonders für den Be- 
hälter- und Kuppelbau wichtigen Rotationsschalen 
unter symmetrischer untl unsymmetrischer Bela- 
stung, die Zylinderschalen, die in letzter Zeit z. B. 
bei den dünnen Eisenbeton-Tonnendächern vielfache 
Anwendung gefunden haben, weiter die Vieleck- 
kuppeln- und Translationstlächen. Den Abschluß der 
Ausführungen über die Schalen ohne Biegungsstei- 
firkeit bildet eine Betrachtung der aftinen Membran- 
schalen. 

Sodann folgen die entsprechenden Untersuchungen 
für biegungssteife Kreiszylinderschalen (Flüssigkeits- 
behälter, Tonnendächer), Rotationsschalen (Kugel-, 
Kerel- und beliebige Rotationsschalen), sowie die 
Theorie der Faltwerke mit und ohne Biegungssteifig- 
keit. 

Ganz besondere Beachtung verdient der Abschnitt 
über den heutigen Stand der Theorie der Kniekung 
von Schalen. Es werden betrachtet: Kreiszylinder 
unter Normalspannungen und Torsion, die Ursachen 
des Bierungsbruches von Schalen, sowie Kugel- 
schalen unter radialem Druck. 

Zum Schluß wird noch eine kurze Darstellung der 
Theorie der Schwingungen drehsymmetrischer Scha- 
len zereben, die an die grundlegenden Arbeiten von 
Rayleigch anknüpft. Die gut ausgewählte Lite- 
raturübersicht dürfte vielen erwünscht sein, 

Aus vorstehender Übersicht erkennt man den rei- 
chen Inhalt des Buches von Flügge. Ein solches 
klares, knappes und doch gut verständliches Buch, 


ipzig | 
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in dem die vielen Originalarbeiten auf dem Gebiete 
(ler Sehalentheorie in einheitlicher Weise verarbeitet 
sind, wurde von Praxis und Forschung schon längst 
gewünscht. Der Verfasser blieb jedoch nicht bei 
der zusammenfassenden Bearbeitung stehen, sondern 
eab auch viel Eigenes. Man findet fast in allen Ab- 
schnitten wesentliche neue Beiträge des Verfassers 
zur Schalentheorie. Das Buch wird sich daher bald 
als unentbehrlich erweisen. 


Hannover, Ferd. Sehleieher. 509 


Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 


bleibt vorbehalten): 


Dr.-Inz. PAUL WERNERS, Enereie-Über- 
tragunge und -Umwandlung mit Wech- 
selstrom. 204 8. Leipzie und Berlin 1055, Ver- 
lag B. 6. Teubner. Preis geb. 18 M. 


Dr.-Ing. F. MOELLER und Dipl.-Ing. TH. WERR, 
Llektrotechnik. (Ein Leitfaden für Studium 
und Praxis, hrsg. von Dr.-Ing. G. Bolz, Berlin, Dr.- 
Ing. F. Moeller. Berlin, Dipl.-Ing. Th. Werr,. München.) 
Bd. II, Teil 3 und 4. Gleieh- und Wechsel- 
strommaschinen,. VIII + 125 S. m. 176 
Leipzig und Berlin 1955. Verlag B. G. Teubner. 
Preis geb. 10 M. 


Dr. phil. LUDWIG HÄNERT, Marine-Studienrat in 
Mürwik,. Geschützund Schuß. Einführung in 
die Gesehützmeehanik und Ballistik. 2. verb. Aufl. 
VI #+370 8. m. 161 Textabb. Berlin 1935, Verlag 
Julius Springer. Preis geb. 27 M. 


W, KRULL, Iılealtheorie, (Ergebnisse (der 
Mathematik und ihrer Grenzgebiete, herauszereben 
von der Sehriftleitune des .Zentralblatt für Mathe- 
matik“. 4. Bd.. 3.) V+152S8. Berlin 1935. Verlag 
Julius Springer. Preis brosch. 17.50 M. 


LAMPE, Studiendirektor der Oberrealschule 
.isterschloß in Elsterwerda. und Dr. A. WAGNER, 
Studienassessor in Jena. MathematikundWehr- 


sport.  (Mathematisch-physikalische Bibliothek. 
teihe I. Bd. 80.) 51 8. m. 20 Fig. Leipzig uni 
Berlin 1935. Verlage B. G. Teubner. Preis kart. 


1.20. M. 


ERNST AUGUST WEISS, a. 0. Prof. a. d. Rhein. 
riedrieh-Wilhelm-Universität Bonn. Einführung 
indie LinienzgeometrieundKinematik. 
(Teubners Mathematische Leitfäden, Bd. 41.) 
VI+122S. Leipzig und Berlin 1935, Verlag B. G. 
Teubner. Preis kart. 7.60 M. 


Dr. WERNER HEINTZE, Studienrat a. d. Hinden- 
burg-Oberrealschule in Eberswalde Kristall- 
projektion im Vergleich mit entspre- 
chenden Erdkarten und mit einer Än- 
wendune auf die Laue-Aufnahmen. 
Mathematisch-physikalische Bibliothek, Reihe 1. 
Bd. 82.) IV m. 27 Abb. Leipzig und Berlin 
1434, Verlag B. G. Teubner, Preis kart. 1.20 M. 


Dr.-Ing. Dr. E. h. WILHELM NUSSELT, 0. Prof. 
a. d. Teehn. Hochseh, München. Technische 
Thermodynamik. (Samml. Göschen Bd. 1084.) 
144 S. m. 65 Abb. Berlin und Leipzie 1934, Verlag 
Walter de Gruyter & Co. Preis geb. 1.62 M. 


FUCHS-HOPF-SEEWALD, \erodynamik. 
Bd. 2. 2. Aufl. Theorie der Luftkräfte von 
R. Fuchs. VITL + 310 S. 224 Abb. und einem An- 
hang. Berlin 1935. Verlag Julius Springer. Preis 
M. 


Dr. GEORG JOOS, 0. ö. Prof. a. d. 


hehrbuceh der Theoretischen 


Univ. ‚Jena. 
Phvsik. 


2. Aufl. XVI+ 676 S. m. 164 Fig. Leipzig 1934, 
Verlag Akademische Verlarszesellschaft m.b.H. 


Preis 24 M. 


NACHRICHTEN 


Martin Grübler 

Am 31. Mai starb in Dresden der emeritierte Pro- 
fessor an der Teehnisehen Hochschule, Geh. Hofrat 
P’rof. Dr, phil. h. e.. Dr. art. ing. h. e. Martin 
(‚rübler. (‚rübler wurde am 19. Dezember 1851 in 
Meerane in Sachsen zeboren, wo sein Vater Apo- 
tleker war. Nachdem er die höhere Gewerbeschule 
in Chemnitz mit Auszeichnung absolviert hatte, be- 
zor er die Teehnische Hochschule in Dresden, um 
von 1870-73 daselbst Bauingenieurwesen zu stu- 
‚lieren. Nach einer zweijährigen praktischen Tätig- 
keit führte ihn der Wunsch. sieh dem Lehrfache zu 
widmen, an die Hochschule zurück. Er stwlerte 
von 1875-1880, seinen Lebensunterhalt durch Er- 
teilunz von Privatstunden verdienend, in Dresien 
und in Leipzie Mathematik und Naturwissenschaften 
und bestand im Sommer 1880 die technische Ober- 
lehrerprüfunge in Dresden. Im Herbst desselben 
‚Jahres zinz er auf Veranlassunz von Zeuner nach 
Zürich. wo er sieh habilitierte und wo er gleich- 
zeitig eine Assistentenstelle bei Herzog (Maschinen- 
technik) und Flieener (Mechanik) bekleidete In 
(diese Zeit fallen seine ersten selbständigen wis- 
sensehaftliehen Arbeiten auf dem Gebiete der Kine- 
matik, die teils im Zivilingenieur, teils in der Zeit- 
schrift für Mathematik und Physik erschienen. Die 
anerkennende Beachtung, die diese Arbeiten fanden. 
drückte sieh darin aus, daß Grübler, nachdem er 
sich im Herbst des Jahres 1885 nach Dresden um- 
habilitiert hatte, im Frühjahr 1886 als Ordinarius 
für Mechanik und Festigekeitslehre an das Polytech- 
nikum nach Riga berufen wurde. Dort hat er bis 


zum Herbst (es Jahres 1895 zewirkt: dann veran- 
laßten ihn die politischen Zustände, einer Aufforde- 


runz Müller-Breslaus folgend nach Berlin über- 
zusiedeln,. wo er bis zum Jahre 1900 Mechanik und 
Festiekeitslehre für Bauingenieure vortrug. Als im 
‚Jahre 1900 das ihm versprochene Ordinariat endlich 
errichtet werden sollte, bekam er gleichzeitir einen 
Ruf, als Nachfolger von Otto Mohr, an die Tech- 
nische Hochschule in Dresden. Er folgte (diesem 
Ruf und hat bis zu seiner Emeritierung 1920 an der 
Dresdner Hochschule zewirkt. 

Eine große Anzahl schöner Publikationen leren 
von der Begeisterung Zeugnis ab. mit der Grübler 


(len wissenschaftlichen Problemen nacheing. Neben 
den schon erwähnten Arbeiten zur Kinematik be- 
schäftieten ihn schon von Anfang an Unter- 


suchungen zur Strömungslehre, die in dem Referat 
über hydraulische Maschinen in (der Enzyklopädie 
der mathematischen Wissenschaften mit verarbeitet 
sind. Anrerungen. die aus seiner Lehrtätiekeit ent- 
sprangen, führten ihn später auf andere Gebiete der 
teehnischen Mechanik. Sehr bekannt geworden sind 
seine Arbeiten zur Festiekeitslehre, insbesondere die 
Untersuchungen über die Festiekeit rotierender 
Scheiben. 

Eine vierzigjährige, fruchtbare Lehrtätigkeit an 
vier Hochschulen. die in seinem Lehrbuch der tech- 
nischen Mechanik ihren Niederschlag gefunden hat. 
hat ihm die Dankbarkeit zahlreicher Schüler einge- 
tragen. Als äußere Anerkennung seines unermükd- 
lichen Wirkens ernannte ihn die Universität Gieben 
1022 zum Dr, phil. h. e.. die lettländische Universität 
zu Riga im Jahre 1929 zum Dr. art. ine. h. e. Die 
(sesellschaft für angewandte Mathematik und 
Mechanik wird dem treffliehen Manne ein ehrendes 
(edächtnis bewahren. Trefftz. 526 
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Gesellschaft für angewandte Mathematik 
und Mechanik. 


Hauptversammlung in Stuttgart. 


Die Gesellschaft für angewandte Mathematik und 
Mechanik hielt ihre ordentliche Hauptversammlung 
und wissenschaftliche Tagung vom 26. bis 28. Sep- 
tember in Stuttgart ab. und zwar im Rahmen des 
11. deutschen Mathematiker- und Physikertages, der 
vom 22. bis 29. September «daselbst stattfand. — Es 
wurden die foleenden Vorträze gehalten: 


Donnerstag, 26. 9: 1. W. Müller. Aachen: An- 
laufströmunge im Kreisring und Kreisrohr. 2. 0. 
Weber, Dresden: Kanalwellen bei geringer Wel- 
Ienhöhe. 3. M. Schuler, Göttingen: Erzeurzung 
von Oberflächenwellen durch schwingende Körper. 
Weber, Dresden: Theoretischer Ansatz für 
lie Sehulersehen Umsehlagwellen. 


Freitag. 27. 9:1. G. Weinblum.,. Berlin: Mo- 
Jellversuche zur Ausbildung günstigster Schiffs- 
(formen. 2. B. Eck. Köln: Experimentelle Behand- 
lung der Strömungslehre. Abresehen von diesen 
beiden Vorträgen war (dieser freigehalten für 
den Besuch der Vorträge über Schwineungsfragen, 
die von «den Physikern gemeinsam mit dem VDI 
veranstaltet wurden. 


Sonnabend, 28. 9: 1. Th. Pöschl. Karlsruhe: 
Über eine direkte Methode zur Lösung von Auf- 
oaben der Plastizitätstheorie. 2. G. Vorelpohl., 


ZUSCHRIFTEN AN 


über den Heaviside-Kalkül. Icrr K. Th. 
Vahlen zitiert in seiner Arbeit ..Über den 
Heaviside-Kalkül“. Band 13. Heft 4, meine Arbeit 
„Über die Lösung von Schwinzungsaufraben mit- 
tels symbolischer Differentialreehnung“, Archiv für 
Elektrotechnik. Band XXIV. Seite 360-382, und 
„war nur insofern, als er an zwei, übrigens nicht 
sehr belangreichen Stellen Unstimmigkeiten zu er- 
kennen Der Zweck der vorliegenden 
Zeilen ist es. zunächst die volle Korrektheit dieser 
beiden Details nachzuweisen. 

Die erste Beanstandung Herrn Vahlens 
(Seite 2094, Fußnote 6 seiner Arbeit) bezieht sieh 
auf die Ableitunz der bekannten (speziellen) 
Heaviside-Gleiehung (entsprechend Gl.  (54b). 
Seite 381 meiner Arbeit): 


n 


r \ 
— Ay Z’ (Am) 


aus der von mir erstmaliz webrachten Verallge- 
meinerung der Heaviside-Gleiehung (entsprechend 

h): 


1 
— Il’ (a, — Z (Gy) 
x 


Die Richtigkeit der beiden Gl. (53b) und (54 b) 
steht nieht zur Diskussion. entsprechen 
übrigens den Gl. (XD (©) bzw. (XII) (6) 
ler Arbeit Vahlens, und es kann sich also nur 
darım handeln, ob die Spezialisierung der allge- 


Berlin: Zur hydrodynamischen Theorie der Lager- 
reibune. 3. H. Neuber, München: Randbedin- 
gungen der Hauptspannungs-Stromlinien. 4. K. 
Hannover: Kniekung von Rechtecks- 
platten mit streifenweise konstanter Dieke. 5. F. 
Schultz-Grunow. Kassel: Anschauliche Scha- 
lentheorie. 6. K. Margzuerre,. Berlin: Thermo- 
elastische Plattengzleiehungen. 7. W. Flügge, 
Göttingen: Bemerkungen zur Theorie statisch un- 
bestimmter Systeme, 


Vor den Vorträgen fand Sonnabend, 28. 9, unter 
lL.eitunze von Herrn Prandtl die ordentliche 
Hauptversammlung der G.A.M.M. statt. Nachdem 
der im Laufe des letzten Jahres verstorbenen Mit- 
elieder ehrend gelacht worden war, erstäattete (der 
Geschäftsführer, Herr C. Weber, Dresien. Bericht 
über Mitgliederzahl und Kassenverhältnisse, Dem 
Vorstand wurde Entlastung erteilt. Gewählt wurden: 
in den wissenschaftlichen Ausschuß an Stelle der 
satzuneseemäb  ausscheidenden Herren Hamel 
(Berlin) und Geekeler (Kiel) die Herren Pröll 
(Hannover) und Funk (Prag) zu Kassenprüfern 
die Herren Betz (Göttineen) und Flügrwe (Göt- 
tineen). Der ausführliche Jahresbericht wird den 
Miteliedern der G.A.M.M. in der üblichen Form 


Vorträge und Hauptversammlung waren erfreulich 
reich besucht. Daß auch der außerwissenschaft- 
liche Teil der Tagung dank der schwäbischen Gast- 
freundschaft einen harmonischen Verlauf nahm. sei 
besonders bemerkt. 


DEN HERAUSGEBER 


meinen Formel (53 b) auf die Formel (54 b) in rieh- 
tiger Weise vorgenommen wurde. Um nun jedes 
Mißverständnis auszuschließen. sei das Wesent- 
liehe der in meiner Arbeit in ein wenig knapper 
Form gebrachten Ableitung mit etwas anderen 
Worten ausführlicher wiederholt: insbesonders sei 
der Übersichtlichkeit halber bei der in meiner 
Arbeit verwendeten, allgemeinen  Betrachtungs- 
weise der „Integration über die Wellenstirn“ eine 
einfache spezielle Form der Störungsfunktion wie 
folgt zugrunde gelezt. 

Spezialisieren wir die Störungsfunktion F (r) 
in Gl. (53b) in der Weise, daß für 2 
verschwindet, für um (den Betrag eins 
springt und für = >.r* den konstanten Wert eins 
besitzt (wobei wir „* >0 voraussetzen wollen), 
so kann man für die Lösung f(r). da anfänglich 
das System in Ruhe ist und bis zum Zeitpunkt 
„niehts geschieht”. sicherlich schreiben: 


d. h. die Lösung ist eine Hauptlösung. Ebenso 


eilt voraussetzungszemäß: 
......... 
61. 53b) läßt sich nun sehr vereinfachen. da die 


Koeffizienten und sämtlich verschwinden, 
und man erhält sofort 


n x 
1 
(ra) =m N X | 
(Am, 
1 0 
n 
= 
Z (— 
1 
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Man erkennt. daß /(r) an der Stelle r = r* 
stetig ist. Dies gilt jedoch im allzemeinen nicht 
für die (n—1) ersten Ableitungen. da 

ist, während (r) 

x<x* 
für ein beliebizes r* im allzemeinen von Null 
verschieden ist. 

Im nun auf den Fall des Heavisideschen 
Einheitsstoßes E(r) überzuzehen. lassen wir .r* 
reren Null konverzieren. Dann ist f(r) für 
r >0 dureh Gl. (54b) bestimmt. Man hat dabei 
zu berücksichtiren,. daß 


\ 1 1 
au) 


1 
ist wie aus Gl. (48) 
m \ — 648 
‚Zi (Am) (Y — Am) 
1 


ohne weiteres hervorgeht, wenn man dort y 0 
setzt. In meiner Arbeit ist nun in diesem Zusam- 
menhang infolge eines Druckfehlers an Stelle der 
(Gl. (48) die Gl. (28) zitiert, und ich halte es für 
möglich, daß gerade dieser Druckfehler zu einer 
mißverständliehen Auffassung meiner Ableitung 
führen kann. Ich zlaube aber. daß aus dem hier 
venau durchgeführten Reehnunzszanz ohne wei- 
teres ersichtlich ist. daß in Wirklichkeit nur 
(1. (48) benützt wird. welche selbstverständlich 
auch für gebrochene Funktionen eilt. 

sildet man aus Gl. (S4b) die (n—1) ersten 
Ableitunzen und setzt man 0, so findet man, 
daß sie im allzemeinen nicht verschwinden, d. h. 
fir) entsprechend (54b) ist keine Hauptlösunge. 
Dies ist jedoch kein Widerspruch zur ursprüng- 
lichen Annahme, da wir den Unstetirkeitspunkt 
„even Null gehen ließen. somit /(r) (0) 
bereehnet aus (54 b) Cr) bedeutet und 

N D>N 
nach dem früher Gesagten im allgemeinen von 
Null verschieden ist, während lim") voraus- 

setzuneszemäß Null ist. Die sind 
(erensatz zur hier vorliexenden Haeviside-Lösung 
nur bei der später zu besprechenden Casper- 
sehen Lösung stetir. 

Sehließlich möchte ieh noch bemerken, dab «die 
Spezialisierung der Gl. zu Gl. (54 b) selbst- 
verständlich auch in anderer Weise vorgenommen 
werden kann. «loch ersehien mir «die gebrachte Ab- 
leitunz durch „Integration über die Wellenstirn“ 
rerade für «den Elektrotechniker besonders an- 
sehanulich. 

Der zweiten Einwendungz des Herrn Vahlen 
(Seite 2097 seiner Arbeit) liert folzender Sachver- 
halt zuerunde: Ich zeige auf Seite 376 meiner 
Arbeit. daß die Lösunz von Casper (Archiv für 
Elektrotechnik 1925, Seite 95) die Eirenschaft hat 
daß die 1) ersten Ableitungen zur Zeit = 0 
verschwinden. d. h. es ist hier lim 

Zu dem selben Resultat kommt Herr Vahlen 
(Seite 297 seiner Arbeit). Ich bemerke nun. dab 
im alleemeinen Fall des Heaviside schen 
Problems (Beschreibung «der Vorzänzre in einem 
mehrmaschigen Netz beim plötzlichen Anlegen 
der Einheitsspannung an das früher in Ruhe be- 
tindliche System) die Caspersche Lösung nicht 
herangezogen werden darf. Es ist nämlich f” (0) 
ein Maß für die anfänzlich in einer Masche wir- 
kende Spannung, und diese wird sowohl bei ein- 
fachen als auch bei zusammengesetzten Netzen 
mindestens bei einer Masche von Null verschieden 
sein. Dies steht, wie ich selbst ausdrücklich 
hervorhebe. in keinem Widerspruch zur Casper- 


schen Lösunz für den Fall eines einfachen Strom- 
kreises, da hier n = 2 jst und die obenzenannte 
Bedinzunz also nur verlanzt, daß die Funktion 
selbst und der erste Differentialquotient nicht aber 
der zweite Differentialquotient für 2 = ver- 
schwindet. Für den allgemeinen Fall hinzeren ist 
n > 2, daher verschwindet bei der Casperschen 
lösung auch der zweite Differentialquotient, wäh- 
rend bei der Lösunz des Heavisideschen 
Problems, wie erwähnt. der zweite Differential- 
quotient mindestens bei einer Masche von Null 
verschieden sein muß. 

Herr Vahlen brinet nun auf Seite 297 den 
Fall eines einfachen Stromkreises und zitiert 
überdies das Buch Carsons (Elektrische Aus- 
rleichsvorzänge und Operatorenreehnung, erw. 
deutsche Ausgabe von F. Ollendorf und 
K. Pohlhausen. Springer 1929) Seite 32, wo 
ebenfalls das Beispiel eines einfachen Stromkreises 
behandelt wird. Er hebt hervor. daß hier der 
zweite Differentialguotient für 70 nicht ver- 
schwindet und bemerkt hierzu. daß dies auch 
nicht bei der Casperschen Lösung (für n = 2) 
vefordert wird. Es bestehe hier also (im Gegen- 
satz zu meiner Behauptung) keinerlei Widerspruch. 
Es scheint also Herr Vahlen zu übersehen, daß 
ieh. wie in meiner Arbeit ausdrücklich angegeben. 
an der beanstandeten Stelle den allzeemeinen 
Fall (a>2) behanllle. 

Ich möchte die Gelegenheit benützen, um fest- 
zustellen, daß der Beweis (der Heaviside-Gleichun:z 
auszehend vom Ansatz /(D)y g(D)z. also un- 
abhängeie von ‚der bekannten  funktionstheore- 
tischen Betrachtunz K. W. Wazners. meines 
W'ssens erstmaliz in meiner erwähnten Arbeit 
ancereben wurde, wobei sich auch eine Verallge- 
meinerunz der Heavisideschen Lösung auf 
einen beliebieen Anfangeszustand des Systems bei 
Anlezen einer Spannung beliebiger Form (ent- 
sprechend den Formeln (53) bzw. (53a) und (53b) 
auf Seite 350 und 381 meiner Arbeit) ergibt. Ich 
halte diesen Hinweis für eine wünschenswerte Er- 
eänzun®e der Zitate meiner Arbeit durch Herrn 
Vahlen,. welehe «die Prioritätsfrage in keiner 
Weise berühren. 

Wien. W.Gauster 591 


Erwiderung. Die Differentialgleichung (Bd. 13. 
S. 293) 


hat die zu r=0 gehörige Hauptlösung 


1 \ | 


Die Gleichung 


«+ 


hat die Lösung 


(0) 
’ 
Dies ist nicht die Hauptlösung, denn sie ver- 
schwindet zwar für aber die n— 1 ersten 
Ableitungen werden für r=0 gleich 


und diese Summen verschwinden nur sämtlich. 
wenn g vom Grade OÖ ist. Man kann die Lösung 
(6) von (4) auch dahin ausdrücken, daß man sagt, 
die Lösung (3) von (1) gilt auch, wenn statt der 


ganzen Funktion } die gebrochene — gesetzt wird. 


» 

# 
h 
7 
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Die Gleichung 


I(Dy=:z (1’) 
mit (2)=0 für 
für 
hat die Hauptlösung 
für \ 
(3°). 
y=(3) für 
Die Gleichung 
Dy=y(D:z .. (4’) 
mit für 
2()=1 für 
hat die Lösung 
für r<V 
| (6). 
y=(6) für r>0 


Die Hauptlösung von (4) ist 
\ 790) 


(8), 


die von (4’) ist 
y=0 für z< 


y=(8 für r>0 (8 


(6) ist die Heaviside-Lösung,. (8) ist die 
Caspar-Lösung von (4). 

Herr Gauster leitet unter der Annahme (8.381. 
7.9): „Das System sei ursprünglich im Ruhezustand“ 
die Heaviside-Lösung her. Die Annahme des 
Ruhezustandes bedeutet doch die Annahme der 
Anfangsbedingungen 


oder das Verschwinden der Summen (7), also das 
Bestehen der Gleichungen (28) S.371 bei Gauster. 

Herr Gauster bezeichnet die Caspar-Lösung 
als .falsch” (S. 376). Das ist unberechtigt. Als 
Heaviside-Aufgabe muß man (4’) ansehen, als 
zwei partikuläre Lösungen (6’) und (8). Letztere 
ist die Hauptlösung, für erstere sind Anfangs- 
bedingungen nieht gegeben, durch die sie gekenn- 
zeichnet werden könnte. 

Herr Gauster will die Unrichtigekeit der 
Caspar-Lösung daraus herleiten, „daß diese 
Randbedingungen (9) nicht allgemein zutreffen 
können” (S. 376, Z. 19). „Man sieht dies schon 
etwa unmittelbar daraus, daß bei Reihenschaltung 
von LU und R(=w) bei Anlegen des Einheits- 
stoßes die gesamte Spannung vorerst an der 
Selbstinduktion liegt, so daß der 2. Differential- 
quotient der Ladung 9 durch die Gleichung 


bestimmt ist, also z. Zt. £=0 nicht verschwindet.“ 
Herr Gauster greift also auf den einfachsten 
Fall einer einzigen Masche, also auf die Gleichung 
(S. 361) 


1 d d? 

zurück, in welchem n=2 ist, also (9) nur das 
l 
Ve'schwinden von g und für fordert. 
Die Randbedingungen (9) sind also erfüllbar. Aber 
auch für n>2 sind sie im Fall (1’) erfüllbar, da 
hier die Heaviside- Lösung die Hauptlösung ist. 
Das Nichterfülltsein der Randbedingungen (9) tritt 
erst bei dem allgemeinen Fall (4’) ein; bier scheidet 
sich Heaviside- und Caspar- Lösung. 

Berlin. Vahlen. 53la 


Die Verwendung nomographischer 
Methoden zur numerischen Integration 
von Differentialgleichungen. dieser 
Arbeit des Herrn H. Mehlig!') erlaube ich mir, 
folgendes zu bemerken: 


1. Daselbst wird eine Vertafelung des folgenden 
Systems (S) von simultanen Funktionsbeziehungen 
gegeben: 


dp. 
da (a) P,(a) 2m W, 
Pı 
dp, Ps \l d Pi 
(a) Ps m 


Der günstige formale Bau des Systems ermöglicht 
die folgenden weiteren drei Lösungen, deren 
schematische Herleitung vielleicht von einigem 
Interesse sein dürfte: 


a) Unter Benutzung des Gedankens von H. 
Mehlig. Zwischenreehnungen zuzulassen, ergibt 
sich folgende erste Lösung: Eine anfängliche, 
zwei Zwischenrechnungen, einfach kotierte 
Punkte, ein unterbrochener Fluchtlinienzug. 


b) Eine zweite Lösung ergibt folgende Tafel: 
Eine anfängliche, keine Zwischenreehnungen, 
einfach kotierte Punkte ein ununter- 
brochener Fluchtlinienzug. 


e) Die dritte Lösung erfordert schließlich 
keine Nebenreehnungen. hat „Doppelkoten- 
punkte“, ferner ebenfalls einen einzigen un- 
unterbrochenen Fluehtlinienzug. 

Die Lösung a) erfordert die geringsten, die 
Lösung ce) die größten Vorkenntnisse aus der 
Nomographie zu ihrem Entwurf; alle Lösungen 
beruhen jedoch bloß auf durchaus bekannten ge- 
bräuchlichen Verfahren. 

2, Erste Lösung: Wie leichte Umformungen 
zeigen, kann (5) auf folgende allgemeine Gestalt 
gebracht werden: 


+ 
@,=F,+F, ) 
Hierin ist: 

Pa 

Fr, 


9 la) 
F, 


Ps 


Unter Anwendung desselben Verfahrens wie a.a. 0, 
ergibt sich die in Abb. 1 schematisch dargestellte 
Tafel. Wie ersichtlich, ist dieselbe gegenüber der- 
jenigen von H. Mehlig durch größere innere 
Symmetrie ausgezeichnet, auch ist es möglich, ohne 
Zapfenlinien auszukommen. Auch hier können, 


1) Diese Z. 14 (1934), H. 2, S. 122. 
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wie a.a. 0%, die Vorteile der logarithmischen 
Teilungen verwertet werden. 

Die Größe (»,/7,) wird vor Benutzung. die 
(Größen G, und G, während derselben durch 
einfache Nebenreehnung zebildet. Aus dem in der 
Abb. 1 eingetragenen Linienzue 1....17 dürfte 
wohl alles hervorgehen. 


Er, dr, Anz 
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Abb. 3. 


>». Zweite Lösung: Dieselbe geht von folgender 
des Systems (8) aus: 


Hierin ist: 


A,=F,(pı) G, (a) 

B,  F,(p)) 6, (a) 
B, F,(p) 6, (a) 


Das entsprechende Rechenbild zeigt schematisch 
Abb. 2. Es ist vielleicht beachtenswert, daß es 
ebenfalls bloß geradlinige Leitern besitzt. ohne daß 
jedoch, wie a. a. O0. und auch vorstehend. L.og- 
arithmen herangezoren werden müßten. Die Größe 
(Ps/ P,) wird, wie bereitserwähnt. anfänglich gebildet. 
4. Dritte Lösung: Dieselbe ergibt sich auf 
Grund folgender Form des gezebenen Systems: 


A, (Pr. Pa) 


(Q,). 
B, yı (Pı, 22) + B, (9; = 1 | 
Hierin ist: 
A,=F, (a) 

A» F,(a) 
da 

dp, (B3). 
B,=G,(a) 
dp: 


Das entsprechende Rechenbild zeigt Abb. >. 

5. Da die Leitern für die Veränderliche «@ eben- 
falls in allen drei Lösungen bloß auf den gerad- 
linieen Trägern vorkommen. können also bei 
Änderung der Größen f. &7;.&- die Tafeln wie a.a.0,. 
rasch den geänderten Verhältnissen angepaßt 
werden. — Über den praktischen Wert der 
Tafeln kann natürlich erst ihr tatsächlicher Ent- 
wurf entscheiden: in erundsätzlicher Hinsicht 
dürften sie jedoch wohl den Vorzug vor der Lösung 
von H. Mehlig beanspruchen. 


Prag. Alexander Fischer. 488 


Erwiderung. In meiner oben angeführten 
kurzen Mitteilung wies ich auf die Verwendbarkeit 
nomographischer Methoden bei der numerischen 
Integration von Differentialgleicehungen hin und 
erläuterte die Brauchbarkeit dieses Verfahrens an 
einem einfachen Beispiel. Herr Fischer hat sich 
nun die Mühe gemacht. dieses Gleichungssystem 
auch noch mit anderen Methoden der Nomographie 
in erundsätzlicher Form darzustellen, ohne indes 
auf die allein entscheidende praktische Eignung 
dieser Darstellungen einzugehen. 

Da sich meine Ausführungen nicht auf die Form 
der nomographischen Wiedergabe des Gleichungs- 
systems, sondern auf die Erleichterung der zahlen- 
mäßiren Integration von Differentialgleichungen 
unter Zuhiltenahme der Nomographie bezogen, so 
habe ich den Ausführungen des Herrn Fischer 
niehts hinzuzufügen. zumal seine Darstellungen 
sieh allgemein bekannter nomographischer Methoden 
bedienen. 


Hannover. Hans Mehlig. 488a 


Druckfehlerberichtigung. In dem Nachruf 


ke für Prof. Neuendorff ist ein Druckfehler unter- 
B | laufen: das Geburtsjahr des Verstorbenen ist 1877, 
Ap nicht 1873. 
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